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SUMMARY 


A modified V-2 was fired October 26, 1950, at White Sands 
Proving Ground to determine whether, under favorable condi- 
tions, laminar flow at high Reynolds Numbers at supersonic 
speeds was possible. The test area was the surface of a special 
polished 20° included angle conical nose more than 8 ft. long, 
which replaced the standard V-2 warhead. The rocket reached 
M = 3 at 25,000 ft. in approximately horizontal flight before it 
exploded from unknown causes. With the large thermal inertia 
of the 0.18-in. thick skin, the skin temperature lagged well 
behind the boundary-layer recovery temperature. Consequently, 
during most of the supersonic portion of the flight, the heat 
transfer from the air to the surface was greater than that required 
to damp small boundary-layer oscillations according to Lees’ 
theory. Measurements of the nose skin temperature and of 
total pressures in the boundary layer were telemetered in order 
to establish the laminar or turbulent state of the boundary layer. 
The pressure and temperature data indicate that the cone 
boundary layer was laminar at an equivalent flat plate Reynolds 
Number of 5 X 107, which was the largest Reynolds Number 
obtained on the cone. 


Presented at the Aerodynamics Session, Twentieth Annual 
Meeting, I.A.S., New York, January 28-February 1, 1952. 

* The author wishes to acknowledge the assistance of a number 
of people on various phases of this experiment. The study to 
determine the feasibility of a low trajectory for the V-2 was 
made jointly with Dr. A. S. Galbraith, of the Ballistic Research 
Laboratories. F.C. Fay and W. R. Norman, of the General 
Electric Company, performed the trajectory calculations re- 
quired to establish the program cam data. The telemetering 
system and trajectory data were provided by the White Sands 
Annex of the Ballistic Research Laboratories. Dr. T. L. Smith, 
of the Ballistic Research Laboratories, proposed the use of re- 
sistance elements for measuring the skin temperatures and 
Supervised the installation of the temperature and pressure in- 
Strumentation. Miss M. J. Cox carried out many of the numeri- 
cal calculations required. The author also wishes to acknowledge 
the cooperation of the General Electric personnel at White Sands 
under the direction of L. D. White and the White Sands Proving 
Ground Technical Division under the direction of H. L. Karsch. 

t Chief. Supersonic Wind Tunnels Branch. 


A Free-Flight Investigation of the Possibility 
of High Reynolds Number Supersonic 
Laminar Boundary Layers’ 
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Ballistic Research Laboratories, Aberdeen Proving Ground 


NOTATION 


= viscosity 
= heat conduction coeflicient 
Cp = air specific heat at constant pressure 
Pr = Prandtl Number = Cpyu/Ad 


= temperature, °C. 
= air density 


p 

x = distance from leading edge or nose of body 

y = normal distance from cone surface 

l = velocity, ft. per sec. 

Pr = total pressure, Ibs. per sq. in. 

Prs = total pressure just outside of the boundary layer 

M = Mach Number of rocket 

M, = Mach Number just outside boundary layer 

6 = boundary-layer thickness 

0 | = boundary-layer momentum thickness 

Ro, Rz = Reynolds Number based on lengths 0, x 

R.,. min. = Stability limit 

To = stagnation temperature for flow, M, 

Tins, = wall temperature for insulated surface 

r = recovery factor = (Tins. — 7.)/(To — Ts) 

q = heat-transfer rate = cal./cm.? sec. 

h = heat-transfer coefficient = q/(Tin, — Te) = 
cal./em.? sec. °K. 

T = local shearing stress 

C; = local friction coefficient = 7/(1/2)psUs? 

Subscripts 
w = wall conditions 
s = conditions just outside boundary layer 
INTRODUCTION 


HIS PAPER, Lees! prepDicTs that, under certain 
conditions, very large extents of laminar boundary- 
layer flow may be attainable at supersonic speeds. 
According to his calculations, the stabilizing effect 
of heat transfer from the boundary layer to the surface 
should be great enough to damp small oscillation in the 
boundary layer for flow over a flat plate with zero 
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Fic. 1. Critical wall temperature ratios 


longitudinal pressure gradient and constant surface 
temperature. In the light of subsonic work, it would 
be expected that pressure gradients would also have 
important effects on boundary-layer stability at 
supersonic speeds. However, the analyses of Weil® 
and Lees* indicate that pressure gradients only have a 
Above 
M = 38, a favorable pressure gradient seems to have 
little effect, and heat stabilization should be the more 
important. 

The required wall to free-stream temperature ratios 
for damping small disturbances are shown in Fig. 1. 
Lees applied his stability criterion to velocity and tem- 
perature profiles obtained by assuming yu proportional to 
T and Pr = 1. More recently, Van Driest has applied 
the Lees stability criterion to boundary-layer profiles, 
which should be more accurate for air (viscosity varying 
with temperature according to Sutherland’s law and 
Pr = 0.75). Also included in Fig. | 1 is a stability 
limit curve obtained by J. V. Lewis and the author using 
boundary-layer profiles for which Pr = 1 and uy is 
proportional to 7°76 Assuming that Van Driest’s 
curve is more nearly correct,* it is clear that a large 
amount of cooling is required to keep below the stability 


significant influence at low supersonic speeds. 


limit curve, though presumably any amount of cooling 
The R,_ at 


will have some stabilizing effect. 
cr. min. 

which amplification of boundary-layer oscillations can 

begin will be increased and the amplification rate will 


be decreased for R, > R, , since, according to the 


cr. min, 


theory, the amplification rate is 


V1/R, 


Ry 


proportional to 
yr mm However, it appears that the variation of 


with decreasing wall temperature is fairly 


cr, min, 


gradual and that R, 


cy. min, TSS rapidly as the wall 
*In the Readers’ Forum, JOURNAL OF THE AERONAUTICAL 
ScIENCES, February, 1952, C. C. Lin points out certain restric- 
tions in the stability theory which may invalidate these com- 
putations at high Mach Numbers. 


‘if not immediate transition. 
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temperature approaches the stability limit. 


Fig. 2 
shows a typical curve obtained by J. V. Lewis and the 


author (also see reference 4). Therefore, it may pe 
necessary to cool fairly close to the stability limit cree 
in order to obtain spectacular results. 

An experimental check of Lees’ theory would require 
measurements of boundary-layer oscillations similar 
to the work of Schubauer and Skramstad. Presumably, 
instrumentation developments now under way will 
make possible such measurements in supersonic 
boundary layers. However, as an alternative, it is of 
interest to investigate the effects of heat stabilization 
on the location of the boundary-layer transition. If 
large boundary-layer disturbances can be eliminated 
when the heat transfer is large enough, the theory 
predicts the continuance of laminar flow. 

Surface roughness or surface discontinuities may 
constitute large disturbances. Also, it is well known 
that the fluctuations in a turbulent boundary layer 
are large disturbances and result in the phenomenon 
known as transverse contamination and that the “‘con- 
tamination”’ can upset laminar boundary layers at 
very low Reynolds Numbers. Some information is 
available on the effect of shock waves on laminar bound- 
ary layers.» Even weak shock waves produce large 
changes in the boundary-layer profiles, apparently 
always causing boundary-layer separation at the wall 
It is unlikely then that 
laminar flow on the body of a missile could be expected 
to persist downstream of a wing or control surface. 
Even if the wing shock proved inadequate, the turbulent 
wing wake, if contiguous with the body boundary layer, 
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Fic. 2. Variation of Réer, min. with wall temperature ratio 
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should serve as a source of contamination. The 
production of lift on bodies or wings may have a large 
if not dominating effect on the state of the boundary 
laver. For instance, on a body of revolution at in- 
cidence, transition on the windward side is delayed and 
on the leeward side is promoted. Wind-tunnel tests 
at the Ballistic Research Laboratories on a 20° included 
angle cone at J = 3.2 show strong effects even at small 
incidence. At zero angle of attack, transition occurred 
at R, of about 3 X 10° at about the same axial position 
for the top and bottom elements of the cone. At 3° 
angle of attack, transition on the top element occurred 
at R, = | X 10°, while on the bottom element it did 
not occur on the cone and so was at least greater than 
R, = 3 X 10% The photographic evidence available 
suggests that the boundary layer on the bottom surface 
is thinned while the boundary layer on the top is thick- 
ened, probably due to the cross flow of low momentum 
boundary-layer air under the action of the transverse 
pressure gradient. If the laminar boundary-layer 
profile is also significantly altered when the body is at an 
angle of attack, early transition may occur even when 
the heat transfer is favorable. 

Wind-tunnel experiments designed to investigate 
the effect of heat transfer on boundary-layer transition 
have been reported by various investigators.*~> How- 
ever, all of these experiments were limited in that the 
maximum R, was about 5 X 10°. Free-flight Reyn- 
olds Numbers that are likely to be of interest may be 
10 to 100 times greater. Also, in some of these experi- 
ments,’ disturbances of uncertain magnitude were 
introduced in order to have transition occur at a con- 
venient place on the model. In the case of Eber’s 
experiment,® the disturbance of the corner of the cone 
cylinder on the boundary layer is probably significant. 
Wind-tunnel experiments have other uncertainties 
that may or may not be important. Little is known 
about the nature or effect of the air turbulence in the 
test section. Also, it is necessary to test in a region 
free of shock waves originating at surface imperfections 
of the nozzle or test section. 

Free-flight experiments where the heat transfer 
and other factors favor the maintenance of laminar 
boundary layers at high R, appear to be limited in 
number. In an experiment designed to obtain heat- 
transfer data, using a modified V-2 rocket, reported by 
Norris and Fischer,® the wall temperatures on a 45° 
polished cone were below the recovery temperature; 
however, they did not approach the stability limit until 
the R, had decreased to about 3 & 10°. 


PRESENT EXPERIMENT 


(a) Test Surface 

In this paper, a free-flight experiment is described in 
which a V-2 rocket was used to determine whether, 
under favorable conditions, laminar flow was possible 


V-2 TRAJECTORY 


CUTOFF | | 
62 SEC | | 
/| | t M R s0° 
/ | | | | 22 38xI0® 
+ 
Q / NORMAL TRAJECTORY | 37 2 
2 WHITE SANDS | 
co | 493 
40 NEW TRAJ | EXPECTED CUTOFF] 
Yee | | | 


20 40 60 80 00 20 KO. 


RANGE , THOUSANDS OF FEET 
Fic. 3. V-2 No. 61 trajectory. 


up to high R,. The test area was the surface of a 20° 
included angle conical nose more than § ft. long, which 
replaced the standard V-2 warhead. A conical nose 
was chosen because it is believed that the transition 
results on a cone can be related to the theoretical calcu- 
lations on a flat plate. Hantzsche and Wendt" have 
shown that the axial velocity and temperature distri- 
bution for the laminary boundary layer on a cone and 
flat plate should be similar but that the velocity normal 
to the surface in the boundary layer would be different. 
However, the normal velocity is neglected in the sta- 
bility analysis, so that, assuming that disturbances 
equivalent to the two-dimensional disturbances for 
the flat plate exist, the flat plate stability results may be 
applicable. 


(b) Trajectory 


In the normal trajectory used at White Sands for 
V-2 firings, the rocket is launched vertically and gradu- 
ally turns toward the horizontal, so that at the end of 
the propulsion stage it is traveling at about JJ = 5 at 
100,000 ft. Because of the rapid decrease of air density 
with altitude, the Reynolds Number for an 8-ft. length 
would reach a maximum of 4 X 10’ at J/ = 1.4 and then 
decrease to 2.6 X 10®at JJ = 5. On the final leg of the 
trajectory large Reynolds Numbers are attained, but 
the flight conditions are changing rapidly as the rocket 
falls to the ground so that the interpretation of measure- 
ments would be relatively difficult. Also, as previously 
mentioned, it is desirable to keep the angle of attack of 
the nose small, preferably less than 1°. Oscillations 
of the rocket necessarily occur when it re-enters the 
atmosphere, and the oscillations might still be un- 
acceptably large in the lower atmosphere where the 
measurements would be of interest. 

For these reasons, the program control was altered 
to produce the trajectory shown in Fig. 3. For safety 
reasons, the turn toward horizontal could only be 
started after 4 sec. Also the rocket was limited, for 
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structural reasons, to 1g normal acceleration, so that of 
the aerodynamic normal force was not a significant of 
30 factor in turning the rocket. Rather, Fig. 4 shows that Cl 
a jet thrust in the early flight stages and gravity over most , 
Ss of the trajectory were the important factors in turnin 
o ROCKET WEIGHT —— the rocket. JJ = 1 was to be reached at a Pet cet 
flight time with the rocket about 30° from the hori- 
S20 JET THRUST—> zontal, and the rest of the turn to horizontal flight ap 
| was to be produced by gravity with the angle of attack 
programmed to be zero. After burnout, the rocket 
4 GRAVITY ~~ was expected to follow a ballistic trajectory for an ad- - 
4 ditional 30 sec. with roll control continued to impact. 7 
Filo As will be indicated later, after the jet cutoff, the rapid “ 
z decelleration of the rocket with the correspondingly ™ 
T rapid decrease in the boundary-layer recovery tempera- 
2 ture would have caused the heat-transfer conditions to a 
4 \ become unfavorable for the maintenance of laminar i 
—— flow so that the location of the boundary layer transition a 
O vs 10 20 3.0 for a sizable range of heat-transfer conditions would a 
— NUMBER, M have been obtained. Unfortunately, the rocket ex- a 
JET VANE LIFT | | ploded after approximately 50 sec. of flight; therefore, 
Fic. 4. Force normal to trajectory. data were obtained only during the propulsion stage. dr 
The rocket position and velocity throughout flight the 
were measured with the Ballistic Research Laboratories’ the 
Dovap system. These measurements showed that the lay 
n25 predicted course was being followed up to the time of 
ud the explosion. Four static pressure holes were placed the 
4 ° 90° apart around the cone circumference. Differential wa 
w20 TT pressures in the pitch and yaw plane were then meas- the 
= (| ae ee CALCULATION ured, and for JJ > 1.2 the cone angle of attack was cot 
f15 “ computed using Kopal’s values for cones at small yaw. lay 
us | if Fig. 5 shows that for the last 20 sec. of the flight, where dis 
rm 10 : ny M, > 1.25, the cone angle of attack was less than 1°, to 
INDIGATED BY DIFFERENTIAL PRESSURE be small not to bound on 
a 1 9g ary-layer transition. There was no indication from nu 
v 5 ns the rocket motion or the telemetered data that anything cot 
. eer L was amiss, and it has not been possible to establish sig 
O : 10 20 30 the cause of the explosion. cor 
MACH NUMBER, M wit 
Fic. 5. Nose angle of attack. (c) Nose Details 
The first 12 in. of the nose was solid. The rest of the = 
cone was made in four segments that were formed with 
a forming brake to the desired shape, then held in a rig, 
welded together, and finally welded to the solid nose. 
The cone received a final polish with emory paper and 
a fine-steel-wool polishing along the cone elements. 
Photomicrographs were taken of small specimens pre- Ni 
pared this way to show the depth and distribution r 
of the polishing scratches. The specimens were electro- 
plated with a hard chromium plate to preserve the con- 
SURFACE dition of the edge and then were cut perpendicular to the 
lishing direction. Fig. 6, which is a typical photo- 
7 «that the amplitude of the short wave-length surface 
His irregularities is about 60 microin. (see Fig. 6). 
— “a cept at the welds, where bead was left on the inside, = 
Fic. 6. Photomicrograph showing cross section of surface. the skin was fairly uniform in thickness with a tolerance 
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of +5 percent. Fig. 7 shows the important dimensions 
of the conical nose, which was made for us by The 
Glenn L. Martin Company. 

In order to keep the c.g. of the V-2 rocket within ac- 
ceptable limits, it was necessary to put lead ballast in 
the nose. Also, instruments, circuitry, and auxiliary 
apparatus required for the pressure and temperature 
measurements that were made at various positions 
on the nose had to be located within the nose volume. 
Accordingly, an internal structure was fabricated which 
was tied to the base ring supporting the conical nose 
(see Fig. 8). The top of the structure shown contained 
the lead ballast and the bottom held batteries and 
other instrumentation. In order to have access to this 
internal structure, it was necessary to cut a door in the 
conical nose (Fig. 9). It was felt that, even if the door 
were carefully fitted to minimize the surface break, 
there would always be an uncertainty about the im- 
portance of the disturbance there. No attempt was 
made to minimize the surface discontinuity, which 
was about one hundredth of an inch, and a few hun- 
dredths gap was left between the door and the rest of 
the nose to facilitate door removal. It was expected 
that, whatever happened near the nose, the boundary 
layer would be turbulent on the door. 


However, there were two ways the door could affect 
the boundary layer elsewhere on the cone. Any shock 
wave generated at the door would propagate around 
the cone surface and, depending on the Mach Number, 
could affect more or less of the cone surface boundary 
layer. (The curve representing the propagation of a 
disturbance that maintains a fixed angle with respect 
to the cone element at any point is a logarithmic spiral 
on the developed cone.) As will be described later, a 
number of measurements were made near the base of the 
cone directly opposite the door. Assuming that a 
significant disturbance might originate at the front 
corner of the door, we would have some interference 
with the boundary layer for 1/J, < 1.7. For boundary- 


Fic. 8. Internal structure. 


sectionA A 


Fic. 7. Structural details of conical nose. 


Fig. 9. Conical nose. 
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Fic. 10. Skin temperatures for the programmed trajectory. 


layer measurements, a rake of tubes was mounted on 
the surface ahead of the door and the shock wave from 
the rake would interfere for M, < 1.9. Also, trans- 
verse contamination from a turbulent boundary layer 
on the door would spread along the cone surface. 
However, available information indicates this angle 
might be about 10°, whereas it would have to be ap- 
proximately 25° to influence any of the measuring 
stations. In summary then, no known disturbances 
were present for the first 48 in. of the cone; for the 
measuring stations between 48 and 90 in. from the nose, 
disturbances created at the door or by boundary- 
layer tubes just ahead of the door could affect the flow 


'on the smooth side, at least for the most rearward 


stations, for M, < 1.9. It should be noted that the 
cone surface Mach Number was greater than 1.9 for 
the last 13 sec. of the flight. 


CHOICE OF SKIN THICKNESS 


The wall temperatures representing the stability 
limit for the programmed trajectory for the V-2, as 
well as the boundary-layer recovery temperature 
and some calculated skin temperatures, are shown in 
Fig. 10. At the time the skin thickness for this experi- 
ment had to be chosen (about 2 years ago), only Lees’ 
stability curve was available. Since the stability 
results were obviously sensitive to the details of the 
velocity and temperature distributions, it seemed pos- 
sible that the simplifying assumptions of » « 7 and 
Pr = 1 might be important. It was a relatively simple 
matter to obtain velocity and temperature profiles 
using the method of Karman-Tsien and profiles so 
obtained were used to obtain the curve labeled uw « 7°:76, 
Pr = 1. It was not obvious that this curve was closer 
to the truth than Lees’ curve, since there was still the 
Pr = 1 assumption; however, since the curve was 
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more conservative than Lees, it was used as a basis for 
choosing the appropriate skin thickness. 

With the programmed flight conditions, the rocket 
was expected to heat up rapidly, though lagging wey 
behind the boundary-layer recovery temperature. 
The skin temperatures in Fig. 10 showed that, with the 
heat capacity of 0.185-in. thick skin, the skin tempera. 
ture would fall below the limit curves during the acceler. 
ation phase. 


DETERMINATION OF THE STATE OF THE 
BOUNDARY LAYER ON THE CONE 


The identification of the state of the boundary layer 
in this experiment is based on telemetered measure- 
ments of the time rate of heating of the nose and of the 
total pressure distribution close to the cone surface. 
At large flight Reynolds Numbers, the turbulent 
heat-transfer rate is many times larger than the laminar 
heat-transfer rate. By comparing the measured and 
predicted surface temperature rise, the two boundary- 
layer régimes can be differentiated. Also, because of 
the large R,’s, the turbulent boundary layer should be 
much thicker than the laminar boundary layer. Total 
pressure rakes were positioned close to the surface 
so that, based on estimates for the boundary-layer 
thickness at the rake, the total pressure tubes would 
either be within the boundary layer if it was turbulent 
or outside of the boundary layer if it was laminar. 
The total pressure value outside the boundary layer 
can be calculated using the tracking and meteorological 
data for the flight. Also as an experimental check, 
an additional total pressure tube was placed | in. 
from the surface, well outside of the predicted turbulent 
boundary-layer thicknesses. For laminar flow, the 
rake pressures should be the same as, and equal to, the 
total pressure referred to above. For turbulent flow, 
the pressures should decrease going toward the surface, 
and the actual pressures can be compared with approxi- 
mate calculations for the total pressure distribution 
in a turbulent boundary layer at the flight conditions. 
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The use of boundary-layer rakes was necessarily 
restricted because of the resulting flow disturbance. 
Only two stations were used, one 5 in. ahead of the door 
and the other on the opposite side of the cone almost 
at the cone base, 95 in. from the cone tip. Termpera- 
ture measurements were also made at these two stations, 
as well as at a number of other locations on the cone 


surface (see Fig. 11). 
TEMPERATURE MEASUREMENTS 


The temperature elements were resistance grids 
about | by 2 cm. in size, made of 0.001-diameter 
Baleo wire which has a high temperature coefficient 
of resistivity (see Fig. 12). By balancing each grid 
against a fixed Manganin resistor in a battery supplied 
bridge, temperature changes of the grid wire were 
converted to voltages from 0 to 5 volts as required by 
the Hermes telemetering system. There were actually 
three grids or temperature elements arranged side by 
side at most of the stations shown on Fig. 11. In order 
to have sufficient accuracy over the wide range of ex- 
pected temperatures, the resistances of the grids, 
the bridge voltages, and reference resistances were 
chosen so that each of the gages at a given station 
covered a different temperature range as follows: 


A gage 20° to 100° C. for 0-5 volts 
B gage 20° to 270° C. for 0-5 volts 
C gage 20° to 700° C. for 0-5 volts 


Using Sauereisen Cement, the gages were first 
cemented to an !/-in. thick asbestos pad, for conven- 
ience in handling, and then were cemented to the inner 
surface of the skin. Because of the early termination 
of the flight, only one station reached temperatures 
above the low (A) temperature range. Since the 
telemetering accuracy is believed to be 2 or 3 per cent 
of full scale, only the recordings from the two lowest 
ranges (A and B) are useful for this firing. 

Obviously, the failure of a gage to read correctly 
might be due to trouble at any one of a number of 
places in the system between the gage and the film 
record of the telemetered data. Reliability for a given 
gage was established by comparing the initial tempera- 
ture values (at take-off) with the initial temperature 
values of the other gages and by comparing the tempera- 
ture changes indicated by the A and B gages at the 
same temperature station. If the A and B gages at a 
given station did not satisfy these criteria, then neither 
of them was used. In all, there were 31 temperature 
gages and 60 per cent of them gave a satisfactory signal 
on this basis. 

Unfortunately, the temperature gages have a signifi- 
cant effect on the skin temperature being measured. 
The heat input to the gage resistance wire was large 
enough to heat the plate in the neighborhood of the 
gage, at least for the A and B gages (where the heat 
input was approximately 2.9 watts and 1.3 watts, 
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Fic. 12. Typical temperature gage. 


respectively). In the case of the A gage, this heat 
input would be sufficient to heat the plate at a rate of 
0.6° C. per sec. if the heat were restricted to the plate 
over the area of the gage. However, heat conduction 
in the plate greatly reduces the plate temperature rise. 
Laboratory experiments showed that, after voltage 
was applied to an A gage, the plate temperature in- 
creased about 3.5° in the first minute and then at a 
rate of about 1° C. every 3 min. The indicated plate 
temperatures on the V-2 cone were about 10° C. 
above the air temperature at the start of the flight. 
By comparing the total heat input from the gages be- 
fore the start of the flight (power was on about 5 min.) 
with the indicated plate temperature rise at the gages 
and the heat capacity of the nose skin, it was clear that 
there were warm spots in the neighborhood of the gages 
and that most of the skin was unaffected. 

On the other hand, where there was a rapid rise of 
skin temperature during the flight, the asbestos pad 
acted as a heat sink, since the pad has about 30 per cent 
of the plate heat capacity for the same surface area. 
All of this heat capacity is not effective because, as the 
skin temperature rises, the temperature of the asbestos 
pad lags far behind. Calculations have been made 
for the temperature distributions in a pad (one-dimen- 
sional) when the temperature of one surface rises at a 
rate comparable with the experimental rates. After 
the surface temperature has increased about 50° C., 
the heat flow into the pad from the heated surface 
is about 15 per cent of the heat flow into the metal and 
is approximately equal to the heat input from the A 
gage. The computations for the skin temperature rise 
due to aerodynamic heating start with the initial indi- 
cated skin temperature, and therefore they include the 
effect of the heat input from the gages. 

Laboratory tests showed that there were other as 
pects of the performance of the temperature gages 
which had to be taken into account. It was found that 
the variation of wire resistance with temperature, with 
the wire embedded in the Sauereisen Cement, was 
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about 10 per cent less than for the wire in air. Also, 
the layer of cement between the resistance wire and the 
plate surface is a good enough insulator to maintain 
a temperature differential between the wire and the 
plate when the heat generated in the wire flows into the 
plate. When the plate is only heated by the gage, 
almost all of the heat input flows into the plate and the 
wire temperature is about 7° C. higher than the plate 


temperature. However, when the plate is heated at a 


representative rate (with a burner) the temperature 
differential between the wire and the plate decreases, so 
that, when the plate temperature has risen about 70° C., 
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the indicated temperatures are the same. As noted 
before, as the plate heats up, the gage heat input anq 
the heat flow to the asbestos pad approach equality, | 
appears, then, that as the plate temperature rises 4 
rapidly increasing percentage of the heat input from 
the wire flows into the pad rather than into the plate 
so that the wire plate temperature differential js re. 
duced. The indicated temperatures have been cor. 
rected on the basis of these laboratory experiments, 
It should be noted that, at high heat-transfer rates, 
there will be a measurable temperature variation 
through the skin. Estimates indicate that, at the 
highest heat-transfer rates reached during the flight, 
there should have been about a 6° C. differential be. 
tween the inner and outer surface. Possibly then, 
the measured temperatures of the inner skin may 
have been several degrees lower than the mean tempera: 
ture of the skin, which was the quantity calculated. 


BOUNDARY-LAYER MEASUREMENTS 


Fig. 13 shows the details of the boundary-layer 
rakes. Two of the three-tube rakes, each covering 
a different pressure range, were placed side by side at 
Station 95S. The following ranges were obtained on 
Giannini Pressure Transmitters: (1) 0-75 Ibs. per sq.in, 
for 0-5 volts, and (2) 0-300 Ibs. per sq.in. for 0-5 volts. 
The four-tube rake was placed at Station 4S8N just 
ahead of the door. The tube 0.1 in. from the surface 
transmitted pressures from 0-300 Ibs. per sq.in.; 
the other three tubes covered the range 0-75 Ibs. per 
sq.in. Fig. 14 shows how the tube distances from the 
surface compare with various estimates of the bound- 
ary-layer thickness through the supersonic portion of 
the flight. The possible extent of the boundary layer, 
if laminar, is the most important dimension involved, 
and, fortunately, the estimate for the laminar case 
should be fairly reliable. The turbulent boundary- 
layer thicknesses involved more guesswork, and it 
was hoped that the thickness would actually be large 
enough to include all of the total pressure tubes. 


SKIN TEMPERATURE PREDICTIONS 


(a) Laminar Heat Transfer 


The theory of the compressible laminar boundary 
layer has been extensively developed by many authors. 
However, at high Mach Numbers where there are large 
temperature variations through the boundary layer, 
some error may arise from the fact that the Pr is 
assumed to be constant through the boundary layer. 
On the basis of available wind-tunnel data, the recovery 
factor seems to be represented by the theoretical solu- 
tion r = V Pr if the Pr is based on the wall tempera- 
ture. (However, the results are for a fairly limited 
temperature range.) A similar difficulty exists in the 
prediction of the theoretical dependence of _ the 
heat transfer, since the ratio of the local heat transfer 


728 
125 
i git 
sh 
wal 
062 OD. X 0451. TUBING 
est! 
-—40 
8 ma 
Sot 
fac 
eff 
0. 
| 
va 
ift 
| ap 
cas 
an 
lat 
me 
tu 
the 
| 
co 
go 
th 
ha 
be 
tu 
bo 


S noted 
put and 
lity. [ft 
€ rises 
ut from 
the plate 
ial is re- 
een cor. 
Timents, 
er rates, 
ariation 
at the 
€ flight, 
itial be- 
y then, 
in may 
em pera- 
ed, 


‘y-layer 
overing 
side at 
ned on 
r sq.in. 
) volts. 
N just 
surface 
sq.in.; 
IS. per 
the 
ound- 
ion of 
layer, 
olved, 


case 
\dary- 
nd it 

large 


to the local skin friction theoretically should depend 
on the (Pr). There appear to be insufficient data to 
show whether using the Pr in this relation, based on the 
wall temperature, improves or worsens agreement with 
experiment. However, even relatively large errors in 
estimating the laminar heat-transfer rate are not impor- 
tant in this experiment, and the calculations have been 
made as follows: results of reference 11 (including correc- 
tions to cone from flat plate) with Pr based on wall con- 


ditions. 


(o) Turbulent Heat Transfer 


The theory for the compressible turbulent boundary 
layer appears to be in a much more primitive state. 
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Some investigators have predicted that the recovery 

factory should vary as Pr’’ and be independent of 

Mach Number, and others predict a Mach Number 
effect. However, such calculations necessarily involve 

a large number of arbitrary assumptions. Measure- 

ments on a flat plate and cones'*~' and some recent 
unpublished measurements on a cone at the Ballistic 

Research Laboratories give a value for r of about 
0.88, with surface temperatures of about 270° K. 
up to Mach Numbers of 3.5. In view of the un- 
certainty of the theories at this stage, the wind-tunnel 
value may be the most reliable value to be used even 
if the temperatures involved are somewhat different. 

Reynolds analogy for obtaining the heat transfer 
apparently works well for incompressible turbulent 
flow, and because of this it has been hopefully applied 
to compressible flow. Some theoretical justification 
for using the analogy is given by Van Driest’® for the 
case when a turbulent Pr (based on an eddy viscosity 
and eddy heat conductivity) is set equal to one. The 
dependence of the temperature on the velocity through 
the boundary layer is found to be the same for both the 
laminar and the turbulent boundary layers provided 
mean temperatures and mean velocities are used in the 
turbulent case. As yet, there does not seem to be a 
theoretical basis for the dependence of the heat transfer 
on (Pr)**, which dependence applies for incompressible 
turbulent flow and for compressible laminar flow. 

Assuming that Reynolds analogy with a (Pr)’ 
correction can be used, the problem reduces to pre- 
dicting the local friction coefficient C,. Several authors 
predict a sizable decrease of the local friction coefficient 
with Mach Number, defining the coefficient for free- 
stream air properties. Rubesin, et al.,'® have found 
good agreement with the predictions of what is called 
the extended Frankl-Voishel method, and their results 
have been used for the insulated flat plate Mach Num- 
ber effect on Cy. 

By using the von Karman integral method for the 
turbulent boundary layer on a cone, the size of the 
cone boundary layer can be related to the flat plate 
boundary layer. The turbulent profile on the cone 
is assumed to be the same as for the flat plate, and the 


coefficients. 


same law for the dependence of the local shear on Ry is 
assumed to apply. For the case of this V-2, the situa- 
tion is complicated by the fact that the effective start- 
ing point of the turbulent layer is probably well back 
on the cone, though its location is uncertain. This is 
especially true for the flow over the door where the 
boundary-layer momentum losses produced by the 
rake or the leading edge of the door are unknown. 
The correction from a flat plate to a cone when the 
turbulent boundary layer starts from the tip amounts 
to an increase of 17 per cent in the local heat-transfer 
rate.” On the other hand, Van Driest predicts some 
reduction of the local friction coefficient when there 
is a large degree of cooling. Because of the uncertain- 
ties involved, the heat-transfer values used in the com- 
putations have not been adjusted either for the cor- 
rection from the flat plate to the cone or for the effect 
of a large amount of cooling. This source of error 
in the heat-transfer values is probably not more than 
15 per cent. In summary then, the following ex- 
pressions have been used for the turbulent heat-trans- 
fer calculations on the V-2 nose: (a) r = 0.88; (b) 
Reynolds analogy with Pr correction, 


h/Cpp (Pr)** = C,/2 


(C;/Cyy, from page 29, N.A.C.A. T.N. No. 2305), 


Pr evaluated at wall temperature, and 7, temperature 
basis used for other terms. 

It should be noted that, on the basis of theoretical 
results for Pr = 1, the temperature basis for N,, C;, 
and R should be the same. The temperature basis 
for Pr is arbitrary. Fig. 15 shows the comparison be- 
tween the laminar and turbulent values for 4, where 


h = q/AT = q/(Tr — Tw) 
at two stations on the cone. Near the end of the flight, 
(Tr — Tw) for the turbulent boundary layer is about 
80% of (Tr — Tw) for the laminar boundary layer 
because of the higher skin temperatures in the turbulent 
case; therefore, the effect of the spread between the h 
values in Fig. 15 is somewhat lessened. 
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SKIN CONDUCTION AND RADIATION 


Conduction along the skin presumably caused signifi- 
cant temperature changes where large temperature 
gradients existed. The solid nose cone and the sub- 
stantial base ring and associated structure probably 
remained fairly close to the initial temperature through 
the flight. Simplified heat-conduction calculations 
indicated that neither the first station, 6 in. from the 
solid tip, or the last station, 9 in. upstream of the base 
structure, should have been appreciably affected by 
heat-conduction effects. When transition occurred 
on the cone, it is to be expected that surface tempera- 
tures within several inches of the transition region 
will be modified. 

The rocket was fired late in the afternoon, so radia- 
tion from the sun was negligible. During the day the 
rocket nose was shielded from the sun by a sun shade 
(see Fig. 16). 

Calculations show that, at the temperatures reached, 
radiation from the skin was negligible; if the flight had 
lasted longer, radiation would have been included 
(see Fig. 10). 

To avoid radiation from the skin to the internal 
structure, the internal structure was covered with an 
asbestos pad over which was cemented aluminum foil. 


RESULTS 


The pressure and temperature structure of the 
atmosphere obtained using meteorological 
balloons. Runs were made before and after the firing 


both at the launching site and about 12 miles up- 
range. The data obtained showed that atmospheric 
conditions were quite stable during the afternoon of the 
firing. Fig. 17 shows the computed values for the 
Mach Number and the R, (for Station 90) based on the 
trajectory and meteorological data. Using available 
conical flow tables and the above data, the supposed 
static pressures on the surface and total pressures near 
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Fic. 18. Total and static pressure near cone surface. 


acct 
The 
shor 
the 
ture 
for 


F 
for 
fron 
tota 
rans 
tubs 
thar 
anot 
fron 
0.07 
0.1 
pres 
mal! 
are 
g006 
usin 
1,25 
eral 
18S 
Stat 
data 
The 
shou 
coul 
corre 


pera 


at 
\ 2 at 
\ leve 
} \ 
Fig: 
asst 
a re Fic. 16. Rocket at launching site M, 


up- 
dheric 
of the 
r the 
the 
ilable 
posed 
Near 


Reo x 


? 


MC TINULUDS NUINVIBER 


ACE 


SUPERSONIC LAMINAR BOUNDARY LAYERS 


the cone surface were computed. A comparison of the 
computed and telemetered pressures shows agree- 
ment within 2 or 3 per cent which is about the order 
of the telemetering accuracy (see Fig. 18). 


Boundary-layer transition probably occurs well 
forward on the cone during the subsonic portion of the 
flight so that, for a transition R, of say 3 X 10°, when 
M = 0.8, transition would be about 10 in. from the 
cone tip. However, predictions of subsonic transition 
may be very much in error. The rocket had a sizable 
angle of attack up to J/ of approximately 0.8 and 
through the transonic region from M, = 1 to 1.35, 
at which point a = 1°. Also, there is a high sound 
level from the motor which can act on the cone bound- 
ary layer when M <1. Therefore, it seems reason- 
able to believe that laminar flow over a large portion 
of the cone surface would not be realized until A/ > 1. 


Computations were made using a_ step-by-step 
integration procedure, and the predicted skin tempera- 
tures for the two stations where there were both tem- 
perature and pressure instrumentation are shown in 
Figs. 19 and 20. For the calculations, it was arbitrarily 
assumed that the boundary layer is turbulent up to 
M, = 1.25 in all cases. Actually, up to this point, 
the total amount of heat transfer is small so that the 
accuracy of this assumption is not too important. 
The temperature measurements on the door are also 
shown on Fig. 19 for comparison purposes. Since 
there is little difference between the predicted tempera- 
ture variation for Stations 43N and 90S, no prediction 
for Station 65N is shown. 


Figure 19 shows the pressure and temperature data 
for Station 43N. The tubes at 0.07 and 0.14 in. 
from the surface show excellent agreement with the 
total pressure for \/, > 1.25 up to the limit of their 
range (75 Ibs. per sq.in.). The pressures for the 
tubes 0.24 and 0.1 in. from the surface are much lower 
than the total pressure and do not agree with one 
another. It is hard to see why a total pressure 0.24 in. 
from the surface should be much less than pressures 
0.07 and 0.14 in. from the surface or why pressure 
0.1 in. from the surface should be less than the static 
pressure. We believe that something in the system 
malfunctioned for these two pressures and that they 
are not correct. The measured temperatures are in 
good agreement with the predicted wall temperatures 
using laminar heat transfer from the time when 1/, = 
1.25. The temperatures at Station 43N also agree gen- 
erally with the temperature measurements at Stations 
18S and 57S (see Fig. 21). The temperature curve for 
Station 65N on the door is in marked contrast to these 
data and is fairly close to the turbulent predictions. 
The approximate nature of the turbulent predictions 
should be emphasized. The heat-transfer estimates 
could easily be off by 20-30 per cent, which would 
correspond to about 15°~-22° C. in the final wall tem- 
perature. However, even with an error of 50 per cent, 
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Fic. 22. Pressure data at cone base. 


the difference between the laminar and turbulent curyes 
would still be well outside of the estimated experimental} 
error. 

Figs. 20 and 22 show the pressure and temperature 
data for Station 90S. The total pressures for the rake 
on the coarse range, 0-300 Ibs. per sq.in., show a 
systematic separation up to WM, = 2.3, after which they 
coalesce and follow the total pressure curve to the flight 
termination. The pressures from the other rake, range 
0-75 lbs. per sq.in., show the same separation and 
coalescence before the pressure range limit is reached 
(accuracy ~ 1.5 lbs. per sq.in.). 

Assuming a turbulent Pr = 1", the variation of the 
total pressure through a turbulent boundary layer 
for J = 1.43 and 2.17 has been computed, and the 
results are shown in Fig. 23. The velocity was as- 
sumed to depend on one-eighth the power of the dis- 
tance from the wall. 


u/U, = (y/6)"* 


and the experimental wall temperatures were used in 
obtaining the mean temperature variation across the 
boundary layer. The experimental results could be 
compared with these calculations by using predicted 
boundary-layer thicknesses and the tube distances 
from the surface for calculating the experimental 
y/6 values. However, because of the crude nature 
of the boundary-layer thickness estimates previously 
indicated, an alternate procedure has been used. The 
ratio of y/6 for the tube nearest the surface has been 
arbitrarily set equal to the computed curve. The 
boundary-layer thicknesses obtained by this procedure 
are shown on Fig. 14, and they are of the expected 
magnitude. The y/6 values for the other two pressure 
tubes can then be computed, and they are shown on 
Fig. 23. It appears that with the pressure nearest 
the surface (or either of the other two points) arbitrarily 
set on the curve, the other two pressures are reasonably 
consistent with the total pressure variation to be ex- 
pected through a turbulent boundary layer. 
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Fic. 23. Turbulent boundary-layer total pressure distributions. 
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SUPERSONIC LAMINAR BOUNDARY LAYERS 


The measured temperatures at Station 90S clearly 
do not follow the laminar curve and, within the limts 
previously discussed, correspond to the turbulent 
curve. There appears to be a break in the measured 
temperatures when M, = 2.3, although the flight did 
not last long enough to make this break definitive. 
For comparison, a computation assuming laminar 
heat transfer from after 1/, = 2.3 has been included 
on Fig. 20. We can only say that, as far as they go, 
the data are consistent with a change from turbulent 
to laminar heat transfer about M, = 2.3. 

Reference to Figs. 19 and 20 shows that, for the super- 
sonic portion of the flight, the wall temperatures were 
well below the recovery temperatures, after 1J, = 1.8 
falling below the Van Driest stability limit curve. 


CONCLUSIONS 


The possibility of maintaining supersonic laminar 
boundary-layer flow at high R has been investigated 
using the surface of a polished cone more than § ft. 
long as the test area. A modified V-2, with the cone 
in place of the standard V-2 warhead, was fired in a low 
trajectory to reach the desired high Reynolds Numbers 
at supersonic speeds. During most of the supersonic 
portion of the flight, the heat transfer from the air to 
the surface was greater than that required to damp 
small laminar boundary-layer oscillations according to 
the boundary-layer stability theory. 

Telemetered skin temperatures at 43 in. from the 
tip are in good agreement with skin temperature 
predictions for laminar flow. Total pressure readings 
0.07 and 0.14 in. from the surface indicate that the 
boundary layer was less than 0.07 in. thick, several 
times smaller than thickness estimates for a turbulent 
boundary layer. It is concluded that the flow was 
laminar at Station 43N from R, = 20 X 10% at M, = 
1.25 to R, = 40 X 10° at M = 2.7. 

Temperature measurements at Stations 18S and 
57S also follow the lamjnar predictions. 

The measured temperatures at Station 65N are 
considered close enough to the calculated values for 
turbulent flow to indicate that the boundary layer at 
Station 65N was turbulent. The turbulent state of the 
flow is attributed to the disturbances caused by the 
total pressure rake ahead of the door with help from 
the surface discontinuity at the door’s leading edge. 

For Station 90S, on the basis of the total pressure 
measurements with supporting evidence from the skin 
temperature measurements, it is concluded that the 
boundary layer was laminar for the last 6 sec. of the 
flight, during which time the Mach Number increased 
from 2.3 to 2.7 and the Reynolds Number increased 
from 75 X 10° to 90 X 10°. According to theory, the 
laminar boundary layer on a cone at some R, should 
be 0.58 of the corresponding laminar boundary-layer 
thickness on a flat plate. On this basis, laminar flow 


733 


was obtained at 1, = 2.7 at an equivalent flat plate 
R, of 50 X 10°. 
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On the Flow Behind 


SUMMARY 


The flow field behind a curved shock is investigated. The es- 
sential feature of this flow pattern is that the flow is rotational 
with nonuniformity of entropy. The characteristic differential 
equation of the stream function for this flow pattern in both the 
two-dimensional and the axially symmetrical case is derived 
and then linearized by an extended form of the Ackeret iteration 
method. 

The first-order linearized differential equation of the stream 
function, the so-called Sears equation, is solved for the flow 
over a plane ogive behind an attached curved shock. The flow 
patterns are divided into two groups—one is designated as ‘‘weak 
shock,”’ behind which the flow is still supersonic, and the other 
is designated as ‘‘strong shock,’’ behind which the flow is sub- 
sonic. They are investigated separately. Analytic expressions 
for the flow field for both cases are obtained. The velocity field 
can be divided into three components—i.e., one is due to the basic 
disturbances from the surface of the ogive; the second is due to 
the reflected disturbances from the shock front; and the third 
is due to vorticity. Their relative contributions are indicated. 
The first-order linearized theory fails at two points: One is 
that point where the change of stream function due to the change 
of velocity and that due to the change of density just cancel each 
other so that the velocity distribution is indeterminate from the 
stream function. The other point is that at which the reflecting 
disturbance, itself, satisfies the Rankine-Hugoniot relations; 
then, according to the first-order linearized theory, the boundary 
condition at the shock front cannot be satisfied. Both of these 
points occur in the weak shock region. Excluding the imme- 
diate neighborhood of these points, the first-order linearized 
theory shows the general behavior of the flow pattern behind the 
curved shock. Within the neighborhood of these two singular 
points, it is suggested that the exact differential equation should 
be used which can be solved by the method of characteristics 
because the flow behind the attached shock is still supersonic. 


List OF SYMBOLS 


= local sound speed 


a = 
A = (S — So)/Co 

b = factor defined by Eq. (45) 

"4 = constant 

Cy = specific heat at constant pressure 
C, = specific heat at constant volume 
d, = small height 

= arbitrary function 

h = enthalpy 

M = Mach Number 

p = pressure 

q = magnitude of velocity 

R = gas constant 
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an Attached Curved 
Shock’ 


PAIt 


Unwerstty of Maryland 


Ry = reflection factor defined by Eq. (44) 

5 = entropy 

ig = absolute temperature 

U, V = velocity components along x- and y-axis, respectively 

u,v = perturbation flow velocity components along x- and 
y-axis, respectively 

x,y = reference axes 

a = shock angle with respect to U; 

Cp/Cy 

6 = shock angle with respect to Us 

€ = 0 for two-dimensional flow 

€ = 1 for axially symmetrical flow 

n = variable similar to y 

= 

p = density 

o = arc length along the shock front 

y = stream function 

w = vorticity 

Subscripts 


1 = value ahead of shock except as otherwise specified 
2 = value behind the shock except as otherwise specified 
0 = some reference value 


(I) INTRODUCTION 


_ FLOWS INVOLVING a curved shock are of 
great theoretical interest and practical impor- 
tance, our knowledge is limited. One of the essential 
characteristics of flows behind a curved shock is that 
the flow is rotational and the entropy distribution 
varies normal to the stream lines. Furthermore, the 
distribution of entropy is in general not known a 
priort and depends on the solution of the problem. 
This is one of the difficulties in solving the problem of 
flow behind a curved shock. Another difficulty is the 
nonlinearity of the differential equation that governs 
the flow of a compressible fluid. In the present paper, 
an iteration method is suggested to study the whole 
pattern of flow behind the curved shock. 

The differential equation of the stream function in 
rotational flow is linearized by an extended form of the 
Ackeret iteration method.* This linearization is ac- 
complished by assuming the development of the stream 
function y to be of the form 


= y) + vy) + vy) + y) +.-- 
(1) 


where yYo(x, y) is a certain known solution for irrota- 
tional motion corresponding to the problem considered. 
y is the first-order perturbation of Y. For the pur- 
pose of defining the iteration procedure, the function 
¥n+1 is regarded as small compared with the preceding 
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FLOW 


function Y,, and the derivatives have similar relation- 
ships. Then, the total index decides the order of the 
terms; for example, yz is of the same order as y,° or 

In order to solve the linearized equations of y,'s, we 
have to know the entropy distribution in the flow 
which is not known a priori but should be consistent 
with the stream lines of the given problem. Thus the 
equation of y,’s is the so-called “‘se/f-consistent’’ equa- 
tion that often occurs in many problems in theoretical 
physics and which can be solved by a method of iter- 
ation. 

Finally, we investigate the flow over a plane ogive 
behind the attached curved shock by the method de- 
veloped above. The flow field has been analyzed in 
detail in order to determine the relative contribution of 
rotationality and other factors. 

When the flow behind the curved shock is supersonic, 
the method of characteristics can be used to calculate 
the flow pattern numerically. One of the disadvantages 
of an all-numerical method is that it is difficult to ob- 
tain general conclusions. Furthermore, if the flow 
behind the curved shock is subsonic, the method of 
characteristics cannot be used. In the present method, 
both the supersonic and the subsonic flows behind the 
curved shock can be investigated and analytic expres- 
sions can be obtained so that general conclusions may 
be drawn. 

In this paper, the flow is considered as steady, invis- 
cid, and adiabatic. Furthermore, the flow in front of 
the curved shock is assumed to be uniform—.e., isen- 
tropic and irrotational—hence, the whole flow field is 
isoenergetic. 


(II) FUNDAMENTAL EQUATIONS 


Consider two-dimensional or axisymmetric flow with 
velocity components U/ and V along the x and y direc- 
tion. A stream function may be defined in order to 
satisfy the equation of continuity such that 


oy oy 
— = pUy', — = (2) 
oy Ox 

where « = 0 for two-dimensional flow and « = 1 for 


axisymmetric flow. 
For steady isoenergetic flow, the vorticity w for the 
present cases may be written as 
_OV OW pds, (3) 
Ox Oy Rdy 


where .S is the entropy. 


The differential equation for the stream function y 
of the rotational flow! 


72\ 22 IV 32 72\ 22 
a® Ox? Ox Oy a’ } Oy? 


* 


BEHIND AN ATTACHED CURVED 


SHOCK 


If the vorticity w or the variation of entropy with stream 
lines dS/dy vanishes, Eq. (4) reduces to the well-known 
equation of stream function for irrotational flow. 


The equation of state for the adiabatic flow is 
p p” e Cy = fy) (5) 


For isoenergetic flow, if 7) is the stagnation temper- 
ature, we have 


T» Ox oy 


Therefore, together with Eqs. (5) and (6), Eq. (4) is the 
differential equation for the stream function in rota- 
tional flow. It can be investigated if wor dS /dy is spec- 
ified. However, in the case of flow behind a curved 
shock, w or d.S/dy is specified by the shape of the curved 
shock which itself has to be fixed by matching the flow 
calculated by Eq. (4) to that given by the shock-wave 
equation. The problem is thus rather complicated. 
However, we shall show that it can be solved by the 
method of iteration. 


(III) LINEARIZATION OF THE DIFFERENTIATION 
EQUATION OF STREAM FUNCTION 


In order to linearize Eq. (4), we assume that the 
stream function y can be developed into the form of 
Eq. (1). wo(x, y) is a certain known solution for irro- 
tational motion corresponding to the problem consid- 
ered. The simplest form of Yo will be that of uniform 
flow with velocity ly—i.e., 


Wo [polo/(1 + e) te 


Another simple form of Yo will be the solution for flow 
over an unyawed cone which is useful in the case of 
axisymmetric flow. The order of the terms is decided 
by the total index, as explained before. 


In order to find the equations for the y,’s, we have 
to calculate the density ratio in terms of y, and the 
entropy perturbation because, in the present problem, 
the perturbed velocity is of the same order of magni- 
tude as the perturbed entropy. 


Write 
Ovo = pol Ove = —pol 
dy ox (6a) 
go” U.? + Vo, Mo qo ‘do 


The density ratio po/p obtained from the energy equa- 
tion is as follows: 


Po p =1+ (7) 
where 
1— M,* pogo"y* OV (1 — Mo*) pogo*y* 
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which is of the same order of magnitude as y. Ox" ay? Ox Oy 
The vorticity w in iteration form is (: ve ) € 
Substituting Eqs. (1) and (10) into Eq. (4) and with This is the equation of the stream function for irrota- 
the help of Eq. (7), and collecting terms of same order tional flow. Hence, our assumption that y is a soly- 
of magnitude, we have tion for irrotational flow is satisfied. ( 
The first-order differential equation for the perturbed stream function y, is 
Ox? ay? Oy y Oy ay? Ox? 
dy + ay? Ox Oy & dx 
— (Kt + A*| = —poyt w [1 — 1) (12 
a? dy" 0 VoL) ox ai [1 + (¥ ) (12) For | 
whic 
where strai; 
K = (y + 1)Mo?/(1 — Mo?) pogo?y* strea 
= poy (13) 
A* = [2A/(y — 1) (1 — Mo’)] 11 + [(y — 1)/2] MS Fina 
Su 
If the basic flow is a uniform flow of velocity Uy) = constant and I’) = 0, the first-order equation of y, becomes hh ai 
— = —poy’ on [1 + (y — 1)Me? 14 
( Dx? dy? on [ (y (14) 
This equation was first derived by Sears.” This 
For the case of uniform basic flow, the second-order differential equation of perturbed stream function for rota- until 
tional flow is W 
Uo? € Oye U0? ) ov, - over 
( Ox? *> Ov? oy ao? Ox? oy 
— poy'(—for1 + we) [1 + (y — 1) —po wily —1) Mo? med 
9 
Es + KU + as (15) (IV) 
poy”“Go° oy W 
In order to solve Eq. (12), (14), or (15), we have to know the distribution of vorticity or dS/dy in the flow. The pass 
entropy S is a function of the shock angle a of the curved shock with respect to the uniform stream VW, in front of of th 
the curved shock, which is shoc 
Behi 
log | M2 sin? a — + - =<) | (16) unife 
y¥+1 + y+ 1 (y 1)My? sin? knor 
are 
strot 
From Eq. (16) we have for most of the practically interesting problems, the shoe 
, shape of the shock is not known a priori but should be stro! 
dS 1 dSda ; case: 
— = —— — (17) consistent with the stream lines of the given problem— _— 
dy sin 8 pqy' da do e., it depends on the solution of the equation of the men 
where 6 is the shock angle with respect to the flow be- Stream function. This is the so-called ‘‘self-consistent” me 
hind the curved shock and a is the arc length along the equation, which can be solved by the method of iter- and 
shock front. ation as follows: T 
If the shape of the curved shock da/do is given, we First, we assume that the flow is irrotational and ob- ma 
can calculate dS/dy by Eq. (17), and then the equation — tain the solution y, for the given boundary conditions. 8 
of the stream function can be investigated. However, say a 
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FLOW BEHIND AN ATTACHED CURVED SHOCK 


«SHOCK CURVE 


SURFACE OF OGIVE 
v4 


BASIC WEDGE 


SURFACE OF OGIVE 
Fic. 1. Uniform supersonic flow over plane ogive. 


vi = fo(x, (18) 


For the irrotational flow, the original shock is straight, 
which location is known. At the position of the original 
straight shock, we can determine the curvature of the 
stream line behind the shock from Eq. (18), and then 
we obtain da/do along the locus of the straight shock. 
Finally, we can calculate dS/dy by Eq. (17). 

Substituting dS/dy so obtained into the equation of 
y, and solving it for the given boundary conditions 
again, we obtain a new solution, say 


= filx, y) (19) 


This process can be repeated until /, is equal to /,_; or 
until the desired accuracy is obtained. 

We are going to apply this method to study the flow 
over a plane ogive behind an attached curved shock. 
Only the first-order perturbation theory is carried out 
indetail. Extension to the higher order theory may be 
made by a similar procedure. 


(IV) UNIFORM SUPERSONIC FLOW OVER PLANE OGIVE 


We consider a uniform supersonic flow of ideal gas 
passing over a plane ogive. If the leading-edge angle 
of the ogive is smaller than the critical angle, a curved 
shock front is attached to the vertex of the ogive. 
Behind the curved shock, the flow is rotational with non- 
uniform entropy. From the shock polar it is well 
known that, for a given angle of deflection of flow, there 
are two possible configurations of shock angles—the 
strong shock and the weak shock. Behind the weak 
shock, the flow is usually supersonic; while, behind the 
strong shock, the flow is always subsonic. These two 
cases should be treated separately, because the funda- 
mental differential equations for the stream function 
are of different nature—i.e., one is of hyperbolic type 
and the other is of elliptical type. 

The basic flow in this case will be that passing over a 
wedge that will be closely approximated to the ogive 
(Fig. 1). Since the flow in front of the shock is super- 
sonic and the curved shock is attached to the ogive, the 
flow over the upper half of the ogive is independent of 
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that over the lower half. One may tlius consider the 
half plane only. We choose the x-axis along the sur- 
face of the wedge and the y-axis perpendicular to the 
x-axis. Hence, the surface of the ogive will be 
y = d, Y(x) (20) 

where d; is a small quantity of same order of magni- 
tude as y;. V(x) gives the geometrical shape of the 
ogive. 

The location of the curved shock will be approxi- 
mately on the line 


y = x tan do (21) 


The basic flow is now a uniform flow of velocity U», 
and density ps. In the rest of this paper we will use a 
nondimensional form for the stream function, which is 


defined as 

U u 

Oy ple pp 

Oy ox = —(pV 
where U = Ul, + wand V = vare the resultant velocity 
components behind the curved shock. Hence, the rela- 
tions between the perturbed stream function y, and the 
perturbed velocity components u, v and density Ap = 
— pe are 


(22) 


oY, ‘oy (Ap/ pz) + (u (23) 
OY /Ox = —(v/ U2) f 
where 
Ap u S— 
(24) 
pz R 


The boundary conditions of the present problems 
are as follows: 
The condition that the surface of the ogive is a stream 
line gives 
y[x, Y(x)] = 0 (25) 


For the first-order theory, Eq. (25) becomes 
0) = —d, Y(x) (26) 


This is one of the boundary conditions for y(x, y), 
which is to be satisfied at the axis y = 0 but not at the 
profile. 

Let the shock angle with respect to the uniform 
flow in front of the curve be ao + a where a is the shock 
angle of the basic flow with respect to l, and a is the 
first-order perturbed shock angle. 

Under the assumption of small perturbation, it is 
possible to find the relations of a to the perturbed ve- 
locity components w,, v, immediately behind the curved 
shock as follows: 


= Ua-a, v7. = Va-a (27) 


where 
[(pr/p2) — 1] sin (ao + bo) 


Ua = = [ 2S) 
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The 
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Va = U, Vay = 
[(pi/p2) — 1] sin ap sin 69 — cos (ao + | 

(29) 
and U; and p; are the velocity and the density in front 
of the shock, respectively. 69 is the shock angle of the 
basic flow with respect to L’>. 


y¥—1\o 


(Ap/ po); = KeUou, (31) 


Kp = 


Eq. (31) shows the relation between the perturbation 


density and x-wise velocity component immediately 
behind the shock, which is obtained from the energy 
equation. Fig. 2 shows a typical variation of K,[ 
atagiven M, = U;/a. 

The variation of the shock angle a is obtained from 
Eqs. (23), (27), and (31) as follows: 


Us (; l Oy, 
l a 1 + Kl 9”) Oy 


at y = x tan do. 
The first-order vorticity «; is 


(say), where 


ds 27M," sin 2 aol 1)M,? sin? ao + 2] = 4y cot ao[2yM? sin” (¥ 1)] (34 


C, da 


Finally, our problem is to solve the nondimensional 
differential equation for 
‘ 
— 1) = + (y — 1) = 
Ox? Oy? 
wol(x, vy) (85) 


with the boundary conditions 


¥i(x, 0) = —d, Y(x) (36) 
OYi(x, x tan bo) Ua, OW—i(x, x tan 5p) 
Oy Vay Ox 
(37) 


The last boundary condition shows that the perturbed 
velocity components must match those given by Ran- 
kine-Hugoniot relations —i.e., Eqs. (27). At the shock 
line, y = x tan do. 


(V) WEAK ATTACHED SHOCK——-SUPERSONIC FLOW 
BEHIND THE CURVED SHOCK 


When the strength of the basic shock is not large 
enough, the flow behind the shock is still supersonic 
i.e., > 1. Eq. (35) is of hyperbolic type. 

The inhomogeneous partial differential equation 
[Eq. (35)] with the boundary conditions (36) and (37) 
can be solved by the following procedure: Write 


= fo + Wn (38) 
where and are functions such that 


— = O 
Vio(x, 0) = —d, Y(x) 


ra) , x tan by Ua 
oy 
OYw(x, tan bu) 
Ox 
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[(y — sin? ao + 2) sin® ao — (y — 1)] 


where = — | and 


Ox? Oy" 
Yu (x, 0) = 0 (40 
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+ Vi Vay OYio + Vu 
= — K,U.?) 


aty = x tan dp. 
The solution of Eq. (39) is 


vio = —d,[f(x — uy) + f, (x + wy)] (41) 


where f(x — wy) represents the basic disturbance orig- 
inating from the surface Y(x) and f, (x + wy) repre- 
sents the reflected disturbances from the shock front. 
The functions f and f, can be determined from the 
boundary conditions. From the first boundary condi- 


tion of Eq. (39), we have 
Y(x) (42) 


From the second boundary condition of Eq. (39), we 


f(x) + f,(x) = 


have 


f(x) = (R,/b) (43) 


where 
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b = (1 — w tan /(1 + tan bo) (45) 


Fig. 3 shows a typical variation of R, at a given M,. 
Substituting Eq. (43) into Eq. (42), we have 


+ (Ry,/b) f(bx) = Y(x) (46) 


The expression of the basic disturbance function f(x) 
can be determined from Eq. (46). Hence, pw is ob- 


tained. 
In order to solve Eq. (40), we must first know the ex- 


pression w’. Since w,’ depends on the solution of Eq. 
(35), it can only be estimated by the method of iter- 
ation. For first approximation, we estimate w,’ by 
putting yi = Wo; hence, 


U2 l 
«= (7) (47) 
U; Ua, (1 + KpU;?) oy 
at y = x tan dy. From Eqs. (47) and (33), we have 
1f’ — 
wo(x, y) = dy = wo(y) (48) 
dy 
where 
K = 1— R, (1 + (y — 1) 
— 1) Me? (1 + 2?) Uae 
Fle: - (1—R)K,’ (49) 
l 1 te da 


Fig. 4 shows a typical variation of K,,’ at a given J). 
In order to solve Eq. (40), we express w(x, y) in 
terms of Fourier integral of x as follows: 


= anf vy) cos A(E— x) dE (50) 
The solution yy, is then of the following form: 
Vu Vi(y) cos Ax dx + sinAxddy (51) 
0 0 
Substituting Eqs. (48), (50), and (51) into Eq. (40) 
we have 
Y,"(y) + wd? = (52) 
where 7 = 1 or 2, 
Y,” = d*Y,(y)/dy’ 
and 
ww(y) sin (Ay tan do) 
Ai(y) =- 
(53 
wo(¥) cos (Ay tan dp) 
FTo(y) = 
The solution is 


y 
Y, = H,(y) sin \u(y — dn (54) 
An So 


Finally, the solution Wu is 


| 
/ 
| 
| 
§ 
iv 
j 


vu 


y 
dn sin A(x — tan bo) X 
o o 


sin wA(y — n) dd (55) 


Since we consider the region behind the curved shock 
only, Eq. (55) can be integrated with respect to A, and 
we then have 


Yulx, y) = yW(0) W(n) dn (56) 
where 
dW(n)/dn = w10(n) (57) 


It is evident that ¥,(x, 0) = 0. For the second bound- 
ary condition of Eq. (40), we substitute Eqs. (56) and 
(41) into the second boundary condition 


142. 1+R 

(58) 

Substituting Eq. (58) into Eq. (42), we have 

K. | fibx) = Y¥(x) + 


[+ 


R 
+ 
& K.f'(O)x (59) 


Eq. (59) gives the functional relation of the basic dis- 
turbance function f. Since the vorticity function de- 
pends mainly on the conditions at the shock, it turns 
out that the vorticity function and the reflected function 
f, are of similar form and that they can be combined 
together.t After the basic function / is determined, we 


then have 
vi = —d,[f(x — wy) + f(x + + vn 


where both /, and Yj; can be expressed in terms of f with 
the help of Eqs. (58), (56), and (48). 


(60) 


(VI) SrronGc ATTACHED SHOCK—SuBSONIC FLOW 


BEHIND THE CURVED SHOCK 


When the strength of the shock is large, the flow 
behind the shock will be subsonic—i.e., 1/2 <1. Eq. 
(35) is then of the elliptical type. 

By a similar procedure as in the case of weak shock, 
we write 


= Yo + vu (61) 
where Yo and Yu are now functions such that 
+ = 0 


at y = x tan do, where np? = 1 — M,? and 


¢ This simplified relation was first pointed out to the author 
by Dr. F. W. Geiger, of Cornell Aeronautical Laboratory, Inc. 
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Vu (x, 0) = () 
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(63) 


at y = x tan 6p. 

In order to solve Eq. (62), it is convenient to change 
the independent variables x and y into r and @ such 
that 


V x? + 6 = (uy/x) (64) 


Eq. (62) then becomes 


or? Or r2 
Yio(r, O) = —d, Y(r) (65) 
K, 455) = K, 960) 
or r 00 


where 
03, = tan—! tan do) 
K, = [usin 65, + (Um/Va,) (1 + KoU2") cos 65] 
Kg = [((Ua,/Va,) (1 + KeU 2?) sin 05, — cos 05,] (66) 
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bo 


(66) 


BEHIND 


FLOW 


The solution of Eq. (65) is 


—d, >> [A, cos + B, sin (67) 


n=0 


= 


The coefficients A, and B, can be determined from the 
boundary conditions. 
If we develop Y(r) into power series, we have 


Y(r) = > a,r* (68) 
n=0 
The first boundary condition of Eq. (65) gives 
(69) 
The second boundary condition of Eq. (65) gives 
K, sin n05, + K, cos 
B, = a, "= b, (70) 
Ky, cos n65, — A, sin n 05, 
(say). Hence, the expression of Yio is 
dy (x? + pry)” /2 4 
n=0 
uy 
E cos tan-* — ) + b,, sin (» (71) 
x x 


In order to solve Eq. (63), we must estimate the shape 
of the curved shock and the vorticity function w(x, y), 
which is found as follows: 


y) = dwk, n(n — 
n=2 
[—a, sin (w — 1)05, + 5, cos (n — = wo(y) 
(72) 


when x > y cot 69 and w(x, vy) = O when x < ycot 6) and 


we 
= + cot? do = yk, (73) 


(say). 
The solution of Eq. (63) with the boundary condition 
¥u(x, 0) = 0 is simply as follows: 


Vulx, y) = w10(7) dy log 
4 Te 0 n cot 50 


(x — (y+)? 
(74) 

In order to satisfy the boundary condition at the 
shock, we must modify the solution of Yo of Eq. (71) 
so that the last equation of Eqs. (63) is satisfied. We 


= (re — py) — R, + + K| 
2 rot. 


The corresponding formula for irrotational flow is 


(80) 
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find that A, = a, still holds true, but the following ex- 
pression for B, should be used: 


C. 1 
B,, b, nk," = by, (75) 
(say), where 
Ua, 
Rin (1 + Kol 9”) Ro» 
Ovi (y cot do, y) 
Vu y co = g,(y) = y” (76) 


oy n=0 


cot do, vy) 
Ox 


Finally, we have 


= —d, > (a, cos nd + b,, sin nO) + (77) 
n=0 
This is the first approximation to y,. An iteration 
method can be used to obtain more accurate results. 


(VII) VELocITY AND PRESSURE FORMULAS 


A knowledge of the velocity and the pressure distribu- 
tion over the surface of an aerodynamic body at high 
speed is of primary importance in the design of high- 
speed aircraft. In the present problem, it is advisable 
to show the difference between the results for rota- 
tional flow and those for the corresponding irrotational 
flow. 


(a) Velocity Distribution 


Let g be the resultant velocity at any point con- 
sidered. In the region behind the attached curved 


shock, we have 


q= + u)? + = + u + higher order 

terms (78) 
Hence, the distribution of « may be used to represent the 
distribution of the whole velocity field. 

(1) The Weak Shock Case.—From Eqs. (23), (24), 
(35), (48), and (60) we have, for the rotational flow 
behind the attached curved shock, the distribution of 
u as follows: 

+ This formula was first derived by Dr. F. W. Geiger when he 
carried out some numerical examples for the present theory, 
which will be published by Cornell Aeronautical Laboratory, Inc. 


if” [(cot — — 7(0)}) (79) 


The ratio of the velocity distributions of Eqs. (79) and 
(80) may be written as follows: 


Us)ier, = Up + Up +, (say) (81) 


| 
63 
‘hange 
1+ (y¥ IAM, 
(u/U2)ie. = + — Y'(x% — py) 
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The first term 
(up = — wy)/YV'(x — py) 


shows the variation of basic disturbance from the irro- 
tational solution. If the reflection factor R, ap- 
proaches zero and the vorticity factor also approaches 
zero, u, Will approach unity, which represents the case 
when the shock strength is extremely weak. 

The second term 


+ 


(S2 
— py) 


uy = —R,; 
shows the variation of reflected disturbance from the 
shock front. This term depends mainly on the re- 
flection factor R,. A typical curve for R, at a given 
M, is shown in Fig. 3. It shows that, when the shock 
strength is large, this term is extremely important. 
Of course, this analysis breaks down when R,; = o, 
because at this point, according to the first order linear- 
ized theory, the reflecting disturbance itself satisfies 
the Rankine-Hugoniot relation; we cannot use it to 
balance the basic disturbance in order that the resultant 
velocity components satisfy the Rankine-Hugoniot re- 
lations. However, excluding the immediate neighbor- 
hood of R,; = ©, this term will show the general be- 
havior of the reflected disturbance. 

Although the general functional form of the reflected 
disturbance is the same as that of the basic disturb- 
ance, there is a shift in phase between these two dis- 
turbances as is expected. 

The third term gives the variation of velocity due to 
vorticity—i.e., 

K 1— M,? I’ {(cot do — uy] = f'(0) 
(y — DM? V’(x — py) 
(83) 


On the surface of the ogive, y is small in our analysis; 
hence, this term , is generally small there in comparison 
with the other two terms. The main effect of this term 
is to modify the functional form of f. 

(2) The Strong Shock Case.—lIt is possible to derive 
a formula similar to Eq. (79) for the (#/U’2),.4. for the 
strong case. Since there is no corresponding solution 
for irrotational flow in this case, we are not going to 
write the formula down here for comparison. How- 
ever, from the expression for ¥, [Eq. (77)], we see that 
the distribution of u/l’, may also be divided into three 
parts: the basic disturbance that depends on the a,,’s, 
the reflected disturbance that depends on the 6,,’s and 
the vorticity term that is given by Wn. 


(b) Pressure Distribution 


By definition, the pressure coefficient is 


C, = (pb — (84) 
for first-order linearized theory, we have 
= —(2u U2) [20S S1) (5) 


The difference in pressure coefficient between the re. 
sults for the rotational flow and those for the irrota. 
tional flow are twofold. The first is due to the difference 
in velocity distribution, and the second is due to the 
change of entropy (S — S;). 


(VIII) LimiraTIONS OF THE PRESENT THEORY 


From the above analysis, we see that there are two 
singular points at which the present theory breaks 
down. 

The first point is that at which the reflection factor 
R, equals infinity, which has been discussed in See- 
tion (VII). 

The second point is that where Apl.’ = —1. At 
this point the change of stream function immediately 
behind the shock due to the change of x-wise velocity 
and that due to the change of density cancel each 
other. Hence, we cannot determine the vorticity 
terms. 

These two singular points are simply mathematical 
defects of the first-order linearized theory. From nu- 
merical calculations it is found that both of them occur 
when \/, > 1. Hence, it is suggested that, within the 
immediate neighborhood of these two points, the exact 
differential equation for the stream function should be 
used, and, since 1/, > 1, the method of characteristics 
can be used to solve the exact differential equation for 
the stream function numerically. 


Excluding the immediate neighborhood of these two 
singular points, the first-order linearized theory will 
show the general behavior of the flow pattern behind 
the attached curved shock. 
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An Impulse-Momentum Method for Calcu- 
lating Landing-Gear Impact Conditions 


in Unsymmetrical Landings’ 


ROBERT T. YNTEMAT 
Langley Aeronautical Laboratery, N.A.CA. 


ABSTRACT 


A brief review of the landing process is made, and the impor- 
tance of the unsymmetrical landing condition is considered. A 
simple impulse-momentum method that takes account of ap- 
proach conditions and airplane geometric and inertia character- 
istics is given for calculating landing-gear contact conditions for 
successive impacts in unsymmetrical landings. Calcuiated re- 
sults for a large airplane are used to compare the severity of 
landing-gear impacts in unsymmetrical and symmetrical land- 
ings and to show the effects of such factors as sinking speed, 
wing lift, and side drift on impact severity in the unsymmetrical 
case. The importance of landing-gear location and configuration 
is also shown by a comparison of the contact conditions for two 
similar airplanes—one having a conventional tricycle arrange- 
ment of gears and the other having a quadricycle arrangement. 


SYMBOLS 


I, = vertical landing-gear impulse, positive upward 
iF = drag impulse, positive aft 
2 = vertical height of the airplane c.g. above the 
ground, positive up 
6 = angle of pitch measured between the fuselage 
reference axis and the ground plane 
= angle of roll, measured between the spanwise 


¢ = 
reference axis and the ground plane 

z = vertical velocity of any point on the airplane, 
positive upward 

A = pitching velocity of the airplane, positive nose 
up looking forward 

¢ = rolling velocity of the airplane, positive clock- 
wise 

V = vertical velocity for the first gear to impact in 


a landing, positive down (used only in figures, 
corresponds to 

V. = vertical velocity for the second gear to impact 
in a landing, positive down (used only in 


figures, corresponds to Zi,) 


E = energy present for a particular gear at impact 

M = total mass of the airplane 

M, = effective mass acting on a landing gear during 
an impact 

Tea = rolling moment of inertia of the airplane 

Ty, = pitching moment of inertia of the airplane 

Kk, = wing lift factor, ratio of aerodynamic lift force 


to total airplane weight 

Presented at the Structures Session, Nineteenth Annual 
Meeting, 1.A.S., New York, January 29-February 1, 1951. 
Revised and received April 17, 1952. 

*The author is indebted to members of the Dynamic Loads 
Division and particularly to Benjamin Milwitzky and John 
Houbolt, who rendered valuable assistance in the preparation 
of this paper. 

t Aeronautical Research Scientist. 


K, = ratio of side impulse to vertical impulse, defined 
by the average coefficient of lateral skidding 
friction between a tire and the runway 


t = time 
Nr = energy dissipation efficiency of a landing gear 
determines rebound characteristics of the gear 
A;, Bi, Ci, 
D;, Aj, B; = combined constants, moment arms about the 
c.g. 
aj, bj, ¢; = constants that define the location of a landing 


gear relative to the airplane c.g. 


Subscripts 
0 = initial conditions at the beginning of a particular landing 
impulse 
¢ = terminal conditions at the end of a particular impulse, 


represents initial conditions for the free-body portion 
of the analysis 

f = final conditions for the free-body motion, corresponds to 
initial conditions for the next impuise 

g = the center of gravity of the airplane 

1 = the landing gear or gears whose contact initiates a given 
stage of the motion 

j = the landing gear whose contact terminates the free-body 
motion, 7 becomes 7 in the next impulse solution 


INTRODUCTION 


IT RECENT YEARS, the increased size and landing 
speed of airplanes and the relatively large weight 
penalty associated with landing gears have stimulated 
considerable interest in landing problems and in asso- 
ciated design requirements. Although considerable 
effort has been directed toward the solution of some 
of these problems, significant gaps still exist in the 
general knowledge pertaining to the landing phenom- 
enon. The present paper touches briefly on the state 
of knowledge concerning various aspects of the landing 
process and has the particular purpose of evaluating 
the importance of unsymmetrical landings from an 
impact standpoint. The remaining introductory mate- 
rial discusses these aims further and is, for convenience, 
divided into three sections. The first of these sections 
reviews the landing process, the second discusses the 
design requirement for the unsymmetrical landing 
condition, and the third indicates the significance of 
the unsymmetrical landing by means of a simplified 
example. The example furnishes the basis for a more 
complete subsequent analysis permitting consideration 
of additional factors involved in landing. 
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The Landing Process 


A normal landing of an airplane may logically be 
divided into three parts: (1) the approach during 
which airplane motion is largely a function of pilot 
control; (2) the transition from an air-borne to a 
ground-borne condition, during which the vertical 
motion of the airplane is rapidly decelerated and pilot 
control has little immediate effect, particularly in severe 
landings; and (3) the ground run, during which the 
pilot again regains control and maneuvers the airplane 
to a stopping point. 

During the approach, the pilot attempts to make a 
symmetrical landing with a low vertical and _ side 
velocity. Pilot error and unpredictable aerodynamic 
effects such as gusts, and even unusual conditions such 
as forced landings may enter the picture, however, and 
tend to prevent this ideal type of landing. Since these 
situations are next to impossible to predict analytically, 
it appears that limits and ranges for airplane velocities 
and attitudes at initial touchdown are best determined 
from actual situations—that is, from statistical data. 
Unfortunately, only meager data concerning these 
touchdown conditions are currently available. 

During the touchdown or transition period, large 
ground reactions build up rapidly and tend to make 
this period the most critical part of a landing. The 
action of these ground loads on the structure is more 
amenable to analytical treatment and has perhaps 
received the most attention by investigators. Most of 
the investigations, however, have been largely con- 
cerned with some particular aspect of the problem such 
as the determination of a component of the ground 
reaction, consideration of landing-gear behavior, or the 
evaluation of dynamic structural effects on the assump- 
tion that landing-gear forcing functions are known. 
It appears that little consideration has been given the 
tie-in between the actual approach conditions and the 
loads that develop in the structure. As a consequence, 
a rational link between the approach conditions and 
the design impact conditions is missing, particularly 
for the unsymmetrical condition. It is with this link 
that the present paper is largely concerned. 

After the transition period, airplane motions are 
largely a function of pilot action and the loads imposed 
on the airplane structure become those associated with 
braking, turning, and surface irregularities. During 
this period, unstable oscillations such as shimmy may 
also arise. Of the problems associated with the runout 
phase of landing, the shimmy problem has perhaps 
received the most attention, but this attention has not 
always been realistic. There is still much to be desired 
with regard to this and other problems. 


Design Requirements for Unsymmetrical Landings 


Since a unified knowledge of the landing process and 
the associated problems has not been available, current 
design requirements are based largely on experience 


rather than on statistical data or rational analytical] 
treatments. Thus they include hidden factors that 
tend to compensate for the combined effect of many 
conditions not rationally considered in detail. As an 
example, the present requirements specify the same 
impact velocities and landing-gear reactions in yp. 
symmetrical as in symmetrical landings. Experience 
and calculations indicate, however, that landing-gear 
impacts can be appreciably more severe in the un- 
symmetrical case. On the other hand, limited experi- 
mental evidence indicates that the design sinking speed 
is considerably higher than the maximum descent 
velocity encountered in normal landings of many air- 
craft.? In other words, a hidden factor that covers the 
ignorance associated with unsymmetrical landings seems 
to be included in the sinking speed. While the present 
requirements generally permit the design of reasonably 
satisfactory landing gears, at least for aircraft similar 
to those that have given satisfactory service, they may 
be inadequate or perhaps unnecessarily stringent for 
gear design on new types of aircraft. Unfortunately, 
the present lack of data regarding approach conditions 
during landing prohibits the design of landing gears on 
a completely rational basis. 


Consideration of Unsymmetrical Impacts 


The present paper is concerned with the severity of 
initial and subsequent impacts in unsymmetrical land- 
ings as compared with the severity of landing-gear im- 
pacts in symmetrical landings. The treatment con- 
siders as the criterion for the severity of impacts the 
impact velocity and energy rather than the maximum 
loads or the detailed variation of the loads developed 
during the impact. In order to assess rapidly the 
effect that an unsymmetrical landing attitude might 
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LANDING-GEAR IMPACT 


have on the severity of landing-gear impacts, a brief 
study of the simplified two degree of freedom system 
indicated by the airplane sketch in Fig. 1 was made. 
The simplified studies provided the motivations for 
the more complete subsequent analysis, which takes 
into account additional factors involved in landing. 
In the simplified case, the airplane has freedom in roll 
and in vertical translation and is assumed to impact 
on one of the main gears. It is further assumed that 
the landing gears are 100 per cent efficient in dissipating 
energy, that airplane weight is exactly balanced by 
lift, that there are no side or drag loads, and that there 
jsno structural flexibility or other complicating effects. 
The impact velocity and energy for initial and second 
impacts is determined using simple impulse-momentum 
relationships. This approach is used because it re- 
quires no knowledge of the load time history or of the 
maximum load factor developed by the gear. While 
the simplified system does not represent an actual 
landing, it can be used to point out some interesting 
and not commonly considered effects. 

Nondimensional results for the simplified case are 
presented in Fig. 1, where a comparison is made of the 
that landing-gear position can have on the 
severity of landing-gear impacts in unsymmetrical and 
in symmetrical landings. The significant parameter of 
this plot, the abscissa, is the ratio of the outboard 
distance of each landing gear a to the radius of gyration 
of the airplane in roll A,; at the value of one, the 
landing gears are thus located at the radius of gyration 
position. The ordinates are 2/1, the ratio of the 
vertical velocity of the second gear to contact to the 
vertical velocity of the first gear to contact; M,/Mr 
the ratio of the effective mass (the portion of the total 
mass which may be considered to act on a gear during 
an impact) to the total mass and in this case is the same 
for either main gear; and /,/F7 and F2/E7, the ratios 
of the energies absorbed, respectively, by the first and 
second gears to contact to the total initial energy of 
the airplane. The dashed straight lines show results 
for the symmetrical case, while the solid curves repre- 
sent the unsymmetrical case. The curves show a 
number of effects, two of which are particularly sig- 
nificant: (1) The contact velocity for the second im- 
pact in an unsymmetrical landing may be appreciably 
larger than the initial descent velocity of the airplane 
if the landing gears are located outboard of the rolling 
radius of gyration of the airplane—that is, a/K, > 1; 
(2) one of the impacts in an unsymmetrical landing 
must be at least as severe as, and will usually be more 
severe than, each landing-gear impact in a symmetrical 
landing (compare energy curve ABC and line DF). 
The preceding observations are even more significant 
when it is realized that such factors as unbalanced wing 
lift and side drift could result in impacts of much 
greater severity than is indicated by the simplified 
results. Other preliminary studies have indicated that 
unsymmetrical landings may be especially important 
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for certain types of unconventional landing-gear con- 
figurations. 

In the present paper, this simplified impulse-momen- 
tum approach has been broadened and improved to 
permit consideration of additional factors such as 
freedom to pitch, side drift, reduced wing lift or the 
action of gravity, angular velocities and attitude angles 
at initial contact, and landing-gear energy dissipation 
efficiency. In this analysis the airplane is again con- 
sidered a rigid body, since flexibility is expected to 
have little effect on the free-body motions of the air- 
plane and on shock-strut behavior. Adequate con- 
sideration of flexibility effects would, of course, greatly 
complicate the analysis. In the interest of brevity, 
mathematical detail has been kept at a minimum in 
the present paper, although a more detailed develop- 
ment has been made. The method is applicable to 
airplanes with any landing-gear configuration and 
permits determination of contact conditions for suc- 
cessive impacts in unsymmetrical landings in which 
landing gear impacts either do not overlap appreciably 
or occur simultaneously. Unfortunately, the extended 
analysis is not readily adaptable to generalized trend 
studies of the type made for the simplified case. Some 
applications of the method to specific airplanes have 
been made, however, and these results are presented. 


ANALYSIS 


The analysis given in the present paper has two basic 
objectives: (1) the determination of the contact veloc- 
ity for landing-gear impacts that occur subsequent to 
initial touchdown and (2) the determination of the 
effective mass acting at a landing gear during a landing- 
gear impact throughout which the resultant ground 
force does not necessarily maintain a fixed direction. 
For convenience, the analysis is divided into three 
parts entitled Impulse, Free-Body Motion, and Effec- 
tive Mass. In each of these parts, general equations 
applicable to landings in which any combination of 
gears contact the ground simultaneously are first writ- 
ten; these equations are then adapted to the particular 
case of an unsymmetrical landing in which initial 
contact occurs on one gear. 


Impulse 

The impulse analysis considers changes in airplane 
linear and angular velocities resulting from the initial 
landing impact. These changes are determined from 
impulse momentum relationships that assume that the 
time duration of vertical, drag, and side impulses acting 
on a landing gear during an impact is short and, thus, 
that the attitude of the airplane remains substantially 
fixed during the impulse of these ground forces. No 
knowledge of the vertical force time history or of the 
maximum load factor developed by the gear is required 
in the solution of the impulse-momentum equations. 
The vertical impulse is included in these relations as an 
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unknown, the side impulse (associated with a drift or 
yawed landing) is assumed to be related to the vertical 
impulse by some average coefficient of friction, and the 
drag impulse is taken equal to the change in angular 
momentum of the wheel or wheels attached to a landing 
gear as they are spun up to ground speed. 

Immediately following impact, the contacting gear is 
assumed to rebound with a vertical velocity that is a 
function of the impact velocity and the energy dissi- 
pation efficiency of the landing gear. Since no time 
delay is considered in the impulse analysis, the effect 
of reduced wing lift (unbalanced airplane weight) on 
the landirg-gear impulses and therefore on the linear 
and angular velocities of the airplane at rebound must 
be neglected. This error tends to make results of the 
impulse solution somewhat unconservative in cases 
where the wing lift is less than the weight of the air- 
plane. 

General Equations. The pertinent impulse equations 
for an airplane making an unsymmetrical landing in 
which » landing gears contact the runway simultane- 
ously may be written as follows: 

For the vertical motion, 


= — 24) (1) 
For the pitching motion, 


For the rolling motion, 


Taal — = B; + K, I, D; (3) 


where the subscript 7 denotes any of the ” landing gears 
making simultaneous contact; the subscripts 0 and ¢ 
refer, respectively, to conditions at initial contact and 
at rebound; and A,, B,, C,, and D; are appropriate 
moment arms about the airplane center of gravity for 
the landing-gear impulses. Since these three equations 
contain ” + 3 unknowns-—namely, the three velocities 
5, 6,, and ¢, and the » vertical landing-gear impulses— 
they are insufficient in number to permit solution. 
Additional equations may, however, be written for the 
kinematic relationship that must exist between the 
vertical velocity of rebound of each gear and the verti- 
cal, pitching, and rolling velocities of the airplane. 
Thus, for each gear that contacts simultaneously during 
a landing impact, the following relationship must be 
true at rebound: 


Zot + A 6, Sig (4) 


where 4;, is the rebound velocity of each landing gear 7, 
which may be computed from the known initial veloci- 
ties by the following equation: 


by = by V1 — 12, = (&, + — Bio) V1 — 2, (5) 


In this equation n, is the energy dissipation efficiency 


of the landing gear which can be determined from drop 
test data or estimated from the energy dissipation 
characteristics of similar gears. For most gears a valye 
of 80-96 per cent may be used. 

Eqs. (1)-(4) now form a set of (” + 3) equations 
which in any particular case may be solved simul. 
taneously to determine the velocities of the airplane at 
the end of the impulse and, if desired, the vertical 
impulse acting at each of the gears making simultaneoys 
contact. 

One Wheel Landing.—-For the particular case of an 
unsymmetrical landing on any one gear, say m, Eqs, 
(1)—(4) may be reduced to the following set: 


+ A 16, Sms (8) 


Simultaneous solution of Eqs. (6)—(8) yields the vertical 
and angular velocities of the airplane 4,,, 6,, and ¢, im- 
mediately following a landing impact on any gear m. 


Free-Body Motion 


The free-body motion analysis considers changes in 
airplane motion and attitude during the interval be- 
tween rebound and the next impact. Since the impulse 
of the ground forces is assumed to occur instantane- 
ously, the interval actually begins immediately follow- 
ing a landing-gear impact. During this interval the 
airplane is considered a free body under the action of a 
constant gravitational and wing lift force, and the 
translational and angular velocities and the attitude 
of the airplane just before the next impact are deter- 
mined. Aerodynamic moments are neglected in the 
analysis. Since these damping moments generally tend 
to reduce the angular velocities of the airplane, the 
free-body motion analysis is somewhat conservative 
and would be expected to yield somewhat higher 
second-gear contact velocities than actually occur. 
In landings with wing lift less than airplane weight, 
however, the previously mentioned unconservatism of 
results of the impulse analysis tends to balance this 
conservatism. A limited evaluation of some of these 
effects is presented in reference 2. 

The free-body motion equations defining the vertical 
velocity, the height of the airplane center of gravity, 
and the pitching and rolling attitude of the airplane at 
the instant the next impact begins may be written, 
respectively, as follows: 


Sor = Sy, — — Kx) (9 

= — — (g/2)t7(1 — Kx) (10) 
0, = 0 + (11) 

= got (12) 


in which the angular velocities of the airplane (6, and 
¢:) are considered constant during the free-body 
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motion, since aerodynamic moments are neglected. 
In these equations, and the duration of 
free-body motion ¢, are unknown conditions at the 
beginning of the next impact which are to be deter- 
mined. A fifth equation defining the height of the 
airplane center of gravity above the ground in terms 
of airplane geometry and attitude at the instant of 
contact on landing gear j, assumed to be the next gear 
to contact, may also be written as follows: 


a6, + ( 13) 


Zor 


where a;, 6;, and c; are distances defining the location 
of the point of contact of gear j relative to the airplane 
center of gravity. 

Eqs. (9)-(13) may now be solved simultaneously to 
obtain the time duration of the free-body motion, f,, 
in terms of the known conditions at the beginning and 
the computed conditions at the end of the preceding 
impulse. With the value of ¢, obtained, the vertical 
velocity, the attitude, and the height of the airplane 
may be determined from Eqs. (9)—(12). With 2,, now 
known and since 6 and ¢ are assumed not to change 
during the free-body motion, the contact velocity for 
landing gear 7 may be computed by means of the 
following equation: 


Big = Say + A 6, — Big, (14) 

If the computed airplane attitude angles following 
free-body motion indicate that gear 7 was not the next 
to contact, the computation must be repeated for the 
landing gear that does actually contact next. Results 
of the preceding impulse solution, however, usually 
give a good indication of which gear will contact next. 


Effective Mass 

The effective mass analysis is concerned with the 
determination of the equivalent or effective mass acting 
ona landing gear during a landing impact in which three 
independently varying components of ground force act 
on the gear. Formulas for effective mass are derived 
on the basis of equal impulse and energy rather than 
equal instantaneous accelerations as are most conven- 
tional formulas.* This approach does not require the 
assumption that the resultant ground force maintains 
a fixed direction throughout the impact, as is necessary 
in the acceleration approach. Thus, the present method 
provides a convenient means of obtaining a mean value 
for the effective mass acting during an impact which is 
useful in simulating or calculating the severity of the 
impact. 

The vertical impulse acting on the gear, as obtained 
in the impulse solution, may also be expressed in terms 
of the change in vertical velocity of the landing gear 
and a fictitious effective mass 1/,, which can be con- 
sidered to act on the gear—that is, K 
(15) 


In = — 


Zio) 


CONDITIONS 


Solving for .\/,, gives 
= 


(16) 


io) 


(S14 


The results of the impulse solution can be used in 
this equation to obtain the effective mass desired. 
With this effective mass and the impact velocity, the 
kinetic energy of each landing gear at contact may be 
computed by means of well-known formulas for kinetic 
energy. 

In cases where the resultant ground reaction does not 
maintain a fixed direction during a landing-gear impact, 
the effective mass defined by Eq. (16) is an average 
value that satisfies overall impulse and energy require- 
ments for the impact; in cases where the ground re- 
action is fixed in direction, the effective mass obtained 
satisfies instantaneous acceleration requirements as 
well. 

It should be noted that aerodynamic damping tends 
to increase the vertical impulse /, on gear 7 and 
therefore tends to increase the effective mass acting 
on this gear. Since this effect is not included in the 
foregoing impulse analysis, the effective mass deter- 
mined from such results will be somewhat smaller than 
if damping effects had been considered. Flexibility can 
also influence the effective mass, but in general its 
effect on the average value of this mass during an 
impact is probably small. 


PRESENTATION AND DISCUSSION OF RESULTS 


Impact Conditions for Main Landing Gears 


The method outlined in the analysis section has been 
used to calculate values of impact velocity and effective 
mass for the main landing gears of a twin-engined cargo 
airplane in unsymmetrical landings. The outboard 
distance from the airplane centerline to each main 
landing gear is about one and one-tenth times the 
rolling radius of gyration of the airplane, and the 
energy dissipation efficiency of each gear was taken as 
80 per cent. In all of the calculations the airplane 
was assumed to be pitched 3° upward and rolled 7° to 
the left at initial contact on the left main gear, and the 
angular velocities of the airplane at initial touchdown 
were assumed to be zero. These initial conditions re- 
sulted in the second impact occurring on the right 
main gear for each different wing lift and sinking speed 
treated. Impact energies, which include the kinetic 
energy based on the vertical impact velocity and the 
effective mass and the potential energy based on the 
unbalanced effective weight [1/.g(1 — A,)] and an 
assumed travel (strut compression and tire deflection) 
of | ft., are also computed for each landing-gear impact. 
In all cases the potential energy was found to be small 
in comparison with the kinetic energy. 

The impact velocities and energies computed for 
landings in which the side-drift velocity of the airplane 
at initial touchdown was assumed to be zero are shown 
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pacts in unsymmetrical landings. 


12 


8 


fps 


Contact velocities V; and V» for first and second im- 


(First impact on one main 


gear; second impact on opposite main gear. ) 


I6 


Fic. 3. 


4 
bale) 


| st impact 
| 2nd impact 


V, , fps 


Impact energy £ for first and second impacts in un- 


symmetrical landings as a function of initial contact velocity Vi. 
(First impact on one main gear; second impact on opposite main 


gear.) 


on opposite main gear.) 


20; 


/ Symmetrical 


O 4 8 


V,, fps 


Fic. 4. Contact velocities V; and V> for first a1 
pacts in unsymmetrical landings with side drift. 


id second im- 
(First impact 


on one main gear; second impact on opposite main gear.) 


4 
16 (x10 | 


| st impact —- 
2nd impact —-: 


O 
V, fps 


Fic. 5. Impact energy £ for first and second in 


ipacts in un- 


symmetrical landings with side drift as a function of initial con- 
tact velocity V;. (First impact on one main gear; second impact 


for 
the 
init 
the 
for 
velc 
wei 
the 
larg 
geal 
int 
imp 
effe 
low 
T 
at 
syn 
pres 
by 
lanc 
be 
ene! 
ings 
san 
met 
inbe 
the 
mor 
met 
I 
con 
sucl 
effe 
late 
obt: 
the 
drif 
left 
puls 
tent 
that 
the 
bet 
fair] 
tion 
C 
shor 
uns 
ciak 
pac 
T 
mor 
are 


| : 748 

20 or 
2 

=| 

a L / 


cal 


12 


im- 
impact 
) 


LANDING-GEAR IMPACT CONDITIONS 749 


for full-wing lift and for two-thirds full-wing lift in 
Figs. 2 and 3. These wing-lift values correspond, re- 
spectively, to gravitational accelerations of zero and 
/g acting during the impact. The vertical velocity 
for first impacts in the cases considered is, of course, 
the same as the descent velocity of the airplane at 
initial contact. The results presented show that, for 
the airplane under consideration, the vertical velocity 
for second impacts can be greater than the contact 
velocity for the first impact. In cases where all the 
weight of the airplane is balanced by wing lift (A, = 1), 
the increase in velocity for second impacts may be 
largely attributed to the outboard location of the main 
gears relative to the radius of gyration of the airplane 
in roll. Reducing the wing lift tends to increase the 
impact velocity for second impacts still more; these 
effects appear to be more pronounced for landings with 
lower sinking speeds. 

The variation of impact energy with vertical velocity 
at initial contact for first and second impacts in un- 
symmetrical landings is shown in Fig. 3. The curves 
presented show that the energy that must be absorbed 
by the second gear to contact in an imsymmetrical 
landing can be considerably greater than the energy to 
be absorbed by the gear making initial contact. The 
energy computed for each gear when symmetrical land- 
ings were considered was found to be practically the 
same as the energy for the first impact in the unsym- 
metrical case. If the main gears had been located 
inboard of the rolling radius of gyration of the airplane, 
the first impact in an unsymmetrical landing would be 
more severe than each landing-gear impact in a sym- 
metrical landing with the same rate of descent. 

It is of interest to consider what effect a side-drift 
velocity or yawed attitude of the airplane at initial 
contact might have on landing-gear impact conditions 
such as those presented in Figs. 2 and 3. To show this 
eflect, the results shown in Figs. 4 and 5 were calcu- 
lated for initial conditions identical to those used in 
obtaining the results presented in Figs. 2 and 3, with 
the exception that the airplane was assumed to be 
drifting to the right when initial contact occurs on the 
left main gear. In these computations, the side im- 
pulse on the first gear to contact was taken as six- 
tenths times the vertical impulse, and it was assunied 
that the lateral velocity of the airplane was zero when 
the second gear to contact touched down. The 0.6 ratio 
between the side and vertical impulse corresponds to a 
fairly large yaw angle on the order of 15° for conven- 
tional airplane tires (less for high-pressure tires). 

Comparison of the results presented in Figs. 2 and 4 
shows that the vertical velocity for second impacts in 
unsymmetrical landings with side drift can be appre- 
ciably higher than the contact velocity for second im- 
pacts in similar landings without side drift. 

The effect of side drift on impact severity is even 
more striking when the impact energies given in Fig. 5 
are compared with the results given in Fig. 3. Such a 


comparison shows that increases in impact energy as 
large as 50 per cent may result from side draft of the 
airplane just prior to initial touchdown. It is of interest 
to note that in most landings with initial side drift the 
airplane would probably also be drifting laterally when 
the second impact occurs. This condition will result in 
a larger effective mass acting during the second impact, 
and therefore the gear would have to absorb even more 
energy than in the cases considered. 

If the airplane had been drifting to the left instead 
of to the right at initial contact on the left gear, the 
vertical veocity for second impacts would have been 
reduced rather than increased. In such a case, how- 
ever, the effective mass acting during the first impact 
would be considerably larger than in the present case 
or than in a landing without side drift, and therefore 
the impact energy would be larger for the first impact. 
Thus, it is evident that an unsymmetrical landing with 
appreciable side drift will usually result in impacts of 
greater severity than those in an unsymmetrical land- 
ing without side drift or in a symmetrical landing. 


Impact Conditions for Auxiliary Landing Gears 


Results have also been computed for landings of an 
airplane in which the initial attitude was such that 
second impacts occur on an auxiliary gear. Two 
similar airplanes having the same inertia characteris- 
tics and the same location of the main gears are con- 
sidered. One of the airplanes has a conventional tri- 
cycle arrangement of landing gears, while the other 
has a quadricycle arrangement. The attitude of the 
airplanes at initial contact was selected in such a 
manner that first impacts occurred on the left main 
gear and second impacts occurred on the nose gear of 
the tricycle and the right front gear of the quadricycle. 

For the case of a positive pitch angle of 2.5° and a 
roll angle of 9° at initial touchdown, the impact velocity 
and energy for the nose gear of the tricycle arrangement 
are compared with similar results for a forward gear 
of the quadricycle arrangement in Figs. 6 and 7. These 
results show that a forward gear of the quadricycle 
arrangement can be subjected to impacts of far greater 
severity than is the nose gear of a conventional tri- 
cycle arrangement. Each of the forward gears must 
therefore have a larger energy capacity and greater 
strength than the nose gear in a conventional arrange- 
ment. The usefulness of the analysis presented is in 
part suggested by this example. In this case, the 
result suggests that there may be an appreciable 
weight penalty associated with the use of a quadricycle 
arrangement of landing gears. 

The resuits presented in the preceding example are 
for a particular landing attitude that perhaps is not 
too common for large airplanes. Other computations 
for landings in which the initial roll and pitch angles 
were assumed to be considerably smaller show the same 
effects, however. For example, reducing the roll angle 
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from 9° to 5° for a landing of the quadricycle airplane 
in which the descent velocity was 12 ft. per sec. reduces 
the vertical velocity for second impacts only about 
1/, ft. per sec. 

It is to be noted that side-drift or angular velocities 
at initial contact can increase the severity of subse- 
quent impacts for forward gears of the quadricycle, 
while, save for pitching velocities, they can have little 
effect on nose gear impacts of the tricycle arrangement. 
In general, these additional factors tend to widen the 
range of initial attitudes which will result in second 
impacts on a forward gear of the quadricycle arrange- 
ment. It may be noted that the forward gears of the 
quadricycle may also be subject to severe initial impacts. 


CONCLUDING REMARKS 


The present paper presents an impulse momentum 
method that permits the calculation of landing-gear 
contact conditions in unsymmetrical landings if the 
approach conditions at initial contact are known. 
Neglect of flexibility effects is not considered deleterious, 
since in many cases these effects have little influence 
on the free-body motions of the airplane and on landing- 
gear behavior, although they may have a significant 
effect on the internal stresses developed in the gear and 
supporting structure. Neglect of aerodynamic damp- 
ing (particularly damping in roll) in the present analysis 
tends to make the results somewhat conservative, par- 
ticularly in landings with wing lift nearly equal to air- 
plane weight. A more detailed analysis, including 
variable aerodynamic effects and flexibility, would, of 
course, be desirable but is probably not warranted at 
the present time because of the lack of data regarding 
the magnitude of approach conditions likely to be 
encountered in service operations. Even without 


adequate information regarding these approach condi- 
tions, however, the method presented provides a basis 
for assessing the effects of airplane characteristics and 
operating conditions on impact severity and may be 
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Fic. 7. Front-wheel impact energy / for unsymmetrical land- 
ings on one rear gear as a function of the initial contact velocity 
Vi. 


used to evaluate the relative merits of various landing- 
gear locations and configurations from the standpoint 
of impact severity. For the cases that were presented 
as examples, the following effects were noted. 

(a) Landing-gear location is an important factor in 
determining the strength and energy requirements for 
landing gears. 

(b) Contact velocities for second impacts in un- 
syminetrical landings can be appreciably higher than 
the sinking speed of the airplane at initial touchdown. 

(c) The impact energy for either the first or second 
impact in an unsymmetrical landing will usually be 
greater than the energy for each gear in a symmetrical 
landing at the same sinking speed. 

(d) Reducing the wing lift that acts during a land- 
ing tends to increase the contact velocity for second 
impacts in unsymmetrical landings. 

(e) Side drift of the airplane at initial contact can 
result in an appreciable increase in the severity of 
either the initial or second impact in an unsymmetrical 
landing. 

(f) In a quadricycle configuration, each of the for- 
ward gears is likely to be subjected to more severe 
impacts than the single nose gear in a conventional tri- 
cycle arrangement. 
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The Subsonic Calculation of Circulatory 


Spanwise Loadings for Oscillating Airfoils 
by Lifting-Line Techniques 


M. A. DENGLER* ann MARTIN GOLAN Dt 
Midwest Research Institute 


SUMMARY 


A theoretical model is proposed to permit the convenient 
calculation of subsonic circulatory spanwise loadings on oscillat- 
ing finite-span airfoils. The basis for the model is a lifting-line 
construction which represents a fairly direct extension of the 
Weissinger steady-state model to the oscillating case. It is 
shown how the present techniques permit treatment of swept 
wings with the same convenience as is normally encountered 
with straight wings. 

To evaluate the performance of the model, a study is made of 
the spanwise loadings on an upswept rectangular wing of aspect 
ratio 5, oscillating with reduced frequency '/; in six different 
spanwise-symmetrical modes of motion; the results of calcula- 
tions are reported on the basis of the present model and on the 
basis of two prior theoretical treatments. An unswept wing is 
dealt with in the evaluation, since a theoretical background is 
available only for this class of oscillating wings. Calculations 
are also reported for a rectangular wing of- aspect ratio 5 and 
45° sweepback, undergoing the same type motions as the up- 
swept wing. 

The considerations of this paper lead to certain general obser- 
vations regarding the character of the Prandtl, Cicala, and 
Weissinger wing theories. An effort is made to show how the 
present ideas stand in relation to current literature on the subject. 

Finally, recommendations for further study are made, aimed 
toward development of a practical design-office procedure for the 
aerodynamic treatment of oscillating wings of general plan form. 


INTRODUCTION 


Spor PRESENT PAPER CONSTITUTES a revision of an 
original draft, presented at the Twentieth Annual 
Meeting of the I.A.S. and circulated as Sherman M. 
Fairchild Publication Fund Preprint No. 347. It is 
hoped that this new version reflects the many helpful 
comments of the discussers of the earlier paper; in 
particular, the authors are indebted to A. H. Flax and 
H. R. Lawrence, both of Cornell Aeronautical Labora- 
tory, Inc., for their criticisms and suggestions. It is 
also possible for the authors to include in this draft 
certain data and ideas not available at the time of the 
original writing. 

The problem dealt with is that of calculating the 
circulatory loading on an oscillating finite-span airfoil 
in subsonic flow. Since subsonic compressibility cor- 
rections are thoroughly explored in the literature, the 

Based on a paper presented at the Aeroelasticity Session, 
Twentieth Annual Meeting, I.A.S., New York, January 28 
February 1, 1952. 

* Research Engineer. 

t Associate Director for Engineering. 


immediate considerations can be restricted to the in- 
compressible flow case. As its principal objective, the 
present study is aimed toward the derivation of a de- 
sign-office procedure adequate for treating the aero- 
dynamics of oscillating wings with strong sweep. If a 
design-office technique is to be deduced, it is obvious 
that simplicity and calculational convenience must be 
stressed. 

For the unswept airfoil in steady-state flow, the 
classic Prandtl lifting-line theory reduces the airfoil to 
an amazingly simple aerodynamic model. In place of 
dealing with the wing as a lifting surface, requiring the 
solution of a partial integrodifferential equation in both 
the chordwise and spanwise coordinates, Prandtl skil- 
fully approximates the solution in terms of an integro- 
differential equation in the spanwise coordinate only. 
The physical basis for the Prandtl approximation is 
discussed later in this paper, and, from a slightly differ- 
ent point of view, is considered in detail by Reissner,’ 
Lawrence,® and Flax and Lawrence.’ Dealing with the 
oscillating unswept airfoil, Cicala’ showed how the 
Prandtl reasoning could be directly extended to the 
unsteady flow problem. For the solution of practical 
problems and subject to the restriction of dealing with 
airfoils with not-too-small aspect ratios, both the 
Prandtl and Cicala techniques afford enormous com- 
putational advantage over any attempt to deal with 
the problem on a more rigorous lifting surface basis. 

For reasons to be made clear later, Weissinger' 
concluded that the Prandtl approach is inapplicable for 
the treatment of sharply swept wings. However, by 
adopting a new aerodynamic model, he was able to 
retain the essential computational advantages of the 
Prandtl method while meeting many of the important 
requirements of swept-wing theory. (For a critical dis- 
cussion of the Weissinger model, see references 7, 8, and 
9, and see the discussion later in this paper.) Weiss- 
inger’s work is restricted to the steady-state problem, 
and the main objective of the present study is to extend 
his ideas to the oscillating airfoil case. 

It should be particularly noted that all of the theories 
considered in this paper suffer from one major dis- 
advantage: Since the objective of each is to reduce the 
lifting surface problem to a singie integrodifferential in 
one variable, they are forced to concentrate attention 
on the /iff distribution over the span of the airfoil; 
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Fic. lb. Composition of the Fig. la wake. 


it must be expected that their consideration of the air- 
foil pitching moment distribution is on a less precise 
basis. A thorough discussion of this question is given 
by Reissner.’? It should also be noted that the present 
study deals only with circulatory aerodynamic forces; 
finite-span corrections for the apparent mass terms that 
enter the unsteady flow problem are not considered. 


Tue Basic APPROACH USED IN ESTABLISHING THE 
CHARACTERISTICS OF THE THEORETICAL WING MODEL 


In connection with the general methodology for for- 
mulating simplified linear airfoil theories, a highly pro- 
ductive series of papers has appeared in the recent 
literature. For the steady-state problem, this includes 
a contribution by Reissner’ and a set of two papers by 
Lawrence’ and Flax and Lawrence.? The oscillating 
unswept airfoil has been treated by Reissner'? and 
Lawrence and Gerber'! on the basis of certain of the 
ideas of references 7 and 58. 

All of these analyses have in common the approach 
of introducing mathematical simplifications into the 
partial integrodifferential equation of lifting surface 
theory. By introducing the simplifications in a sys- 
tematic, though arbitrary, fashion, the authors men- 
tioned are able to duplicate the mathematics of the 
various known simplified theories and are also able to 
point toward new directions for study. 


Fic. la. A typical horseshoe vorticity pattern. 
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In contrast to the essentially mathematical approach 
of Reissner, Lawrence, et al., the present authors de. 
cided to stress a more physical approach to the formy. 
lation of an aerodynamic model of the wing. Consider. 
ing the Prandtl theory first, the physical reasoning 
leading to the associated steady-state aerodynamic 
model is reviewed from the standpoint of relaxing the 
requirements of the accurate surface loading theory 

‘to the least extent. The Cicala extension for the 
oscillating airfoil is considered in a similar light. 

The Weissinger steady-state theory is then treated 
from an analogous physical point of view, thus leading 
to the considerations important to the objectives of 
this paper—.e., the extension of the Weissinger model 
to the oscillating airfoil case. 

For each of the airfoil theories, a concrete physical 
model is deduced; however, the exact mathematical 
relationship of these models to the integrodifferential 
equation of surface loading theory is not explored. For 
the immediate problem, at least, the authors are of the 
opinion that this line of attack affords greatest sim- 
plicity without sacrifice of theoretical power and au- 
thority. 

One advantage of the present approach is that, once 
the appropriate physical model is deduced, convenient 
means for arranging the required wing calculations at 
once suggest themselves. Also, directions for research 
aimed toward improved theories become evident. 

Since the present reasoning places the Prandtl and 
Weissinger theories in a somewhat different perspective 
from the work of references 7 through 9, this paper 
may also be of value in supplying a different approach 
to classification of these basic wing theories. 


ANALYTICAL TOOLS 


While the approach of this paper is essentially non- 
mathematical, it is desirable to fix certain ideas in 
mathematical terms. To permit this, an outline is now 
given of the aerodynainics of a horseshoe vortex pattern 
with a bound vorticity of uniform intensity across the 
horseshoe span but of harmonically varying intensity 
in the time coordinate. The horseshoe pattern can be 
thought of as a small spanwise segment of the wing in 
question. A schematic representation of the pattern 
is shown in Fig. la. 

To nondimensionalize the analysis, let all distances 
be referred to the local semichord 6 of the wing, and let 
the time be referred to b/V, where V is the wing for- 
ward speed (taken as a constant). Without loss of 
generality, b can be given the value of 1. Let 2n be 


the span of the horseshoe pattern, and let k be the 
reduced frequency of the bound vorticity oscillation. 

If the bound vorticity of the horseshoe leads to 4 
bound circulation, constant along the span, of Te’, 
then known principles permit the deduction that the 
two trailing vortices [, must have the variable inten- 
sity Te"-®, while the wake vorticity has the dis 
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CIRCULATORY SPANWISE 
tribution ikTe““’~*. The positive directions for cir- 
culation and vorticity correspond with the arrow direc- 
tions in Fig. 1. 

It should next be noted that the wake vorticity of the 
Fig. la pattern can obviously be represented as the 
sum of the two wake patterns shown in Fig. lb. Con- 
sidering now Fig. 1b, it is seen that the wake effects 
for the horseshoe can be represented in terms of the 
two-dimensional vorticity distribution y,, plus the 
pattern consisting of the two trailing vortices I, and 
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the two spanwise semi-infinite vorticity distributions ys. 
Evidently, in absolute magnitude, the vorticities 7, 
and are all equal. 

Next, it is of interest to calculate the downward 
induction velocities due to the Fig. lb patterns at a 
point (x, y), the coordinates being defined in the figure. 
By well-known methods and after some algebraic 
manipulation and evaluation of integrals, it can be 
shown that the downwash due to the y; distribution is 
where 


A@, = (1) 


G(kx) = e 


—Ci(kx) — + si(kx)! (2) 


The downwash due to the I’, — y2 combination is given by a more complicated expression and is A@,e""’, where 


ry 
(7 + vy) Ri(x) 


and where 


= + (n + 


! 
= 
—kzr Ue 


2n x x 


tke [[{k(n + y)} R(n (3) 


R.(x) { x? + (n — 
(4) 


+ = “gan 


A(RE) 


For the evaluation of the integral /(kf), it is readily shown that, for ¢ > 0, 
(kt) = F(kt) + H(kg) + — (5) 


where F(kf) is the Cicala function? 


Fag) = + 
£ lee 


and H(k¢) is a new function, readily evaluated, 


MRE) () 


For ¢ < 0, the integral /(kf) becomes 


ACRE) 


(ko) = 2G(kx) — (8) 


It is now of interest to show how the Fig. 1 horseshoe 
pattern can be used to construct the flow field about 
an unswept rectangular wing of span 2/ with a bound 
circulation distribution across the span of ['(y); for this 
simple example, it will further be assumed that the 
wing is oscillating in a spanwise-symmetrical mode 
i.e., [(y) is symmetrical about the mid-span. The 
procedure in this case is well known; rather than think 
of the Fig. 1 pattern as a small spanwise segment, it is 
equally convenient to construct the wing from a series 
of Fig. 1 patterns, each symmetrical about the mid- 
span and each of successively larger span from a value 
near zero to a value approaching 2/. The bound circu- 
lation for each such horseshoe is notated AT. 


Since only symmetrical wing modes are considered, it 
is sufficient to consider the flow over the right-hand side 
of the wing 0<y</. Then, for the horseshoe with 
semispan 7, the bound circulation intensity should be 


Ar = (dT, dy)dy}, 


The downwash velocity at a point (x,y) due to the 
entire wake behind the wing can thus be compounded 
as a spanwise integral representing the summation of 
Aw, and Aw, components. If these summations yield 
downwash velocities of w; and we, respectively, and if 
w, is the downwash due to the bound vorticity, then 
the total downwash at point (x, y) is given by 
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W, = Wy + Wi + We (9) 


If point (x, y) is on the wing and if the downward de- 
flection of the wing at any instant is Z,, then 


w, = (0Z,/07r) + (0Z,/0x) (10) 


The combination of Eqs. (9) and (10) leads to the 
integrodifferential equation governing the aerodynamic 
surface loading problem. 

It will be shown later how the Fig. | type of pattern 
-an be used to construct the flow about a wing of arbi- 
trary plan form. For this purpose, each horseshoe is 
used in place of a small spanwise segment of the 
wing. 


THE PRANDTL AND CICALA MODELS 


The classic Prandtl model for a thin wing in steady- 
state unseparated flow and the Cicala extension to the 
oscillating wing case are based on so-called ‘‘lifting- 
line’ considerations. This designation arises from the 
mental image of a wing with a substantially larger span 
than chord, immersed in a large body of fluid. The 
wing dimensions being small compared with those of 
the fluid boundaries, the lift forces induced on the 
wing appear to act over a relatively thin ‘“‘line’’ in 
space. 

For the purposes of this paper, it is desirable to re- 
place this graphic usage of the term ‘‘lifting line’ by 
its strictly literal connotation—.e., as a line vortex 
immersed in a body of fluid. On the basis of such a literal 
interpretation of the term ‘‘lifting line,’’ it is interesting 
to note that the Prandtl-Cicala model does not entail 
lifting-line considerations but rather can be con- 
structed solely on the basis of surface loading argu- 
ments. This will now be demonstrated, using an un- 
swept rectangular wing purely for illustrative purposes 
and following the line of reasoning advanced by von 
Karman and Burgers.® 

The Prandtl-Cicala simplification of exact wing 
theory presupposes an essentially unswept airfoil of 
appreciable aspect ratio. The vorticity pattern asso- 
ciated with the wing aerodynamics consists of the 
bound vorticity, which runs essentially normal to the 
main flow direction, and the wake vorticity. For the 
steady-state case, the wake vorticity trails rearward 
indefinitely and extends only parallel to the main flow; 
if the airfoil is oscillating, a wake vorticity pattern 
running normal to the main flow direction is added 
(see Fig. 1). For the unswept rectangular wing in 
particular, it has already been shown how the flow 
pattern can be constructed from a series of Fig. 1 
horseshoe patterns, each of slightly larger span than 
the next. 

Consider now the fluid downwash induced at a span- 
wise station of the wing, and consider first the com- 
ponent arising from the bound vorticity. Since the 
wing has a span that is large compared with the chord, 


: 


it can be assumed that each spanwise wing station 
“sees’’ only a slowly changing circulation; thus it be. 
comes possible to approximate the bound vorticity jn- 
duction by setting it equal to that induced on a wing 
in two-dimensional flow, with the saime chord and total 
circulation as the wing station in question. It is im. 
portant to note that by this reasoning the bound 
vorticity at each station acts independently of that at 
all other stations—i.e., interaction of the bound vorticity 
at different span stations is neglected. ; 


In terms of a typical horseshoe pattern of span 2y 
from which the unswept rectangular wing vorticity js 
constructed, the above argument can be rephrased as 
follows: In evaluating the induction effects at a point 
on the airfoil due to the bound vorticity associated 
with each horseshoe pattern, the span is assumed to 
be so large compared with the chord that, in the region 
—n <y <7, the induction of the two-dimensional air- 
foil theory is obtained. Outside the horseshoe span, 
the bound vorticity induction at points on the airfoil 
is zero. 

Now, within the horseshoe span, it is at once clear 
that the combination of the induction due to the bound 
vorticity and the y; distribution of Fig. lb is merely 
equal to the total downwash from the two-dimensional 
theory. The two-dimensional analysis shows that if 
the downwash velocity of an uncambered oscillating 


airfoil is measured at the rear quarter-chord point, 
ik? 


then the two-dimensional downwash A@,je“" is given 
by, for b = 1, 
= (AP/2r)y (11) 


where 


k | 
u(k) = bet (12) 


With regard to the calculation of the induction effects 
due to I’, — y2 within the horseshoe span and due to 
v1 — T, — yz outside the span, it is seen that con- 
sistency with the prior assumptions implies that (a) the 
wake flows can be considered to start at any arbitrary 
chordwise point without introducing appreciable error 
into the analysts and (b) the induction due to the wake is 
relatively constant over the chord and, hence, can be 
evaluated at an arbitrary chordwise point. For con- 
venience, it at once appears reasonable both to start 
the wake and to evaluate the downwash along the line 
x = 0. It is important to note that the origin of x 
need no longer be located as in Fig. 1b; its precise loca- 
tion need not be specified and is immaterial, it only 
being required that the line x = 0 pass somewhere 
within the airfoil chord. 

With the model construction now complete, the 
corresponding differential equation is readily formu- 
lated. Adding up the effects of the various horseshoe 
patterns, 
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— tk[Fik(y + — Fik(y — dn — 
4dr Jo dn 


y 
da J, dy 


It follows that @,, which describes the actual wing 
motion and is a function of the spanwise coordinate y, 
should be evaluated at the rear quarter-chord point. 


The Prandtl-Cicala model has thus been established by 
a reasoning process based entirely on surface loading 
considerations. The terininology ‘‘quasi-surface’’ ap- 
pears to delineate the model features more accurately 
than does the conventional lifting-line description. 


THE WEISSINGER MODEL 


In analyzing the steady-state spanwise loading of 
sharply swept wings, Weissinger' concluded that the 
Prandtl theory is not adequate and will not yield re- 
sults of acceptable accuracy. <A principal reason for 
this is the neglect of bound vorticity interaction, which 
is embodied in the Prandtl approach. Thus, for a 
sharply swept wing, it is clear that a particular wing 
station cannot “‘see’’ a bound vorticity induction that 
can be approximated in terms of a straight two-dimen- 
sional (infinite span) wing. To remedy this deficiency, 
Weissinger proposed a new aerodynamic wing model. 
However, unlike the Prandtl model, the Weissinger 
counterpart cannot be constructed solely from surface 
loading considerations but must be based on replace- 
ment of the wing bound vorticity by an actual lifting 
line i.e., a line vortex. 

The reasoning leading to the Weissinger model is as 
follows: A general conclusion of two-dimensional un- 
steady airfoil theory is that the circulation induced 
around a wing is directly proportional to the down- 
wash acting at the rear quarter-chord point. More- 
over, for the steady-state case, if the bound vorticity 
of the airfoil is assumed to be concentrated into a line 
vortex placed at the forward quarter-chore point, then 
the relationship between circulation and downwash for 
the lifting-line model is the same as that from two- 
dimensional wing theory. Note further that the 
center of pressure of the circulatory forces for a two- 
dimensional wing acts at the forward quarter-chord 
point; for the line vortex model, the total life force acts 
at this same point. 

Since Weissinger presupposes a wing of reasonably 
large aspect ratio, in establishing a model for an airfoil 
with variable circulation across the span, he resorts to 
a modification of the two-dimensional line vortex model 
just described. Thus, he first replaces the bound vor- 
ticity by a line vortex running along the forward 
quarter-chord line of the wing; the circulation of the 
line vortex at each station is the same as that acting 


tk[Fik(n + y)} + — dyn (13) 


for the actual airfoil. The wake vorticity pattern, 
which must start at the line vortex and trail to the rear, 
is then fixed in accordance with the known principles of 
vortex continuity. 

The downwash due to this pattern of line and wake 
vorticity is then calculated along the rear quarter- 
chord line of the wing, thus affording a relationship 
between the downwash at each station and the span- 
wise distribution of bound circulation. The downwash 
so calculated is set equal to that required by the airfoil 
motion along the rear quarter-chord line of the plan 
form. It is readily seen that this procedure leads to an 
integrodifferential equation relating the wing motion 
and the spanwise circulation distribution—.e., the 
Weissinger equation. 

It is to be particularly noted that this model incorporates 
the effects of interaction of the bound vorticity all along the 
wing on the downwash at a particular station, a require- 
ment noted earlier as essential to an adequate swept- 
wing theory. 

The final feature of the Weissinger model is the inter- 
pretation of the computed circulation distribution in 
terms of lift and pitching moment. Since the wing is 
of reasonably large aspect ratio, so that each station 
can be assumed to approximate two-dimensional be- 
havior, Weissinger sets the local lift and pitching 
moment for each spanwise station of the wing equal to 
that which would act on a two-dimensional wing with 
the same local chord and with the same local circula- 
tion as the concentrated line vortex. 


THE PRESENT MODEL 


It is now possible to establish the characteristics of a 
model that extends the Weissinger steady-state argu- 
ment to the oscillating airfoil case. As for the Weiss- 
inger model, the bound vorticity of the airfoil is re- 
placed by a concentrated line vortex running along the 
forward quarter-chord line of the wing. This placement 
of the lifting line coincides with the theoretical location 
of the circulatory center of pressure for a two-dimen- 
sional wing in unsteady motion. 

Consider now the rearward trailing wake vorticity for 
the actual airfoil. This consists of a component parallel 
to the main flow and a component normal to the main 
flow shed from the trailing edge of the wing. Both of 
these components are shown in Fig. | for an isolated 
spanwise seginent of the wing, although it is to be noted 
that the present argument is not restricted to an un- 
swept wing but is applicable for a general plan form. 
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Fic. 2. Real and imaginary parts of the function n(k). 


It has earlier been deduced that, if the bound vorticity 
is i. then the trailing vortices in Fig. 1 have 
the value T, = Fe™"-®, and the wake vorticity 
normal to the flow is given by y,. = ikTe*‘"~®; here, 
£ is a rearward coordinate starting at the trailing 
edge. 

Now define a new rearward coordinate £’ with origin 
at the local forward quarter-chord point. For a unit 
local semichord, it follows that &’ = + (3/2). Then 
the wake trailing vortices for a Fig. 1 type wing element 
ean be written as Te*!-*+°/2)]) and the wake 
vorticity normal to written as 


the flow can be 


Returning now to the construction of the lifting-line 
model, it is an obvious computational convenience to 
assume that the wake springs from the line vortex at 
the forward quarter-chord point. Since the wing is 
supposed to be of reasonable aspect ratio, this type of 
construction seems warranted by the same reasoning as 
was used by Weissinger. Hence, as a preliminary trial, 
the model should consist of the concentrated line vortex 
plus the wake specified in the last paragraph, extend- 
ing from the forward quarter-chord point —’ = 0 to 


As in the Weissinger case, the downwash satisfaction 
should be accomplished along the rear quarter-chord 
line of the wing—.e., the actual wing motion should be 
duplicated by the model along this line. 


Now consider the accuracy of this model in terms of 
a two-dimensional wing; it is clear that one require- 
ment of the model should be that it reduces to the 
exact results for this limiting situation. The two- 
dimensional theory leads to the relations Eqs. (11) and 
(12) between the bound circulation and _ the 
quarter-chord downwash. For the present model, it will 
be found necessary to arbitrarily multiply the bound 
vorticity by the semiempirical factor ne'*/”* in order to 
duplicate these results. ‘This is seen from the following 
reasoning : 


rear 


If the lifting line, actually of intensity Te'*’, is as- 
signed the intensity ne'*/?*Te’, then the downwash 


‘wing circulation. 
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induced by it at the rear quarter chord is, for a unit 
semichord, 


Wy = + (14) 


The two-dimensional wake of the model contributes a 
downwash at this same point of 


0 


where G(k) is given by Eq. (2). 
the exact theoretical value, 


Equating w, + w,; to 


Wit = 27] (16) 
vields the result 
+ ikG(k)} = p(k) (17) 


A plot of vs. k is shown in Fig. 2. It is seen that x 
has a value close to unity for small values of k (<04) 
and departs from unity as & increases. 


Based on this discussion, it is now possible to finally 
fix the details of the lifting-line model for the oscillating 
wing: 


(a) The bound circulation of the actual wing is re- 
placed by a line vortex running along the forward 
quarter-chord line of the wing. For a spanwise wing 
station where the actual bound circulation is I and the 
local semichord is }, the line vortex is arbitrarily as- 
signed the intensity 

(b) The wake intensity at a point behind the wing 
is deduced on the basis of accurate surface loading con- 
siderations—4.e., as if the wake is shed from the trailing 
edge. This determines the wake vorticity at a distance 
After this is fixed, in the 
interests of computational convenience, the wake is 
arbitrarily extended from the trailing edge up to the 
lifting line. 


é’ behind the line vortex. 


(c) The boundary conditions relating to the wing 
downwash are then satisfied along the rear quarter- 
chord line of the wing. 


(d) The model described above leads to an integro- 
differential equation for the spanwise distribution of 
After this is solved, the circulation 
at each station is related to the local wing lift by the 
expression from two-dimensional wing theory. The 
chordwise center of pressure at each station is taken at 
the forward quarter-chord point. 


To fix ideas, the application of this model to the un- 
swept rectangular wing in symmetrical motion is now 
considered. Constructing the wing from a series of 
horseshoe patterns, each symmetrical about the mid- 
span and each of slightly larger span than the next, 
the integrodifferential equation for the problem is 
readily formulated as 
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SAMPLE CALCULATIONS 


To illustrate the character of the results yielded by 
the presently proposed model (\M.R.1I. model), a series 
of calculations are now reported for an uncambered, 
utapered wing of aspect ratio 5, oscillating with re- 
duced frequency '/3; both an unswept wing and a 
wing with 45° sweepback are considered. The two 
wings are assuined to be oscillating in a variety of 
spanwise-symimetrical modes having the following 
shapes: (a) uniform translation of the entire wing, 
normal to its own plane, as a rigid body; ‘b) uniform 
pitching of the wing about the forward quarter-chord 
line as a rigid body; (c) flapping motion of the wing 
according to a mode shape that is linear along the span, 
starting with zero at the root and reaching the value 
one at the tip; (d) pitching motion of the wing about 
the forward quarter-chord line according to a linear 
spanwise mode shape, starting with zero at the root 
and reaching the value one at the tip; (e) flapping 
motion of the wing according to a parabolic spanwise 
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Fic. 3. Bound circulation distribution for uniform vertical 
translation (A.R. = 5.0, k = 1/3) 


mode shape, with amplitude zero at the root and one 
at the tip; and (f) pitching motion of the wing accord- 
ing to a parabolic spanwise mode shape, with amplitude 
zero at the root and one at the tip. 

For the unswept wing, to permit comparison of the 
M.R.I. model results with those from other models, 
calculations are also reported on the basis of the Cicala 
theory and on the basis of the Reissner theory.* * 

For the Cicala and the M.R.I. models, the appro- 
priate integrodifferential equations were solved ap- 
proximately by constructing the flow for each wing 
from five discrete horseshoe vortex patterns, each sym- 
metrical about the wing mid-span. Letting y = / cos @, 
the horseshoe spans were determined by @ = 0°, 20°, 
10°, 60°, and 80°. Satisfaction of the downwash con- 
dition was then accomplished at the stations defined by 
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Fic. 4. Bound circulation distribution for uniform pitching 
(A.R. = 5.0, k = 1/3). 
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Fic. 5. Bound circulation distribution for symmetrical linear 
bending (A.R. = 5.0, k = 1/3). 


6 = 14°, 32°, 51°, 70°, and 90°; save for the 6 = 90° 
station, these stations geometrically bisect the dis- 
tances from the extremity of one horseshoe to the next. 
This construction permits the determination of the 
circulation [ at the five stations where the downwash 
conditions are satisfied. 

For the computations by the Reissner method, the 
procedures outlined in reference 4 were followed without 
change. 

The results for all the calculations are shown in 
Figs. 3 to 8. Considering first the M.R.I. model results 
alone, it is to be noted that the effect of sweep is to 
cause a reduction in the relative inboard wing loading; 
this behavior is similar to that encountered in the 
steady-state loading case. 

For the unswept wing, it is to be noted that both 
the Cicala and Reissner results are in approximate 
agreement, while the M.R.I. model yields generally 
lower wing loading magnitudes, particularly in the tip 
regions. This behavior of the M.R.I. model is not un- 
expected, since it coincides with the general perform- 
ance of the Weissinger method for the steady-state case 
(see reference 1, p. 18). Weissinger points out that his 
lifting-line model is useful mainly in the study of swept 
wings and is, in fact, at a disadvantage when applied 
to unswept cases. 

In view of the paucity of reliable experimental data 
for the spanwise loading of oscillating wings, a more 
detailed analysis of the present calculations does not 
appear to be useful. Perhaps the most efficient means 
for checking the validity of the M.R.I. model under 
present experimental circumstances is by comparing 
the results of carefully conducted flutter experiments 
with those from analytical computations. Limited 
experience with the M.R.I. model in such studies indi- 
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cates that it yields fair results for unswept wings and 
fairly good results for sharply swept wings. 


FURTHER DEVELOPMENTS 


As extensions to the present work, it is worthy of 
mention that a continuation of the present reasoning to 
a more advanced wing model might entail the use of 
two discrete line vortices instead of one to replace the 
wing bound vorticity. Holme! discusses the place- 
ment and use of multiple lifting lines in the steady- 
state case. An extension of his work would be required 
for the oscillating airfoil. It is to be noted that the 
employment of multiple lifting lines enables the study 
of shifts in the local chordwise centers of pressure from 
the forward quarter-chord point. 

Also of interest as an extension to the present work is 
the deduction of convenient computational procedures 
for use with the M.R.I. model. For the steady-state 
case, Diederich'*® shows how the flow about a wing of 
completely arbitrary plan form can be constructed 
from a series of small-span unswept horseshoe patterns 
placed side by side. The use of this technique appears 
well suited for the present purposes and affords sub- 
stantial advantages in reducing the amount of basic 
tabulations required for a_ design-office method. 
Through this approach, wings with nonstraight quarter- 
chord lines can be handled with the same facility as un- 
swept or straight-swept wings. 
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On the Temperature Distribution Along a 
Semi-Infinite Sweat-Cooled Plate’ 


NATHAN NESS? 
Polytechnic Institute of Brooklyn 


ABSTRACT 


A theoretical investigation of the temperature distribution 
along a semi-infinite porous flat plate under the condition of uni- 
form injection from the bottom of the plate is made. 

A heat-balance differential equation of the second order, in- 
cluding a term containing the physical parameters of the plate, 
is derived and is employed in conjunction with the well-known 
Prandtl boundary-layer equations of continuity, momentum, and 
energy for a solution of the problem. 

The temperature distribution along the plate is obtained for 
the respective cases of thermal conductivity not equal to, and 
equal to, zero. Results show that the inclusion of the thermal 
conductivity term in the heat-balance equation eliminates the 
infinite temperature gradient at the leading edge. 

Other results when thermal conductivity is considered are: 
(a) The boundary conditions at the leading edge—i.e., the ratio 
of wall to free-stream temperature and the axial temperature 
gradient there—may not be chosen arbitrarily but are related 
mathematically to each other; (b) the ratio of thermal to velocity 
boundary-layer thickness is less than unity at the leading edge for 
Prandtl Numbers greater than 0.35 and is equal to unity at an 
infinite distance downstream; and (c) at an infinite distance 
downstream, the wall temperature approaches the coolant tem- 
perature. 


NOMENCLATURE 


x = distance along the plate measured from the leading edge 
= distance normal to the plate 

u = velocity component parallel to plate 

v = velocity component normal to plate 

p = mass fluid density 


= coefficient of viscosity 
v = p/p, kinematic viscosity 
Cp = coefficient of specific heat at constant pressure 
P = wc,/k,, Prandt] Number 
1 = distance between coolant reservoir and porous plate 
h = thickness of porous plate 
k, = coefficient of heat conduction of the hot fluid 
kn = coefficient of heat conduction of the plate 
T = temperature, deg. abs. 
T. = temperature of the coolant 


a = (he) f T dy, a constant limited to 0 < a < 1 for in- 


jection cooling 
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r = constant ratio of porous to total plate cross-sectiongl 
area 

6 = boundary-layer thickness 

A,» = vw5,/(6v), nondimensional velocity boundary-layer thick. 
ness 

= Pv. 6./(6v), nondimensional thermal boundary-layer 
thickness 

= nondimensional length along the plate 

7 = y/6, nondimensional distance normal to the plate 

L = finite length of plate 

K = — + 7], a constant 

C = —1/(210K), a constant 

= 1 — 7,/T, 

J = &/6, = X/(PX,), boundary-layer thickness ratio 

F is defined by Eq. (13). Prime =: d/dt; double prime = 
d?/dg?. 
Subscripts 

w = hot surface of the plate 

O = leading edge of plate 

o = infinite distance downstream along the plate 

. v = velocity boundary layer 


t = thermal boundary layer 
1 local values at the outer edge of the boundary layer 


INTRODUCTION 


8 tee PROBLEM OF DETERMINING the temperature dis- 
tribution along a semi-infinite porous flat plate, 
subjected to a high-temperature gas flow on one side 
and cooled by means of cold fluid injection from the 
other side of the wall, is considered. This method of 
decreasing wall temperatures is known as ‘‘sweat- 
cooling.” 

The fluid dynamic equations considered are those 
that are valid in the velocity and thermal boundary 
layers. These equations are the well-known Prandtl 
boundary-layer equations of continuity, momentum, 
and energy. The wall temperature is introduced as a 
dependent variable through the derivation of a heat- 
balance differential equation. This equation is of the 
second order; the second-order term (which is multi- 
plied by a factor containing the physical parameters of 
the plate) arises from a consideration of the thermal 
conductivity along the plate and appears not to have 
been treated previously in the literature. 

The assumptions made in this investigation are: 
(1) The fluid is incompressible; (2) the flow is entirely 
laminar; and (3) the hot fluid flowing a'ong the wall 
and the coolant flowing through the pores are the 
same. 

In the solution of the laminar boundary-layer equa- 
tions, the momentum and energy equations are reduced 
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to the von K4rm4n integral form. Velocity and tem- 
perature boundary-layer profiles as polynomials of the 
fourth degree are assumed. The assumption of in- 
compressibility is found to imply that the momentum 
equation is not coupled with the energy equation; there- 
fore, the substitution of the velocity boundary-layer 
profile into the momentum integral equation permits 
the determination of the velocity boundary-layer thick- 
ness without the use of the energy or heat-balance equa- 
tions. 

The temperature distribution along the plate and the 
thermal boundary-layer thickness are calculated, first 
with the thermal conductivity term included in the 
heat-balance equation and then with the thermal con- 
ductivity term neglected. For the first case, the tem- 
perature distribution and the thermal boundary-layer 
thickness are obtained by the numerical solution of two 
simultaneous first-order nonlinear ordinary differential 
equations, which result from a combination of the heat- 
balance and energy equations, and from the energy 
equation itself. For the second case, a combination of 
the simplified heat-balance equation (i.e., the heat- 
balance equation with plate thermal conductivity 
neglected) and the energy equation results in an alge- 
braic equation that is solved graphically for the thermal 
boundary-layer thickness. The substitution, then, of 
the thermal boundary-layer thickness into the simpli- 
fied heat-balance equation gives the wall tempera- 
ture. 


DERIVATION OF HEAT BALANCE EQUATION 


Basic EQUATIONS 


The laminar incompressible two-dimensional bound- 
ary-layer equations for flow over a flat plate parallel 
to the stream are: 

The equation of conservation of momentum: 


u(Ou/Ox) + v(Ou/Oy) = (1) 
The equation of the conservation of energy: 
u(OT/dx) + v(OT/Oy) = (v/P)(0*T/dy?) (2) 
The equation of continuity: 
(Ou/Ox) + = 0 (3) 


It is assumed here that the viscosity and thermal coef- 
ficients are constant; moreover, the term representing 
the viscous dissipation—namely, (v/c,)(Ou/Oy)*—has 
been omitted because of the assumption of incompressi- 
bility. 

The wall temperature is introduced as a variable by 
considering the heat balance of an element [of length 
Ax, of unit width, and of height (h + /)] of the system 
in Fig. 1. In the region between the porous plate and 
the coolant reservoir, it is assumed that O77/Ox is 
negligible and that the coolant is in one-dimensional 
flow normal to the plate. Furthermore, experimental 
results have shown! that the difference between the 
temperature of the cooling medium entering the pores 
and the temperature of the wall decreases to a negligible 
amount within an extremely thin layer inside the wall 
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The resulting heat-balance equation is, then, 


k(OT/Ov)w + ah[(1 — r)km + = 
(p0) — Te) (4) 


where a = 
h 


a <1 for injection cooling; + is the constant ratio of 
porous to total plate cross-sectional area; and the sub- 
script w refers to conditions at the hot surface of the 


Wnts) fi T dy, a constant limited to 0 < 
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compressibility has been made in the derivation of Eq. 
(4), although k, and c, have been assumed con. 
stant. 

The first term of Eq. (4) represents the heat trans. 
ferred from the hot fluid to the plate; the second term, 
which contains the physical parameters of the plate, js 
the heat transmitted along the plate; while the term on 
the right side of the equation is the heat absorbed by 
the coolant in its flow from the coolant reservoir to the 


plate. It should be noted that no assumption of in- hot surface of the plate. 


SOLUTION OF THE Basic EQUATIONS 


A solution of Eqs. (1) 6 (4) is obtained by the Karman-Pohlhausen method of reducing the momentum and 
energy equations to integral form and assuming velocity and temperature profiles as polynomials in y of the fourth 
degree. 

Integrating Eq. (1) between the limits y = 0 and y = 4,, utilizing the equation of continuity and satisfying the 
following boundary conditions: at y = 0, 

u=0, v =v, (a constant) 


at y = 6, 


u= v = v(6,), Ou/Oy = 0 
one obtains 
bv 5p 
u? dy — u dy + = —v(OUu/OY) (5) 
0 


A polynomial of the fourth degree as an approximation to the velocity boundary-layer profile satisfying the fol- 


lowing boundary conditions: at y = 0, 


u = 0,v = vp, 


at y = 6, 


yields u OUu/Oy = O'u/dy? = 0 
u 2 1+ 3), 
uy 1 +A, 1 +X, 


where 7, = ¥/5r; Ay = Vwd,/(6v); and 6, and \, are the dimensional and nondimensional velocity boundary-layer 
thicknesses, respectively 

The substitution of Eq. (6) into Eq. (5) results in 

ie + 108A,2 + 119A, + 37) 

35(1 + A,)?(3A,2 + 3A, + 1) 

where = (v,./%1)?(%x/v) is a nondimensional length. Direct integration of Eq. (7), with the use of the initial 
condition & = 0, A, = 0, gives A, implicitly as a function of é. 


E = 
— T 
2 


dry (7) 


(1 + X,)? 
(3A,2 + Br» + 1) 


The integral form of the energy equation is obtained by integrating Eq. (2), which contains the velocity compo- 
nent u and the temperature 7, across the greater (i.e., velocity or thermal) boundary-layer thickness. An assump- 
tion that the thermal boundary-layer thickness is greater than the velocity boundary-layer thickness, for the 
Prandtl Numbers considered in this paper, can be shown to lead to contradictory results; therefore, Eq. (2) is inte- 
grated between the limits y = 0 and y = 6,._ Use of the equation of continuity and the following boundary condi- 
tions: at y = 0, 


35V3 
tan-! V3(2, + 1) + tan? V3 + 10 (8) 


T= T(x) 


Uv = Vy; 
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Ba. 
led con- 
v = T= 7%, OT/dy = 0 
at trans- 
yields 
pnd term, | 
‘Plate, is f uT dy — udy + — Te) = —(v/P)(OT/Oy) (9) 
term on 0 0 
rrbed by A polynomial of the fourth degree as an approximation to the thermal boundary-layer profile satisfying the fol- 
ir to the lowing boundary conditions: at y = 0, 
T = T,,(x), U(OT/Oy)~ = 
at y = 4, 
T=T7,, OT/dy = 0°7/dy? = 0 
um and results in 


fourth r ( Ta) 2 ( 6X, ( 2+ BA, ( 


ing the (10) 


where n, = 9/6; A, = Pv,6,/(Gv); and 6, and X, are the dimensional and nondimensional thermal boundary-layer 


thicknesses, respectively. 
Substituting Eqs. (6) and (10) into Eq. (9) and introducing gives the energy equation in the following form: 


(rF)’ = 70Pr(1 + 3d, + 3d,7)/AL(1 + Ad) (11) 
where 
= 1 — (T,/T)) (12) 
— 2772 + 7%) + (180 + 252P — 1087 + — 45PJ? + + 
e fol- (288 — 180J + 36J°)A,7]/(1 + Ad (PJ +A) (13) 
J = = (14) 
The heat-balance equation [Eq. (4)], with the use of Eq. (10), when nondimensionalized and solved for r, be- 
comes 
r = 3A(1 + AD{[1 — + + + (15) 
where 
(6) K = — 1) (Rm/ky) +7) 
ayer a constant. 


The substitution of Eq. (15) into the right-hand side (only) of Eq. (11) results in an equation that may be im- 
mediately integrated once with respect to &. Use of the initial condition £ = 0, \, = 0 gives the integrated expres- 


sion in the following form: 


(7) + CFr = — (P/K)[1 — (16) 
tial where the subscript 0 refers to conditions at £ = 0. Eq. (16) replaces the original heat-balance equation [Eq. 
(15)]. If F=F(é) is known, Eq. (16) is an ordinary first-order linear differential equation whose solution, with the 
initial condition £ = 0, r = 7», is 
ak F dé 
(8) r= + ro’ — (P/K)[1 — (T./T) ae) (17) 
0 
0- From a consideration of Eqs. (8), (13), and (14), it can be shown (assuming J,, finite) that the denominator of 
p- Eq. (17) is equal to zero at = ~. Therefore, in order to enforce the condition that r,, be finite, it is necessary 
he that the numerator of Eq. (17) vanish at = ~. Thus, 
li- = (P/K){1 — — ro’ ce 0 dé (18) 
0 } 
Eq. (18) shows that the boundary conditions for temperature and temperature gradient at the leading edge of 
the plate—namely, ro and 79’—are related to each other mathematically and cannot both be arbitrarily chosen. 


764 JOURNAL OF THE AERONAUTICAL SCIENCES-NOVEMBER, 1952 


From Egg. (8), (11), and (14) at & = ©, one obtains J, = 1—.e., the velocity and thermal boundary-layer thick. 
nesses are equal. The application of L’Hdspital’s rule to Eq. (17) at infinity then gives 


=[1 = 1 (T/T) 


Thus, the reasonable physical result is obtained that the wall temperature approaches the coolant temperature 
at infinity. 5 

The thermal boundary-layer thickness and the temperature distribution along the plate are obtained by the 
numerical solution of two simultaneous first-order nonlinear ordinary differential equations. One of these equa- 
tions is the replaced heat-balance equation [Eq. (16)]. The other is obtained by differentiating the energy equa- 
tion [Eq. (11)]. 


21 + A)(PJ + 


2)’ 
JO 


(19) 


where 


W(A,J) = P(I68S — 27/2 + 7J*) + (180 + 252P — 108J + 21J? — 45PJ? + 12PJ*)\, + 
(288 — 180J + 


O(A,/7) = [336P2J — 108P2J* + 35P2J*] + [336P + (540P + 924P2)J — d40PJ? + (140P — 306P2)J* + 
100P2J*]\, + [(540 + 924P) + (—4382 + 1,512P + 504P?)J + (105 — 1,530P)J? + (400P - 
18SO0P?)J* + 60P?J*)\,? + [(1,512 + 504P) + (—1,224 + 864P)J + (300 — 900P)J? + 240PJ*)\; 
+ [864 — 720J + 
= [168P?J — 81P?J* + 28P?J*] + [336P + (180P + 420P")J — 324P/? + (98P — 216P?)J* + 
+ [(860 + S40P) + (—3824 + 468P + 252P?)J + (84 — 864P)J* + (266P — 135P*) 
+ 48P?J*]\,? + [(936 + 504P) + (—864 + 288P)J + (228 — 540P)J? + 168PJ*]\,3 + 
[576 — 5407 + 
The numerical solution starts at the leading edge; therefore, (,°)o’ must be determined. From Eqs. (7) and 
(19), assuming Jo’ finite, it is found that 
= (35/37) P? J? (20) 
where Jy depends on the boundary conditions 7) and 79’. Reference to Eq. (18) indicates that three different 
combinations of 7) and 79’ are possible: (a) to = 0, to’ # 0; (b) ro ¥ 0, to’ = 0; and (c) to ¥ 0, to’ A 0. How- 
ever, the condition 7) = 0 appears to imply physically that the thermal boundary-layer thickness is zero at the 
leading edge [i.e., (A,)o = O]. 
From a consideration of Eqs. (7), (11), and (16), there results the following equation for Jo in case (a) :* 


Jo®(7Jo® — 27Jo? + 168) — [148/(3P)] = 0 (21) 


SOLUTION OF THE PROBLEM NEGLECTING PLATE THERMAL CONDUCTIVITY 


When plate thermal conductivity is neglected, the original heat-balance equation [Eq. (15)] assumes a simplified 
form, 


The elimination of 7 from Eq. (11) results in an expression that can be immediately integrated with respect to 
£. Use of the initial condition & = 0, A, = 0 gives 


{J?rx7(P(168 — 27)? + 7J*) + (180 + 252P — 108J + 21J? — 45PJ? + 12PJ*)X, + 
(288 — 180J + 3d; + 3A7)(PJ + — TOPE = 0 (23) 


* For cases (b) and (c) one would find Jo3(7Jo? — 27/9? + 168) — (148/P) = 0. 


The solution of the algebraic equation is obtained At the leading edge, the equation defining Jo is ob- 
graphically. A value of X, is specified thereby, from tained from Eqs. (7) and (23) and is 

Eq. (8), giving the corresponding ~ Then a graph of 

Eq. (23) vs. A, for arbitrarily assumed values of \, will Jo?(7Jo® — 27Jo? + 168) — (74/P) = 0 (24) 
indicate the correct A, for the value of A, specified. Re- 

course to Eq. (22) yields the corresponding wall tem- Also, at the leading edge, as (A,)o = 0, Eq. (22) indi- 
perature. cates that 7) = 0, while, from a consideration of Eqs. 
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(7) and (22) and use of the conditions that Jy and Jo’ 
be finite, it is found that to’ = ©. 

This infinite temperature gradient at the leading 
edge of the plate, when axial plate thermal conductivity 
is neglected, is in agreement with results obtained by 
previous investigators—e.g., Grootenhuis and Moore” 
and Yuan.* Grootenhuis and Moore, however, con- 
sidered thermal conductivity across the plate. Their 
results give the axial temperature distribution for vari- 
ous injection ratios but with an infinite gradient at the 


leading edge. Yuan neglected wall thermal conduc- 
tivity entirely and assumed that the wall temperature 
is constant. His results give an ‘‘average’’ wall tem- 
perature for varying lengths of plate. For example, for 
a sweat-cooled plateof length L whose free-stream Reyn- 
olds Number (m/v) = 200,000, Yuan’s method 
states that 7\,/7, will be constant all along the plate 
and will be equal to 0.616 when v,,/u; = 0.001 and P = 
1 (see Fig. 5). 

At infinity, from Eq. (22), r, = 1 — (7./T)) 
while, from Eqs. (8) and (23), J, = 1. These results 
are identical to those obtained for the complete heat- 
balance equation and indicate that results at infinity 
are independent of wall thermal conductivity. 


NUMERICAL SOLUTION OF THE PROBLEM 


When wall thermal conductivity is neglected, the 
simplified heat-balance equation gives to = 0, 7’ = 0. 
Therefore, a solution of the general problem will be ob- 
tained for case (a)—i.e., to = 0, ro’ # O—so that the 
effect of axial thermal conductivity on the wall tem- 
perature distribution may be determined. 


Physical data assumed are P = 1, 7.1, = '5, 
v,/u, = 0.001 and 0.01, r = 1/4, hk = '/41n., m = 100 
ft. per sec., vy = 1/6,000 sq.ft. per sec., a = | 2, RyRy = 
3,360. For these values, A = 0.015752 and 15.752 for 


?»/u, = 0.001 and 0.01, respectively. 

To calculate the temperature distribution along the 
plate, it is necessary that 7’ be uniquely determined 
Eq. (18) defines ro’ in terms of F, but F(£) cannot be 
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obtained until 79’ is uniquely determined. Therefore, 
as Eq. (18) cannot be used to evaluate 70’, arbitrary 
values must be assumed until a physically possible 
axial temperature distribution—1.e., approaching the 
coolant temperature in the downstream direction—is 
obtained. It can be shown from the numerical calcu- 
lations that, at any value of &, 7 is a linear function of 
to. This useful result implies that the numerical solu- 
tion of the differential equations {Eqs. (16) and (19) | 
need be applied to only two arbitrarily chosen values of 
vo’. Solutions for other values of 7)’ are then obtained 
by linear interpolation (see Figs. 2 and 3). 

The numerical solution is begun at the leading edge. 
At this point \, = 0, A, = O and for case (a), 7o = 0, 
while 70’ as determined by the above procedure is 4.437 
for v,/u, = 0.001 and 0.20 for v,/mu, = 0.01; moreover, 
from Eq. (21), Jo = 0.67877 for P = 1. The incre- 
ments of £ used are AE = 0.02 and 0.2 for v,,./u, = 0.001 
and 0.01, respectively. 

Solutions for the simplified heat-balance equation are 
obtained for P = 0.72 and 1 and for v,,/m, = 0.001 and 
0.01. From Eq. (24), Jo = 0.878 and 0.782 for P = 
0.72 and 1, respectively. 


CONCLUSIONS 


The inclusion of the axial wall thermal conductivity 
eliminates the infinite temperature gradient at the 
leading edge of the plate. With wall thermal conduc- 
tivity considered, the reasonable assumption that the 
wall temperature at an infinite distance downstream be 
finite results in a mathematical relationship between 
the boundary conditions at the leading edge of the plate. 
Three possible combinations for the boundary condi- 
tions at the leading edge exist. For the case considered 
in this investigation—i.e., the wall-free-stream temper- 
ature ratio equal to unity and the temperature gradient 
not equal to zero—the boundary-layer thickness ratio is 
less than unity at the leading edge for Prandtl Numbers 
greater than 0.35 and is equal to unity at 2n infinite 
distance downstream; furthermore, at an infinite dis- 
tance downstream the wall temperature approaches the 
coolant temperature. 

When wall thermal conductivity is neglected, the re- 
sulting simplified heat-balance equation gives a wall- 
free-stream temperature ratio of unity and an infinite 
temperature gradient at the leading edge. The bound- 
ary-layer thickness ratio is less than unity at the leading 
edge for Prandtl Numbers greater than '/2. The in- 
clusion or omission of wall thermal conductivity does 
not affect conditions at infinity since, in both analyses 
(i.e., complete or simplified heat-belance equation), 


the boundary-layer thickness ratio is equal to unity 
and the wall temperature is equal to the coolant tem. 
perature at an infinite distance downstream. 


Many of the above results are illustrated in the fig- 
ures. Fig. 4 shows the finite temperature gradient at 
the leading edge when wall thermal conductivity is con. 
sidered. The figure also indicates the greater cooling 
effect of a higher fluid injection ratio; for example, a 
wall-free-stream temperature ratio of 1/2 will be ob- 
tained in approximately one-sixth the distance for the 
injection ratio of 0.01 as compared with 0.001. 


Fig. 5 illustrates the infinite temperature gradient at 
the leading edge when the plate thermal conductivity 
is not considered and readily indicates, for the higher 
fluid injection ratio, the asymptotic tendency of the 
temperature distribution. 


Fig. 6 presents a comparison of the temperature dis- 
tribution resulting from the use of various heat-balance 
equations. For the range of local Reynolds Numbers 
considered and for the physical data assumed for the 
plate, the complete heat-balance equation indicates a 
consistently higher wall temperature than the simplified 
heat-balance equation. 


Fig. 7 presents the distribution of the boundary-layer 
thickness ratio along the plate. It will be noted that 
at the leading edge the boundary-layer thickness ratio 
is larger for the simplified heat-balance equation; also, 
the asymptotic tendency of the boundary-layer thick- 
ness ratio is evident from the results of the simplified 
heat-balance equation for the higher fluid injection 
ratio. It is interesting to observe that the velocity 
boundary-layer thickness is here greater than the 
thermal boundary-layer thickness—contrary to results 
usually obtained.* This is due to the inclusion of the 
second-order heat-balance differential equation [Eq. 
(4)], in conjunction with the variable temperature 
along the plate, into the analysis. 
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The Aerodynamic Forces on Low Aspect 
Ratio Wings Oscillating in an 
Incompressible Flow 


H. R. LAWRENCE* anp E. H. 


Cornell Aeronautical Laboratory, Ine. 


SUMMARY 


A method for calculating the aerodynamic forces acting on a 
harmonically oscillating low aspect ratio wing in incompressible 
flow is presented for the class of plan forms with straight trailing 
edges. The aerodynamic lift and moment resulting from oscilla- 
tory vertical translation and rotation are computed for rectangu- 
lar and triangular wings of aspect ratios ranging from 0 to 4 with 
reduced frequencies from zero to unity. The aerodynamic 
forces arising from the harmonic rotation of 25, 37, and 50 per 
cent chord trailing-edge flaps are computed for rectangular wings 
of aspect ratios 1 and 2 for the zero to unity reduced frequency 
range. The results are shown to be in accord with the Garrick 
theory for plan forms of very low aspect ratio and with the infi- 
nite aspect ratio theory for rectangular wings, In the limiting 
case of zero frequency, the results reduce to the steady-state solu- 
tions of reference 9. A rapid method for solving the integral 
equation is developed for cases where the steady-state solution 
has been expressed in matrix form, 


INTRODUCTION 


_ DEVELOPMENT OF METHODS for calculating the 
aerodynamic forces acting on an oscillating wing of 
finite span in an inviscid incompressible flow has oc- 
cupied the attention of numerous investigators. Theo- 
retical developments and calculations for high aspect 
ratio wings have been presented by Cicala,'! Borbely,’ 
Kiissner,* Jones,‘ Reissner,® and others. An extension 
of the Jones low aspect ratio theory® to nonstationary 
flow was presented recently by Garrick.’ The oscil- 
lating wing of circular plan form was ireated by Schade 
and Krienes* by a natural extension of the known 
steady-state solution. However, the writers have been 
unable to find any reference to the subject of nonsteady 
flow for low aspect ratio wings (say, aspect ratio from 
| to 3) with straight trailing edges. 


The starting point for the present investigation is the 
exact two-variable integral equation for an oscillating 
lifting surface in a linearized incompressible flow. An 
approximate single variable integral equation for low 
aspect ratio wings with straight trailing edges in steady 
flow was obtained from this two-variable integral equa- 
tion in reference 9. Solutions of this single-variable 
integral equation have been shown to agree with wind- 
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tunnel test data for small angles of attack. The two- 
variable integral equation for an oscillating wing will 
be simplified by using the same type of approximation 
as in reference 9 to obtain a single-variable integral 
equation for an oscillating wing. 

As in the stationary case, the result will be applicable 


reference 9. 


only to wings with straight trailing edges. The integral 
equation will be solved by an extension of the method of 


List OF SYMBOLS 


Cartesian coordinates of the point of ob- 


x,y, 8 
servation 

gn Cartesian coordinates of source point 

U free-stream velocity 

w frequency of oscillation 

we/U reduced frequency 

2b wing span 

28(x) local wing span 

2¢ wing chord 

p free-stream density 

A.R, aspect ratio 

o(x, y) perturbation velocity potential on the 
upper surface of the airfoil 

u(x, y) y) 

z = h(x, y) equation of the mean camber surface of 
the airfoil at maximum displacement 

—w(x, y) ikh(x, y) + hr(x, y) = local slope of air- 
foil in stream direction 

dL/dS lift per unit area [see Eq. (3)! 

dL/dx lift per unit chord [see Eq. (13)]} 

f(x) [see Eq. (4)] 

a(x) [see Eq. (11)] 


2p U2m'(x) 
20 U%e'(x) 


gx) 
pl k) 


S*(x, B), R(x, 8), 


rolling moment per unit chord |see Eq. 
(17)] 

rolling moment for extremely low aspect 
ratio wing [see Eq. (16)] 

[see Eq. (20)] 

{see Eq. (21)] 


P(x, k, b) [see Eqs. (23), (24), and (25)] 
G. [see Eqs. (26) and (27)] 
PF [see Eq. (27)]; the Fare the F of refer- 


M(x, k, b), N(R, 5) 


ence 9 
[see Eqs. (44) and (45)] 
matrices defined by Eqs. (32) 


matrices defined by Eq. (34) 


| 
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'b 

n’ dyt 
1 b 

dyt 
1 b 

Le’ dyt 


B(x, k) [see Eq. (47) ] 
Ci, St = cosine and sine integrals (see reference 
14) 


Subscripts x, y, or z denote partial differentiation with respect 
to x, y, or z. 


(1) THe INTEGRAL EQUATION FOR A LIFTING SURFACE 
IN UNSTEADY FLOW 


Determination of the aerodynamic forces acting on a 
wing of finite span oscillating in an incompressible 


t The aerodynamic force coefficients per unit span L;’, 1)’, 
La’, etc., are defined as in A.A.F. T.R. 4798. 


1 dnu(é, 7) ( r ) 


where 


linearized flow has been clearly and rigorously formy. 
lated as the solution of an integral equation by Reiss. 


ner.” A slight modification of the Reissner result wil] 
be used as the basis:for further development. 


Consider a thin wing oscillating harmonically about 
the plane z = 0 in a free stream of velocity L’ directed 
along the positive x axis. Here, x, y, 2 are a left. 
handed set of Cartesian coordinates with the y axis 
spanwise and the z axis normal to the wing, as illus- 
trated in Fig. 1. The 2 coordinate of the mean camber 
surface is given at any instant by the real part of the 
expression h(x, ye. On the plane z = 0+, the verti- 
cal perturbation velocity U’¢,(x, y) may be expressed in 
terms of U/u(x, y), the streamwise perturbation velocity, 


as follows’ 
l (a) 1 1 ’ 
dé dnu(é "(4 r 
OY Jao x—¢ 


where r = Vix — &)* + (y — n)* and the integration is 
taken over the entire plane z = 0. It will be noted 
that Eq. (1) is applicable to the case of a wing with 
swept leading edge. 


y) = 


Boundary conditions on the plane z = O* are given 
in Eqs. (49), (54), and (55) of reference 5 as fol- 
lows: 


o:(x, y) = tkh(x, y) + h,(x, y) on R, 
d(x, vy) = 0 on R, 


where the regions R,, R,, and R, represent, respectively, 
the airfoil, the wake, and the remainder of the plane 
z = 0 (see Fig. |). Upon substituting these boundary 
conditions in Eq. (1) with the following definition for 
w(x, y), there results the equation 


° — 1) 
Oy J; = x— 


—w(x, y) = tkh(x, vy) + h,(x, y) (2.1) 


The lift per unit area is given by the expression® 


dL/dS = 2pU*(¢, + ikd) (3) 


where ¢ and ¢, are evaluated on the z = O* plane. 
ishes at the trailing edge x = 1. 
formulation of the problem. 


The Kutta condition is satisfied if the lift per unit area van- 
Eqs. (2) and (3) and the Kutta condition represent a convenient mathematical 


(II) DEVELOPMENT OF AN APPROXIMATE INTEGRAL EQUATION IN ONE VARIABLE 


A single-variable integral equation will be obtained from the two-variable integral equation {Eq. (2)] by the 
same process as that used in reference 9. For a symmetric spanwise lift distribution, the integral Eq. (2) is multt- 
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plied by a weighting factor V 8°(x) — y? and integrated over the span. It will be observed that no assumption 
has been made regarding the spanwise lift distribution. The integral equation is to be satisfied on the average in 
the spanwise direction, and the weighting factor merely determines the relative importance assigned to various 
spanwise locations. The weighting factor V B? — y? was shown by Flax” to be consistent with the variational 
approach to wing theory and leads to an exact result in the limiting case of low aspect ratio. 
With the definition 
(x) = / dy w(x, y) V B(x) — y® (4) 


B(x) 


the following result is obtained: 


| B(x) * 1 ( 7) 
{(x) = dy V — y’ dé (: + ) 
dy V — y* (: +— ) 


After integrating by parts with respect to y and interchanging orders of integration, the above equation may be 


written in the form 


f(x) = dé dnu(é, n) J. — dnd(1, n) J, (5) 


where /, and J, are defined as follows: 


* B(x) dyy r 
J, = (Y — BAR) — (: — BE) <1 < BCE) (6) 
B(x) dyy r 
(¥ — V B(x) — (: x—€& 
The integral /, may be approximated by the expression® 


The result of using this approximation is an integral equation for low aspect ratio wings in steady flow. This 
theory is known to be in good agreement with experiment for a sufficiently small range of angle of attack. By mak- 
ing use of this same approximation, the integral equation for unsteady flow will reduce to this known steady-state 
solution as the frequency goes to zero. Since the integral defined as J, in Eq. (7) is of the same form as J, defined 
in Eq. (6), it seems reasonable to use the approximation 


x—€& 


Upon substituting the values of J, and J, given by Eqs. (8) and (9) in Eq. (5), the following equation 1s obtained : 


= / dé dnulé, n) [ + 
4 — Big) 
b Ix V (x — &)? + 
die dno(1, n) [ + +1+ | (10) 
4S, =" x—& x—& 


By making use of the definition 


Vix 


+14 


See Str 


B(x) 
g(x) = dyo(x, y) (11) 


Eq. (10) may be simplified to the form 


V(x — + =| 


V (x — + B(x) tk _ ik | 
F(x = a(x — 4 1) 1 l 


(12) 
The aerodynamic forces may be expressed in terms of the function g(x) by integrating Eq. (3) across the span. 
B(x) dL dL 
dy — = = 2pl/*(g'(x) + ikg(x)) (13) 
 @S dx 
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where dL /dx represents the lift per unit chord. The Kutta condition can no longer be satisfied at all points along 
the trailing edge. However, an integrated Kutta condition may still be imposed. The total lift per unit chord 
will be required to vanish at the trailing edge or 


g’(1) + tkg(1) = O (14) 


aspect ratio wings with symmetric lift distribution. 


(IIT) THE SINGLE-VARIABLE EQUATION FOR AN ANTISYMMETRIC Lift DISTRIBUTION 


If the function h(x, y) is antisymmetric in y, the lift per unit area will also be antisymmetric and the function 
g(x) will be identically zero. However, by multiplying Eq. (2) by the weighting factor —(1/2)yV/ 8°(x) — y®, inte. 
grating over the span, and approximating the integral on the right-hand side as in Section (II), there results the 
following equation: 


Vix — £)2 + B%(x) 
e(x) = = m(x) + dém'(&) [ + BX) 


k V(x — +B? 
m(1) [ + (15) 
1 — 


where the functions e(v) and m(x) are defined as fol- The functions e(x) and m(x) may be interpreted 


lows: through the defining Eqs. (16) and (17). Thus, 
1 2pU"e'(x) represents the rolling moment per unit chord 
= dyyux, y) V B(x) — for a wing of extremely low aspect ratio and 2pl/*m'(x) 


ete represents the corresponding quantity for a wing of low 

= ff. dyy(x, y) (lv) aspect ratio. It will be noted that Eq. (15) becomes 
identical with Eq. (12) when f(x) and g(x) are substi- 
The Kutta condition is satisfied approximately as a tuted for e(x) and m(x). Hence, most of the numerical 
weighted average by the expression labor required to solve Eq. (12) need not be repeated 
m'(x) + ikm(x) = 0 (18) in solving Eq. (15). For the case of vanishing fre- 
quency, the integral Eq. (15) reduces to the steady- 

state integral Eq. (47) of reference 9. A comparison 


TABLE I 
of calculated results with test data for the steady state’ 
appears to show satisfactory agreement. 
My be = 
(IV) SOLUTION OF THE INTEGRAL EQUATION 
.216- | + 42861 | -1.888 - 1.5521 | - 
The integral Eq. (12) will be solved by expanding 
= 4235! in a series of V terms, with N arbitrary constants. 
| 1.9681 | + | | A collocation at .V points will be used to determine the 
| 0372 = + -3.750 2.7461 | arbitrary constants. 
50 | 1190+ | - -2.3771 «7908 A series expansion of g(x) which satisfies the Kutta- 
| .530- 3.7024 | | 29.384 - | 5.262 - | 1S 
= — 7.(1)(1 x) 19 
S51 .9094 + 41504 1.560 = 2.2624 550 1,273 g(x u(k)gs( (l + ( 
= 4270 173 | 322- Where g,(v) is a point function defined as the series: 
0125 0513 = 6.2004 2396 49.485 -10.623 = 641321 
4-302 = 6.13 g.(0) = (rf — 0)(Ao + A}) + 
25 2560 = 3.0301 §2 ¢ 11.952 - 1,324 = 3.0551 ° 
0630 = 26 2.627% e547 © 11,5172 (a. A,41) —— (20 
1.00 696 6864 380 + .057 Lhe 0363 ya] r 
The angle was introduced by the substitution x = 
| 
025 0302 Su 6.763 = 5.5084 = 3.816 cos 6. The function p(k) is defined as: 
5 $70 «8874 .07 3.652 = 3.0994 L483 = 1.7261 ‘ 
tk 2k? + ik 
6 Tl u(k) = (21) 


4k? +1 


Eqs. (4), (12), (13), and (14) represent the mathematical formulation of an approximate theory for oscillating low 
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1,00 1.00 e221 - 23954 299 = .2621 - 7594 = .0798- .6951 
2,00 25 e2hh 5.4734 271 3.4554 -641 - 9.0821 | -27.385 - 8.0691 
25 226L = 2.7164 274 = 1.7271 ~10.702 4.5964 6.651 = 
50 279 = 0270 = .8594 -2,480 2.2954 2.0081 
| 
00 125 -0481- 8.4894 +357 = 4.5701 -68,097 -11.0151 | -36.242 -10,5234 
4.00 25 0116 = 2339 - “16.636 = 5.6431 ~8.876 ~ 5.198: 
SO 2221 1.9561 1.1361 -3,877 = 2.8914 «1,995 - 2,559: 
| 4.00 1,00 332 - 2331 - = = 1.4524 = 0256 1,258: 
= 
SPECT PE 
k L, Mz La Mg 
225 +119 = 1,992 922.057 = 1.6964 - 1.8694 
25 0120 ~ =. 9974 3.966 = 
1.00 e117 «250i 316 .2324 2230 
+125 155 3.158 +180 = 1,884 $25,216 = 2,9934 “1.720 = 2.8024 
| .25 2166 = 1.570 .9261 ~5,231 1.536 -3.636 1.107! 
50 17h = 67554 = = 67774 .864 
37 | 1.00 3804 0178 220 = 169 = 
| s 125 6 = 3.3L24 227 = 1.30hi -26. #204322 © 3.6534 
25 0258 = 1.6624 228 .65hi 6.55 2.3071 -2.502 
50 50 2266 68204 0226 = 65204 21.542 1,161 
|* 0276 1654 = 0295 = .585i 
AERODYNAMIC COEFFICIENTS FOR RECTANGULAR WINGS OF 
ASPECT RATIO 2 WITH FLAPS 
| 
ord, Ls M, La 
25 | =.0030L- 4.1714 33.397 1.L074 | -19.489 = 2.3074 
| 
12 | 20286 2.0381] - 1.22hi 8.171 - - 
| 
50 20763 = .97Li} 6154 .L87i -1.20h - .5654 
25 2.00 0122 - 23074 280i = 2264 = .280i 
2% = 1.900 230 2.3364 -39.217 = -18,608 = 3.1861 
37 | 0100 = 2.39hi 227 = 1.1734 =9.573 1.6774 1.609 = 1.7351 
| 
7 | .50 0156 = 22 =2.257 = .9224 “1.101 - 
7 2.00 oo - 2 218 = .29 -. 460 02200 
0125 e170 = 5.3154 293 = 1.717i 242.522 = = 
50 | .25 e2ll - 2.6121 290 - “10.418 2.7391 -3.313 2.27894 
50 | .50 = 1.2021] - #2.410 - 0735 = 1.1364 
50 0329 = .280 - .2208 .7584 .0861- .5664 


TABLE V 
VALUES OF THE FUNCTION R(x, 4 ) 
of / 
jxel 
This function was evaluated only for the values of 4 necessary for the 
rectangular and triangular planforms considered, 
x 3/2 0 -1/2 - 
AB R(x, 4) R(x, 4) R(x) | 
201675 4.00000 
2033L9 253202 
206250 1.997L0 
206698 051977 
212500 1.98966 
213397 47611 
18750 6.02306 
#23325 7.62278 
«25000 7.64326 6.04132 1.95950 
~37500 6.08829 
46651 7.969UL 
250000 8.02718 6.1509 1.64901 -.02718 
«75000 6.31083 
«93302 8.8795 
1.00000 9.01170 6.5032h 1.49676 -1.01170 
1.50000 6.94579 
1.86603 | 11,0652 
2.00000 11.39528 =3.39528 
4.00000 16.7678 930688 4.00000 -1. 30688 


Upon substitution of the expansion (19) for g(x) in 
the integral Eq. (12), the resulting equation can be re- 
duced to the form 


f(x) = S*(x, B) — (1/4)u(k)g.(1)[R(x, 8) + 
P(x, b, (22) 


where the functions S*(x, 8), R(x, 8), and P(x, 6, k) are 
defined as follows: 

1 
B) = 5 + 


/(x — §)? + 


(x — &) 


1 


R(x, B) = 21 +x) + 


\ 
gs (€) 1+ 
1 


/ 
1 V (x — &)? + B(x) 
dé}1+- (24) 
1 


V(x — £)? + 
P(x, b, k) = J - 
1 (x f) 


Eq. (22) will be solved by a collocation method, at V 
points given by x, = cos (nmm/N). 

The right-hand side of Eq. (23) is identical in form 
to the right-hand side of Eq. (43) of reference 9. In 
this reference, it is shown that, after substitution of a 
series of the form of Eq. (20) into Eq. (23), the latter 


. 


can be written as 
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TABLE VI 
VALUES OF THE INTEGRAL M(x,b.k) 
k 
M (x, b,k) - 
M (pq) [1-4 } 
Where p = k(l-x) and q = kb. For convenience the integral is tabulated as a function ss : 
lof the two variables p and q, rather than the three variables, x, b, and k. 42 i. ot 8 1.0 1.2 ah 1.6 1.8 ye 
PART 
+2 | 201405] .01673] .01776] .01829] .01862] .01883} .01899} .01911] .01920] 
3 | 205566] 205246] 206627] 06868] 07033] :07153 :07249 ‘o7315| ‘orm 
sh | .2439 | .3025 | .3285 | .328 .3623 | .3656 .3682 | .3703 | | 22139 | | 22981 | 23259 | | 23590] [3703] | 
| | | | | 1.2 | | 22672 | 13377 | 23854 | | tissu | 25073] “size 
28 «3389 | .4585 | -523u4 | .5637 | | -6439 | -6509 | | 23175 | | | 25228 | 25616 | 25922] 26170] 26375] 
160 | 3637 | -5022 | | .6323 | | 1.6 | 22163 | 23638 | | | | 26650 | [7375] | 
1.2 | .3787 | .5295 | .618h | .6770 +7713 | .7894 | .80L0 | .8160 1.8 | .2362} | .5294 | 6223] .7002 | .7605| 28095] 28500] 8839} “9105 
5436 | .6278 | .7007 6039 | .82L0 .8L03 2.0 2571 | .6858 | .7751 | .8453 .9030 9510 | .9916 1026 
1.6 | «3862 | .5470 | | 7065 S120 | 6327 | | .6632 2.2 | .2728| | .6224 | .7393 | .8367 | .9177 | 9832 | 1.0362 | 1.0847 
8 8 L 6350 | .6972 7966 | .8182] .83L0 | .8470 1246 
1.8 | | | .6350 | .697 2.4 | 22852] | .6564 | .782h | .8901 | .9764 1.048 | 1.2092 | 1.1612 | 1. 29¢7 
2.0 | «3790 | | .6756 -7667 | .7638) «797k | .808L 2:6 | | 5117 | 26617 | 281,7 | .287 | 1.0980 | 1.263L | 1.2195 | 1.2678 
2.2 99 | .5122 .5931 | +7199 | .7330] 7430 | .7509 2:8 | 23005 | 28235 | 26986 | 18363 | | 1.1515 | 122001 | 
| | .5629 | 3.0 | 23037 | 25296 | 27075 | 26478 | [9683 ]1.0668 | | 1.2198 | 1:2805 [173339 
2.6 | | .L675 | .5290 | .5628 +5956 | 25972] .5968 | .595 3.2 | .30h2 | 5307 | .7090 | .8497 | .9706 1.069h | 1.152L | 1.2232 | 1.2842 }1.3368 
2.8 | .3336 | | | .5172 | 5295] .5130 | 5271 | :7038| 8430 | 9625 |1.0600 1.1428 | 1.2113 1.2713 1320 
3.0 | .3207 | | .4573 | 24525 | | 24139 326 | [2985 | | 26929 | | | 121392 | 121863 | 
3.2 | | 13938 | | 24260 3818 | 23622] | 23231 | 4293 
3.b | .2967 | | .3898 | .3839 3151 | .2885| .2622 | .2370 162015 1142511 
3.6 | .2863 | .3512 | .3605 | .3L59 22547 | .2216 | 21892 | 21585 | | 11.1977 
2.8 | .277k | .3339 | .3353 | .3132 2025 | 11635] .1258 | 
4.0 | .2701 | .3198 | .31L6 | .2865 01595 | .1157] .07352 | .03369 
q | 22 | 2b | 26 | 2.8 3.6 | 3.8 22 | 2b | 26 | 28 30] 3.2) 361 3.0] 
| .igh2 | .i9k6 | | .1952 | .1959] | .1962 .2 | -01933] .01938] .01942] .019L9] .02952) .01956 .01958] .o1960 
| 3720 | 23735 | | 23757 3782 | .3788 | .3793 | | 207419] .O7L58] .07L9f .07520] .O75L5] .07567] .0758G .07602} .07618] .07632 
| .5277 | .5307 | .5333 | .5355 | | 25431 | | .1590 | .1603 | .161, | .162h | .1632 | .1639] .2645] .1650] .1656 | .1660 
58 | 26568 | | 26659 | 26696 -6778 | .6800] .6819 | .6837 | .2676 | .270h | £2730] .2751 | .2769 | -2786 | .2800] 22812 | | 
1.0 | .7566 | .7633 | .7692 | .7859 | .7889] .7926 | .79L2 1.0 | | .3984 | | | 2133] :l260] 24283] | 
1.2 | .6259 | .63L3 | .6415 | .8478 | .8660] .869, | .8725 1.2 | .5291 | | .5511 | .5566 | .561L | .5658] 25697] .°730] 25763 
1b | 68648 | | .882h | 28895 .9059 | .9102] .91h0 | .9175 1.4 | .6689 | .6810 | .6917] 27010 | .7091 | £7162] 27227] 27281 | 
1.6 | | | .8929 2 | .9215] .9255 | .9292 1.6 | .8067 | .8230 | .8375| .8501 | .8612 | .8710| 28799] 28877] 28947 | 9012 
1.8 | .8579 | .867¢ | .87L9 | .8817 ~8976 | .9017] .9053 29088 1.8 | .9372 | .9581 | .97 29928 |1,0070 1.0197 | 1.0321 | 1.0413 | 1.0504 | 1.0589 
2.0 8250 -8371 8499 8588 2.0 11.0558 |1.0813 |1.1039 | 1.1237 |1.1L12 1.1568 | 1.1709 | 1.1835 | 1.19L8 1.205 
2.2 | 7571 | .7621 | .766L | .7700 27778 | «7797 | .7829 2.2 |1.1592 [1.1890 |1.215L | 1.2386 ]1.2593 |1.2777 | 1.2943 | 1.3092 | 1.3227 |1.3362 
2.4 | .6809 | .6826 | .6838 | .68L6 6855 | .6856 6855 204 [1.2781 ]1.3080 | 1.3343 2.3578 | 1.3787 | 1.3977 | 2.4147 | 1.4301 
2.6 | .5933 | .5908 | .5883 | .5857 5780 | £5758] | 266 ]1.3101 |1.3468 1.3794 | 1.4082 ]1.43L0 | 1.4781 | 1.4968 | 2.5139 |1.5297 
2.8 | | | | 4605 | 2.8 |1.3550 {1.3939 |1.4285 | 1.4592 ]1.4867 | 1.5237 | 1.5538 |1.5721 | 1.5890 
3.0 | .3887 | .3773 3361 | 3300] .3222 | .3152 3.0 ]1.3791 ]21.4192 ]1.4550 | 1.4668 ]1.5152 | 1.5639 | 1.5848 | 1.6038 }1.6213 
3.2 | | .2671 | .2708 | .2556 -2159 | .20h5] .1937 .1839 3.2 ]1.3832 ]1.4235 ]1.4595 1.b914 ]1.5202 ]1.5158 | 1.5691 | 1.5900 |1.6092 | 1.6268 
3.4 | .2129 | .1903 | .169h | .09AAB] .07082 | .05757 3.4 1.4061 | 1.4746 1.5026 |1.5276 | 1.550k | 1.5708 | 1.5895 |1.6067 
3.6 | .1292 .1018 | .07640] .05256 .00929] =.0270L] -.0L319 f-.05918 326 [1.3373 [1.3750 |1.4086 | 1.4383 1.4885 | 1.5100 | 1.5293 | 1.5469 |1.5630 
3.8 | .05627] 21043 |--.12L9 | -.1437 .1622 3.8 |1.2917 ]1.3268 |1.3578 | 1.3851 ]1.4095 | 1.4507 | 1.4682 |1..987 
|-.1319 |-.1585 |-.1833 |-. =.2276 [1.2346 [1.2662 [1.2940 | 1.3183 |1.3398 |1. 3586 [1.3759] 1.3911 
. 
, ‘ . are tabulated in Tables VI and VII. For @ = 0 (x = 
(x, B) = GAo + G,A, (26) 9 
1), Eq. (20) yields 
where the G are defined as follows: g,(x)|,21 = (Ao + Ai) (30) 


Go(8) 
G,(0) 


In the above-mentioned reference, the F are defined 
and a method for calculating them is also outlined. 

After performing the integration indicated in Eq. 
(24), R(x, 8) becomes 


— 


R(x, = 4+ 2x — V(x — 1)? + + 
V(x + 1)? + + 8 sinh~" 
— 
-] | IS 


The function R(x, 8) has been evaluated for some 
special cases, and the results can be found in Table V. 
The integral (25) is evaluated in Appendix (1); the 
result is 
b, 


P(x, b, k) = be (a (b, k) 
Ci[k(1 — x)] + — Sifk(l — » It) (29) 


where J/(x, 6, k) and N(6, k) are integrals defined in Ap- 
pendix (1). J/(x, b, k) and the function P(x, }, k) 


Substituting Eqs. (26) and (30) into Eq. (22) for 
x, = cos (nx/N), we have 


1 
r=1 
u(k) (Ao + [R(Xn, Bn) + P(%n, 6, (31) 


This set of N simultaneous equations will give Ao, Ai, 

.An-1. With these A,, the series expansion (19) 
satisfies exactly the original integral equation at the .V 
points x,. The G,, are evaluated by the method out- 
lined in reference 9; R(x, 8B) and P(x, b, k) are evalu- 
ated by substituting in Eqs. (28) and (29). 

For (x = 1) the functions R(x, B) and P(x, b, k) be- 
come infinite; however, their sum is finite, and an ex- 
pression for this sum is calculated in Appendix (2). 

Eq. (31) can be written as a matrix equation 


F=(G+HM)A=TA (32) 


where the matrices F, G, H, and A are defined as fol- 


lows: 
F = (4 { f(xo), f(x), (x2), (x3), Fea); (x5) (32-1 ) 


A = {Ao, Ai, Ao, A3, As, As} 


(32-2) 


Gat 


The 


whe 
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AERODYNAMIC FORCES O 


G and H] are square matrices, with 
Gon = F\(6,) \ 99 « 
Hon Thy, = — p(k) { R[x», B(Xn) + b, k) } (32-4) 
Hon Hy, = Hay = Hs, = 0 


where x, = cos 0, and 0, = n/N. 
The braces { } in F and A indicate column matrices. 


(V) SOLUTION OF THE MATRIX EQUATION 


For numerical calculations, a six-point collocation 
method was used, with the trailing edge as one of the 
collocation points. As will be seen in Section (X), the 
results obtained by this method for the case of infinite 
aspect ratio are satisfactory. 

Premultiplying Eq. (32) by 7~', we have 


T'F=A (33) 
and defining new matrices Q and J 


T= (G+ = (1+ = 
(+ = JG"! (34) 


Direct inversion of T (a 6 by 6 matrix with complex 
elements) may be avoided by calculating the inverses 
of Gand J. The elements of the matrix G are real, and 
J may be inverted by the following procedure. 

The matrix J = J + Qis of the form 


J can be partitioned along the dotted lines into four 


submatrices 
Ju 0 
J= 


The inverse of J is 


where 


LOW ASPECT RATIO WINGS 


TABLE VII 
VALUES OF THE INTEGRAL P(x,b,k) 
( J 
yr? : 
K / t J 
Where p+k(f-x)and g*kD . For convenience, the integral is tabulated as 
ia function of the two variables p and q, rather than the-three variables x,b and \ 
R_ E A L fF 

rg <2 1 2 
=.3532 -.6030 =-.7696 -.9698 =1,168) -1,2593 
ol -.1513 -.2830 ~.3859 -05278 - 6893 = 7723 
6 ~.06481 -.16L5 -.2317 23339 - .L6L8 = 51,00 
6 205396 -.1069 -.2282 = 402k 

1.0 203701 -.07L55 -.1078 -.16L1 =. «2506 3119 

1.2 =.0267L =05LL8 -.07926 -.122h - .1933 

1.4 =.01986 -.06038 ~ 209102 = 1524 - 2020 

1.6 =.01570 = 03223 -.04713 -.07376 = 1223 - 1668 

1.8 01229 -.02557 203750 -205901 = 09966 1397 

2.0 ~.01013 ~.020L2 =-030L5 04789 08218 - 1183 

2.2 -.01677 203947 .06855 - .1011 

02 ok 1 2 
21258 +2380 03325 4807 27048 
oh 207968 01555 2238 23381 05263 26999 
6 05708 01) 2k 21640 4132 25710 
6 201,376 208578 21989 23343 4778 

1.0 203 206795 21006 1601 24,069 

1.2 202861, 20552 08213 21316 0232h 23513 

1k «02372 201,566 206826 21103 1979 306L, 

1.6 202025 203827 205752 209359 21706 02696 

1.8 2017L6 003269 01,921 -080L0 o1L85 22391 

2.0 001516 202801 004234 206980 22134 

262 201308 202LL0 203696 206119 


— QO 36 
14+Q] 


The inversion of J has been reduced to a substitution 
into Eq. (36) and to the multiplication of J/x by Ju~', 
or that of a 4 by 2 matrix by a 2 by 2 matrix. The in- 
version of G is equivalent to solving the steady-state 
problem. In most cases, this problem will have been 
solved prior to the consideration of the nonsteady case. 
Its availability will reduce considerably the amount of 
work required to compute the inverse of 7. 

it will be observed that the matrix G depends only 
on the plan form, while the matrix J depends both on 
the plan form and on the frequency of oscillations. 
Hence, for a given wing, it will be necessary to calculate 
G~ only once, while /~' will have to be calculated for 
each frequency. For different camber surfaces G7! 
and J~' do not change; hence, the solutions require 
only the multiplication of 7~' by a new F matrix. 


(VI) CALCULATION OF AIR FORCES 


The superscript 7 refers to translatory oscillations; R, to rotational oscillations; and F, to wings with flaps. 


For translatory oscillations of the wing, F” may be taken as F,,’7 = 


lift and moment coefficients L, and M, are 


Ly = —(2iA.R./k)[(Ao” + Ai™)(1 + tk)(1 — + (tk/2)(Ao” — A2”)] (37) 


= 6°(x,)/b*, and if the A satisfy Eq. (33), the 


My = — ag + 2)(1 + + (; + 


(Ao? — As")(ik — 2) + (Ay? Ast (38) 
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For rotational oscillations about the wing quarter chord, F,” = [8?(x,)/b*]x,. The lift and moment coefficients 
L, and M, are 
L. = i L Ao® Aj®)(1 —2 th Ao® Ms 
Ok -R. | (Ao® + Ai®)(1 + 2u) + — (39) 
Mn — A.R. (Ao® + | (1 + th) ( — Qn) + 2 a ~ + 
ik 
(AF — — 2) + 4 — A;") (40) 
It 
the fol 
being 
produc 
Let a denote the position of the hinge point with re- The 
spect to the mid-chord point. a must either be given | witht 
by a = cos (mr/N) ora = cos + 1)/2N]m. For tion is 
7 other values of a, the lift and moment coefficients may 
| be obtained from the ones corresponding to the above a 
| by int lation. 
y interpolation The 
Case (1) (1)! 
The hinge point is midway between two collocation 
fe) 30 60 90 “ 150 points. F,7" will then have the following values: 
=¢05 
HINGE POINT = B?(x,,)/b? > a) 
D1aGRAM 1. © <a) 
(2) 
In Diagram 1, the hinge point is located at 6 = 75° 
(this corresponds to a 37 per cent flap). The solid 
line indicates the approximate value of the angle of at- 
i acy tack, while the dotted line indicates the correct value. 
Case (II) Th 
& | The hinge point is located at one of the collocation a th 
| points. F,"" will then be 
= B?(x,/b*) > a) 
1 B?(x,) 
30 60 le) 120 150 
x —~-  < 4) TI 
HINGE POINT In Diagram 2, the hinge point is located at 6 = 90°, | from 
DIAGRAM 2. This corresponds to a 50 per cent flap. 
The lift and moment coefficients L, and .V/, can be ob- 
For the case of wings with flaps, an approximate pro- tained by substituting A” for A? in Eqs. (37) and (38). 
cedure will be employed. A clue to this procedure may For rotational oscillations about the hinge point, 
be obtained by consideration of the Prandtl lifting-line | Case (1) 
equation. This equation has been solved approxi- FPF = [B(x,)/b?]x,  (X_ > a) 
mately for cases where the angle of attack distribution F.¥ = 0 (x, <a) 
is discontinuous by choosing the hinge point either at a : qT 
collocation point or midway between two collocation Case (II), Fr"? = [8?(x%n)/b*}xn (Xn > a) 
points.'* The results of such approximate calculations 1 B2(x,) 
are found to be sufficiently accurate for most engineering i? * = - (Xn = a) 
calculations. A similar procedure will be employed a 
The coefficients L, and M, can be calculated by substituting in 
Ll, = - (« + L,— AR. + + — (41) 
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AERODYNAMIC FORCES ON LOW ASPECT RATIO WINGS 


(« + ~ 2u)(1 + ik) + 2ik (j : | + 


1 rR 
5 (Ag™® — A™)(ik — 2) + (A,7* — (42) 


(VII) CoMPARISON OF RESULTS WITH LIMITING CASES 


It can be shown (for method, see reference 9) that the integral equation approaches a correct limiting form in 
the following cases: (1) as k approaches 0, or steady-state case;® (2) as the aspect ratio approaches 0, the limit 
being the extended Jones theory;’ and (3) as the aspect ratio approaches ©, the two-dimensional result" is re- 


produced for rectangular wings. 
The results, obtained by solving the integral equation by the method presented i in the preceding sections, agree 


with the limiting cases, as shown below. This indicates that the proposed method of solution of the integral equa- 
tion is a Satisfactory approximation. 
(VIII) REsuLTS FoR VERY Low Aspect RATIO WINGS 


The Jones theory, extended to nonsteady forces, yields for the lift and moment coefficients: 
(1)* Of rectangular wings: 
L, = A.R.[1 — (i/2k)], M, = (A.R./2)[1 + (é/2k)] 


ALR. 7i A.R. 9i 
L ——{1-—---——)} M, = — 
2 ( k? Ta 2 + 2k? 


(2) Of triangular wings: 
L, = A.R.[(1/3) — (4/2k)], M, = A.R.[(1/3) — (51/12k)] 


A.R. [13 1 Qk ALR. 46:10 91 
M= 
7 2 E 8 ( 


These results would be represented on the plots of L and M by straight lines, tangent at the origin to the results 


of the present theory. 


(IX) RESULTS FOR INFINITE ASPECT RATIO 


In Appendix (3) it is shown that the limit of Eq. (31) when 6 becomes infinite is 


4 N -1 
lim ~ f(x) = Ao +2 D> A,x, — w(k)(Ao + Ai)B,(x,k) (43) 
r=1 


b— @ 
The function B is defined in Appendix (3). The matrix G~ for infinite aspect ratio can be obtained analytically 
from the orthogonal properties of finite trigonometric series and 


mm 4 0 4 0 
2 (/3 -—2 3-2 VD 
= l 0 3 0 -3 3 
12 2 —2 "0 -2 4 —2 
0 4 -3 
| 2 (-Y3-2) 3 -2 173-2 


The results for k = 2 are given in the Table (A) for two-dimensional profiles with 50 per cent chord flap. 


TABLE (A) 
——————Real _Im:z iginary Part — 

Two-dimensional Present Two-dimensional Present 

theory!! theory theory!! theory 

Lh 0.9423 0.946 —0.5129 —0.510 
0.5000 0.500 0 0.002 
La 0.1858 0.191 —0.98405 —0.983 
Ma 0.3750 0.376 —0. 5000 —0.498 
L, 0.45278 0.456 —().41971 —0.414 
M, 0.35610 0.349 —(0.15916 —0.154 
Le —0.03048 --0.040 —0.51788 —0.514 
Mp 0.08903 0.084 —0.35610 —0.355 
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(X) RESULTS OF NUMERICAL CALCULATIONS 


The coefficients L,, M,, L,, and \/, were calculated 
for rectangular and triangular wings for an aspect ratio 
ranging from 0 to 4 and for reduced frequencies from 
zero to unity. The results are tabulated in Tables I 
and II, and they are plotted on Figs. 2, 3, 4, and 5. 

Lines joining points with the same aspect ratio or 
with the same frequency of oscillation were drawn. 

The results of the steady-state theory of reference 9 
(limit of the present theory) were included, and it is 
seen that the aspect ratio lines terminate smoothly at 
the correct values. 

For the infinite aspect ratio case, the coefficients L 
and M were calculated for several frequencies of oscilla- 


1952 


SCIENCES 


NOVEMBER, 


tion. As seen from Table (A), the agreement with the 
corresponding values of reference 11 is satisfactory 
Hence, the results of reference 11 were used in plotting 
the graphs. The constant frequency lines terminate 
at the expected points on the infinite aspect ratio lines. 

In these graphs, the Garrick’ theory gives straight 
lines through the origin. In all cases considered, they 
are tangent at the origin to the curves representing the 
present results. 

The coefficients L,, MW/., Ly, and MW, were calculated 
for rectangular wings of aspect ratio | and 2, with 25, 
37, and 50 per cent chord flaps. The results can be 
found in Tables III and IV. 


(XI) COMPARISON WITH HIGH ASPECT RATIO THEORIES 


The present theory will be compared with Reissner's 
theory,'® since the latter is probably the most accurate 
of the high aspect ratio theories. For simplicity, only 
total lift coefficients will be considered. 

In Fig. 6, the values of kL, predicted by both theories 
for rectangular wings of aspect ratio 3 were plotted. 
It will be observed that whereas the real parts agree, 
the imaginary parts differ by about 15 per cent. In 
Fig. 7, the values of k’?L, for the same configuration 
were plotted; in this case the imaginary parts agree, 
while the real parts differ by 15 per cent. Since the 
Reissner theory is an extension of the Prandtl theory 
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whereas the present theory is an extension of reference 
9 and since these steady-state solutions differ by 15 
per cent, the disagreement noted above is not surpris- 


ing. 
(XII) Limits oF VALIDITY 


The present theory is applicable to swept wings for 
all useful frequencies of oscillation, but it is restricted 
to straight trailing edges and to aspect ratios less than 
4, 

Since the present theory agrees with the theory of 
reference 9 for low k, with the Garrick results for low 
aspect ratio and with the two-dimensional results for 
A.R. = © and for k < 2.5, there is good reason to be- 
lieve that satisfactory values are predicted by it for 
(1) the wing lift and moment coefficients and (2) 
the chordwise lift distribution. 
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Appendix (1) 


EVALUATION OF THE INTEGRAL (25) 


P(x, k, b) = dte~ -1) + V(x + 


(25) 


(x — &) 


After making the following change in variable, § = x + (t/k), P can be written as 


— x) © V2 
P(x, k, 6) = J dte-* (: + ) 
k k(1 — x) t 


Adding and subtracting kb/t, 


— x) 
P(z, &, = —— dte~" (1 + 
k(1 — x) 


(v (b, k) — b, k) + Ci — + 14% — Si ot) 


where N(k, 6) is defined in reference 13 as the integral 


= 1 + 


— x) 
M(x, 6, k) = dte-* (1 + _ vi 
0 


and M(x, k, 6) will be defined as 


45) 
t t ‘ 


Ci and Si are the sine and cosine integrals of reference (14). 


Appendix (2) 
EVALUATION OF R(1, 8) + P(1, k, 


As x approaches 1, R(x,8) approaches the following limit: 


lim R(x, 8) =6—b+ V44 824 dlo 
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ghile P(x, k, 6) approaches 


lim P(x, k, b) = 6 | N(k, b) — log “i - a + =| 


l 
gnce Ci x > — log rae with y = 1.78107 (see Reference 14), 


x20 


lim | Ra, 8) + PU, k, | = 68 — bsinh— + + b log 2ybk + 


Appendix (3) 
Limit oF Eg. (31) When the Aspect Ratio Becomes Infinite 
We have 


4 N -1 
lim ~S,*(x, 8) = B4o+2 ¥ cos BA, 


@ us r=] N 
lim R(x,, 8) = Blog [(x, + 1)/(x, — 1)] 
@ 


lim P(x,, k, b) = Cik(l—x) +72 [5 Sik(l — 


b— ow 


Combining these expressions (8 — 


4 y=! 
lim b) = E +2 A, cos r6 — w(k)(Ao + (46) 


b— « 


where B(x, k) is defined as 


B(x, k) = log | + k(1—x) +2 k — Sik(l (47) 
For x = 1, 
4 N-1 
lim ~—f(1, 6) = A +2 >> A, — wl(k)(Ao + Ai) (10g 2yk +1 | (48) 
b— r= 


In Eqs. (46) and (48), let A, be replaced by A,,’. In the expressions for the L and ./ coefficients, we must now 
substitute A,,’/b; however, for a rectangular wing A.R. = 5, and 6 cancels out. 
For a rectangular wing of infinite aspect ratio, it is only necessary to calculate the A’ and substitute them for the 


Ain the expressions for L and M, taking the aspect ratio as 1. 
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The Study of Compressible Flow Through 
Cascades Using a Linear Pressure-Volume 
Relationship 


David C. Prince, Jr. 
General Electric Company, West Lynn, Mass. 
July 14, 1952 


F sere caer AMOUNT OF EFFORT has been put into the 
study of compressible flows through cascades by approxi- 
mating the relationship between pressure and specific volume by 
a Straight line. In particular, the work of Lin! and of Costello? 
may be mentioned. The relations established by Costello appear 
to allow a relatively simple physical interpretation if a slight revi- 
sion in form is made. 

It is convenient to make use of a hodograph velocity parameter 
w, such that 

1 — 
w = In- — 

1+ V1 — 
where q is the local velocity and a, and M, are the sound velocity 
and Mach number at a reference point in the flow, where the 
straight-line pressure-volume relation is taken tangent to the 
true isentrope. For decelerating cascades, the reference point is 
frequently the undisturbed condition upstream from the cascade. 

Using this definition, a suitably defined velocity potential, such 
that 


(1) 


— = (2) 
and a stream function, 


qaV 1 — M2) = — (3) 
satisfy Cauchy-Riemann equations in the hodograph parameters 
w and (—\A), where ) is the angle made by the flow with a reference 
direction such as the normal to the cascade axis. It follows, of 
course, that (@ + 7y) is an analytic function of (wo — iA). (¢@ + 
iy) and (w — id) can also be expressed, if desired, as analytic func- 
tions of some other parameter £ + in, such as position in the pic- 
ture plane used by Weinig* and others. 

The position variable in the physical plane (x + zy) may be 
obtained from ¢, y¥, w, and \ by the transformation 


1 + idy 


dx + idy = — (dg — iav)| (4) 


which is analogous to Eq. (11.51) of reference 5. 
The flow through a cascade must satisfy continuity and lead 
toclosed airfoils. The continuity condition, 
COS Ay = COS Az (5) 


(where the subscript 1 refers to upstream conditions, subscript 
2 refers to downstream conditions, and p is the density) leads to 
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VARIATION OF HODOGRAPH WITH MACH NUMBER 
APPLICATION OF KARMAN- TSIEN METHOD TO CASCADES 


FIGURE I 


cos A; /cosh w; = cos A2/cosh (6) 


‘which means that the downstream velocity vector cannot be 
chosen independently of the upstream vector. The closure 
condition is to be satisfied by integrating Eq. (4) across the flow 
between a pair of adjacent airfoils along paths upstream and 
downstream from the cascade. The result of such integrations 
will be paths parallel to the cascade axis if the complex velocity 
potential (@ + iy) as a function of the complex hodograph vari- 
able (w — 7d) has a logarithmic singularity at the upstream con- 
dition with an angle 8; representing the ratio of vortex strength 
to source strength, such that 


tan 6B, = —tan A,/tanh a, (7) 


and a corresponding downstream condition with 


tan Bo = —tan d./tanh (8) 


These paths of integration upstream and downstream from the 
cascade will be the same length if the source and sink strengths 
of the logarithmic singularities are adjusted so that 


cos B;/sinh w, = cos 8:/sinh (9) 


Eqs. (6), (7), (8), and (9) provide the basis for comparing cas- 
cade performance at different Mach Numbers. For extremely 
low Mach Numbers, tanh w approaches —1, so that d, approaches 
and approaches In other words, 8; and are the up 
stream and downstream flow angles for an incompressible cascade 
with which the compressible cascade may be compared. — This 
suggests that the proper basis for comparison between cascades 
is the possession of identical, or at least geometrically similar, 
hodographs. Actually, geometric similarity is as much as call 
be specified. Fig. I shows the manner in which geometrically 
similar hodographs vary as the reference Mach Number is varied. 
The functional relationship between complex velocity potential 
and complex hodograph parameter is varied only by translation 
and scale and by a slight rotation (0.3° from M, = 0 to MM = 
1.0). Eqs. (6), (7), (8), and (9) are satisfied for all conditions, 
with 8; and 8; held fixed. These particular hodographs are circles 
with diameter one and one-half times the hodograph distance 
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between upstream and downstream conditions and are symmet- 
rically placed. 

The nature of the relationship between complex hodograph 
parameter and complex velocity potential shows that similar 
yelocity variations will occur at similar points in the complex 
yelocity potential field when the latter has been determined in the 
physical plane As the Mach Number is increased, the upstream 
and downstream flow angles, the stagger, and the stream deflec- 
The camber angle of the airfoils, given by the 


tion all decrease 
The cascade 


diameter of the hodograph circle, also decreases. 
solidity will decrease at the same time, since the potential gradient 
at any point in the physical plane is greater than the stream- 
fynction gradient and since the ratio between the two increases 
with increasing Mach Number. 

The comparison between cascades at different Mach Numbers 
described here is for airfoils with cusp ends. Comparison for air- 
foils with rounded leading edges is possible only if the edge effect 
is neglected, since the relationships among angles really imply 
that the included angle of the leading edge will decrease along 
with the other flow angles in the field. 

The relations obtained here appear to imply those obtained 
by Woolard® using the Prandtl-Glauert approach when the limit- 
ing case of small stream deflections is considered. 
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A New Mechanical Analogy for the Flow of 
Compressible Fluid 


Nobuo Inoue 
Faculty of Engineering, Doshisha University, Kyoto, Japan 


July 25, 1952 


"pea A STATICALLY DETERMINATE STATE OF STRESS in a 
perfectly plastic material. For plane plastic flow parallel 
tothe x, y plane, the state of stress at a point is completely defined 
when the normal and shearing stresses oz, oy, and rz, are known. 
For the determination of these three functions of x and y, we have 
the equations of equilibrium 


and the yield condition 
rT = (2) 
where ¢ = (0; + o2)/2, r = (a; — o2)/2, with the algebraically 


greater and smaller principal stresses o; and o2. In addition, the 
three components of stress can be expressed, with the mean prin- 
cipal stress o and the angle a, which the first principal direction 


makes with the x-axis, as follows: 


= + cos 2a; 
oy = — cos 2am (3) 
Try = sin \ 
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Now we make the state of stress correspond to the flow of non- 
viscous compressible fluid by the following relations: 


d l1— 7’ d da 
qd 2r p T 
(when 7’ < 1), or 
dq 1+ 7’ dp da 
q 2r p T 


(when rt’ > —1), where g, @are the polar coordinates of the fluid 
velocity vector in the hodograph plane, a» is the angle that the 
second principal direction makes with the x-axis, and r’ = 
dr/do. Then it will be seen, after some calculations, that the 
equations of equilibrium |Eq. (1)] become the equations of con- 
tinuity and irrotationality 


(pu) + ) 
u“ = 
Ox oy 


(5) 
v/Ox = Ou/Oy \ 


if the yield condition {Eq. (2)] satisfies the differential equation 


4 (y + 4 
do da 


when the equation of velocity-stress correspondence, Eq. (4), is 
d?r 


4 dr 
"de? da 


when the equation of velocity-stress correspondence, Eq. (4’), is 
used, y being the adiabatic exponent of the fluid. 

Therefore, the two-dimensional steady irrotational flow of the 
fluid can be wholly replaced by the statically determinate two- 
dimensional state of stress in the perfectly plastic material, 
stream lines and equipotential lines of the former corresponding 
to isostatics of the latter and Mach lines to slip lines when the 
fluid flow is supersonic. The shock front in supersonic flow of 
the fluid becomes equivalent to the stress-discontinuity surface 
in the perfectly plastic material, and the jump conditions at the 
shock front (which have been obtained by the conservation of 
mass, Momentum, and energy in passing through it) yield the 
continuity of the normal and shearing stresses acting in the stress- 
discontinuity surface when translated into the terminology of the 
plasticity by Eq. (4) or (4’), the normal stress transmitted across 
a surface element normal to the stress-discontinuity surface being 


d?r 


(6) 


used, or 


dr 


Ytorn +7=0 (6’) 
da 


left discontinuous. 

The yield condition of the perfectly plastic material, the state 
of stress of which yields the flow of the nonviscous compressible 
fluid, can be obtained by solving the differential equation, Eq. (6) 

vor (6’), which becomes, by once integrating, 


— yp) — = constant (7) 


or 


r |1 + = constant (7’) 
respectively, and in Fig. 1 there is shown an example of the solu- 
tions that have been solved so as to coincide on both o- and r- 
axes with the yield condition in ‘‘plane stress’’ according to M. I. 
Huber and H. Hencky’s theory of constant elastic strain energy 
of distortion.! The astonishing resemblance between these two 
yield conditions observed in Fig. 1 and, still further, the more 
precise coincidence that will be seen when the comparison is re- 
stricted locally indicate the possibility of the study of the com- 
pressible-fluid flow by the plastic deformation of the ductile 
metals at normal temperature which are considered to obey the 
above-mentioned theory of strength. Since it is well known that 
the slip lines are sometimes visible on the surface of the material 
subjected to plastic flow without any artifice, in other cases they 
can be brought out by etching; while the isostatics are fre- 
quently observed as lines of fracture on the surface of the mate- 
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Fic. 1. Yield conditions. -~—-—-—-Nonviscous compressible 
fluid [Eq. (7) for ¢ > 0 and Eq. (7’) for « < OO]. --- - Theory 
of constant elastic strain energy of distortion (Ellipse). 


rial failed by cleavage, we can see various types of the fluid flow 
by simply pressing the ductile metals in such a way as we like. 
Although we have confined ourselves so far particularly to the 
plane problem, the analogy can be extended to the axially sym- 
metric flow. However, its detailed discussions will be re- 
ported elsewhere. 

When the perfectly plastic material to be used for the analogy 
does not obey the yield condition satisfying Eq. (7) or (7’), we 
are to deal with a kind of hypothetical-fluid flow, the pressure- 
density relation of which is determined by the yield condition of 
the material considered. Concerning the flow of several kinds of 
hypothetical fluids which have been introduced to represent the 
states of stress of the perfectly plastic materials obeying the 
well-known yield conditions proposed so far, reference should be 
made to my previous papers.” 4 
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Reciprocity Relations in Wing-Body Systems 


W. M. Bleakney and G. K. Morikawa* 
Research and Development Laboratories, Hughes Aircraft 
Company, Culver City, Calif. 


June 20, 1952 


RECIPROCITY RELATIONS (reverse-flow theorems) 
in linearized inviscid flow theory, analogous to the well- 
known Maxwell's reciprocal theorems, have recently been dis- 
cussed by several authors—notably, Ursell and Ward,! Jones,? 
Flax,’ and Heaslet and Spreiter. However, with the exception 
of reference 4, to our knowledge the published discussions have 
been unnecessarily confined to planar systems, possibly because 
of difficulty in prescribing reasonable boundary conditions for 
finite nonplanar systems. We wish to emphasize that, for a large 
class of practical wing-body systems, the body can be represented 
as an infinitely long cylinder (or quasi-cylinder) with the body 
axis nominally in the free-stream direction; in reference 4 the 
body is apparently restricted to zero angle of attack. Since 
sufficient background literature exists—e.g., the above references 
an heuristic presentation is given here. 

Consider an infinitely long cylindrical body with the body 
axis (x-axis) nominally in the free-stream direction and finite- 
area planar or curved surfaces (called wings) arbitrarily dis- 
tributed in the finite space, attached or unattached to the body, 
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but nominally parallel to the body axis. Let ag and ap” de. 
note arbitrary angle of attack distributions of the body axis be 
forward and reverse flows respectively. Then, if 


ap” = ag® = 0 


atx = + o,f the integral reciprocity relation, first demonstrate; 
in reference 1 for planar systems, is valid 


ff [er™ + gr gn] dS = 0 (2) 


where the area of integration S includes both body and Wing sur. 
faces. g and ¢) are the perturbation velocity potentials for 
the forward and reverse flows, respectively, and the subscripts 
x and m denote diffcrentiations in the x-direction and surface nor. 
mal direction n, respectively. Further, conditions (1) can be re. 
placed by 


(| = 


3) 
ap? = = ap? |. = +40 
This is evident since the difference between the boundary valye 
problem given by conditions (3) and that given by conditions (|) 
is the lifting (or antisymmetric) problem where the body is at con. 
stant angle of attack ag“ (constant ag) for the reverse flow) and 
the wings have zero angle of attack. This lifting problem js 
further reducible to two problems: one, in which the body js 
at zero angle of attack having wings with apparent twist (body 
upwash) and, two, in which the body alone is at constant a," 
(constant ag for the reverse flow). The first lifting problem 
satisfies conditions (1) and the second lifting problem satisfies 
Eq. (2) identically. 
We give several interesting consequences of this extension for 
lifting wing-body and wing-body-tail systems, where “‘lift”’ is used 


- to mean cross-wind force in the plane of the angle of attack; 


“sideforce”’ is the cross-wind force at right angles to lift; “com. 
ponent angle of attack” is the angle of rotation of a component, 
such as a section of body with attached wings, from zero inci- 
dence with the free stream; ‘‘system angle of attack” is the cor- 
responding angle for the entire wing-body system; and “‘local 
angle of attack” is the angle of incidence of an element of surface 
area to the free stream. 

(1) For a cylindrical body with wings and with nose and base 
beyond the range of wing interference, the magnitude of the lift 
on a component due to angle of attack of that component is un- 
changed by flow reversal, provided nose and base are inter- 
changed at the same time. 

(2) With the same limitations on wings and body, the side- 
force on a component per unit component angle of attack in any 
plane is equal to the sideforce on the component in reversed flow 
per unit component angle of attack at right angles to that plane. 
It may be shown independently of Eq. (1), from symmetry con- 
siderations alone, that this relationship holds without flow rever- 
sal if the configuration and the component possess a common 
plane of symmetry. : 

(3) The ‘“‘components”’ in (1) and (2) may be the entire system 

(4) Subject to the same limitations on wings and body, the mag- 
nitude of the interference lift from any pair of components is u1- 
changed by flow reversal. 

(5) The wing portion of the lift Lw for a wing-body system 
where the body is at a and wing at zero angle of attack is equal to 
the body portion of the lift Zz" for the same wing-body system 
in reverse flow where the body is at zero angle of attack and the 
wing has the component angle of attack a in the same plane. 

(6) It follows immediately fron (5) that the total lift L©@ for 
the wing-body system where the wing is at a and body at zero 
angle of attack can be expressed as 


L® = Lp® + Ly® = Ly + Ly® (4 
t Actually the conditions can be ag) = 0 at x = —© and ag” = 0 
atx = +o, but we see no advantage gained by these weakened conditions 
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Usually, better approximations are available for the wing portion 
of the lift of wing-body systems. 

(7) The spanwise lift distribution on the wing of a wing-body 
system where both body and wing are at a system angle of attack 
gis equal to the spanwise lift-influence function (cf. reference 5) 
for the same wing-body system in reverse flow where both body 
and wing are at zero angle of attack, excepting an infinitesimal 
chordwise strip on the wing which has a component angle of at- 
tack ain the same plane. If the wing is plane and normal to the 
angle of attack plane, the component angle of attack becomes the 


local angle of attack. 


REFERENCES 

1 Ursell, F., and Ward, G. N., On Some General Theorems in the Linearized 
Theory of Compressible Flow, Quarterly Journal of Mechanics and Applied 
Mathematics, Vol. III, Pt. 3, September, 1950. 

Jones, R. T., The Minimum Drag of Thin Wings in Frictionless Flow, 
Journal of the Aeronautical Sciences, Vol. 18, No. 2, p. 75, February, 1951. 

)Flax, A. H., Reverse-Flow and Variational Theorems for Lifting Surfaces 
in Non-Stationary Compressible Flow, Cornell Aeronautical Laboratory Re- 
port, CAL-42, February, 1952. 

‘Heaslet, M. A., and Spreiter, J. R., Reciprocity Relations in Aerodynam 
ics, N.A.C.A. T.N., No. 2700, May, 1952. 

’ Morikawa, G. K., and Puckett, A. E., Equivalence of the Spanwise Lift 
Distribution to the Lift-Influence Function for Slender Wings and Wing-Bodies, 
Readers’ Forum, Journal of the Aeronautical Sciences, Vol. 18, No. 7, p. 


503, July, 1951. 


On the Kirchhoff Solution for an Oscillating 
Wing in Subsonic Compressible Flow 


John W. Miles 

Associate Professor of Engineering, University of California at 
Los Angeles 

August 1, 1952 


I A RECENT REPORT,! Timman has applied Kirchhoff’s method 
to the boundary value problem for a two-dimensional airfoil 
executing a high-frequency (k > 1) oscillation in a subsonic com- 
pressible flow. It is our opinion that, since he fails to include 
the Kutta condition in his formulation, his end results are in- 
complete. Fortunately, this omission is easily remedied along 
the lines indicated below. 

In the interests of brevity, we shall follow Timman’s notation 
and simply state the boundary-value problem to be considered. 
Let exp (7.\/kx)y be the velocity potential, where 7 is the Mach 
Number and & is a reduced frequency parameter (but not the 
usual U.S. definition). Then y satisfies the wave equation 


+ Vyy + = 0 (1) 
to which a solution is to be found subject to the boundary condi- 
tions 


= w(x), xi <a,y =0* (2) 


y =0, x 


< 
ve + (ik/M)y = 0, x>a,y = 0* (4) 


corresponding to a prescribed downwash (w) at the wing, con- 
tinuity of the potential forward of the wing and continuity of 
the pressure across the wake. The (strong) Kutta condition is 
prescribed by invoking Eq. (4) at x = a (the trailing edge), as 
well as for x > a (the wake). 

The classical diffraction problem differs from that defined by 
Eqs. (1)-(4) only in having the simpler boundary condition y = 0 
forx >a. The Kirchhoff approximation then replaces the bound- 
ary condition y = 0 by yy = Oon y = 0, |x| > a. It appears to 
the writer that this approximation may be expected to have the 
same (frequency ) régime of validity in the airfoil problem only if an 
tigensolution is added to-cancel the singularity exhibited by the 
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pressure at the trailing edge, thereby satisfying (in the weak 
sense) the Kutta condition. 
In order to obtain the required eigensolution, we write 


Ux,y) = vo(x,y) + f “exp [(ik/M) (& x)] dy (5) 


where y7 is Timman’s solution, and x is a solution to Eq. (1) 
satisfying the boundary conditions 


xy = 0, x #a,y=0 (6) 
vz + (tk/M)y = + + x x =a,y =0 
(7) 


The requirement (6) implies that, since Wr satisfies Eq. (2), it 
remains satisfied by Eq. (5), while Eq. (7) is the (weak) Kutta 
condition. 
The required solution for x is evidently a source at x = a and 
is given by 
x = CHo® (x + y?|! 2] (8) 
Substituting Eq. (8) in Eq. (7) and thereby determining the 
constant C, we find that the pressure jump corresponding to x is 


given by 
~ = iV(1 — w(a) exp [ivt + 
p z 


iMk(x — a)| He — x)] (9) 


and just cancels the corresponding term in Timman’s result 
[Eq. (4.1), reference 1]. 

It should be remarked that the end solution for the pressure, 
obtained by adding Eq. (9) to Timman’s result, still does not give 
zero pressure at the trailing edge (strong Kutta condition), as 
would be the case for the exact solution. This is due to the fact 
that the Kirchhoff approximation does not satisfy the boundary 
conditions for |x| > a. Nevertheless, it would appear that the 
suggested modification of Timman’s solution furnishes a logical 
extension of the Kirchhoff approximation to the airfoil problem. 
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Experiments with Large Pitot Tubes in a 
Narrow Supersonic Wake 


N. H. Johannesen and W. A. Mair 
Fluid Motion Laboratory, University of Manchester, England 


July 25, 1952 


INTRODUCTION 


A‘ LOW AIR SPEEDS, the behavior of a pitot tube in a transverse 
total-pressure gradient wasinvestigated by Young and Maas! 
in 1936. They found that the point at which the total pressure 
was effectively measured was displaced toward the region of 
higher velocity. For a number of widely different conditions, 
this displacement was equal to 0.18D, where D is the external 
diameter of the pitot tube. 

In the course of an investigation of two-dimensional supersonic 
wakes, experiments were made to investigate the effect of the finite 
size of the pitot tube. In other investigations, the effect of vary- 
ing the tube diameter has been found to be small, but only rela 
tively small diameters have been considered. The present ex- 
periments were extended to include large tube diameters and 
showed that at supersonic speeds there is hardly any effect of 
the kind found by Young and Maas. 


+ 
— 
(4 
= 0 
100s 


across the wake. The condition of syminetry was taken into ac- 
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TABLE 1 In 
External Diameter, Internal Diameter, ’ and 
D, d, s ben 
Tube In In. D/d Re is mt 
A 0.036 0.026 1.4 \ 
B 0.095 0.068 1.4 
0.375 0.268 1.4 0.68 specs 
D 0.720 0.515 1.4 1 be 
I 0.095 0.015 6.3 ‘ r corre 
F 0.720 0.030 24 H r) in at 
0-67 large 
Al 
EXPERIMENTS 
The experiments were made in a small intermittent supersonic x x from 
tunnel with working section 4 by 5 in., free-stream) Mach Number x the s 
1.96, and Reynolds Number 3.3 X 105 per in. A wedge of 20° x * of tl 
apex angle with °/,-in. chord was set at zero incidence, and its 0.65 a carv 
wake was traversed with pitot and static tubes at a distance 2.5 x x stati 
in. downstream of the trailing edge. The static pressure was oo tion 
sensibly constant across the wake. The diameters of the pitot 0.64 quit 
tubes are given in Table 1. 
In Figs. 1 and 2, the dimensionless ratio ?/// of pitot pressure 
to free-stream stagnation pressure is plotted against distance y 0.63 } \ 
4 


count in drawing the curves. 5 ee \ 
The curves for Tubes B and E are not shown in the figures, but 062 -- C-#- : 
they both coincided with the curve for Tube A. 5 x 
The slight differences in the levels of the curves outside the 
wakes are of the same order of magnitude as the day-to-day vari- 0-61 i Cor 
ations of Mach Number in the tunnel (about '/; per cent). The -0.18 -010 -0.05 oO 0.05 O10 Os F 
pressures near the edge of the wake were rather unsteady, and y inches 
average readings were used. It will be seen that the curves have, Fic. 1. FB 
maximums just outside the wakes. This is because the stream Prof 
lines in these regions have passed through the fan-shaped foot of July 
the shock wave behind the trailing edge, where the compression 
is more efficient than in the single shock wave further away from 
the wake. [' 
D 
ISCUSSION OF RESULTS 
Fig. 1 shows that, for a constant value of the ratio D/d (= 1.4), the 
variations of external diameter up to about 0.4 in. have a negligible 0.69 > es tnt 
effect on the apparent width of the wake. However, when D is PP nol ons 
about 0.7 in., the wake profile is considerably distorted. These pooh / A tial 
tube diameters should be compared with the actual width of the 0.68 oe r ; criti 
wake (about 0.1 in.). The difference in the peak values of P/H \ ps 
at the center of the wake for Tubes A and C is perhaps not sig- ' 
nificant and may be caused by a small change in free-stream Mach 0.67 H Thi: 
Number. well 
Fig. 2 shows a comparison between Tubes A and F, having ap- pap 
proximately the same internal diameter but with external diam- 066 not 
eters differing by a factor of 20. The curves show that the effect In t 
of variation of D, with constant d, is small indeed. (In this case, ; whi 
also, the difference in the peak values at the center of the wake - 0465 plat 
may not be significant. ) bein 
Comparison of Tubes D (Fig. 1) and F (Fig. 2) shows that for = 
large values of D the internal diameter is of major importance. 0.64 - 
As already mentioned, the curves for Tubes B and E did not ala 
show any significant difference; thus, the effect of varying d is A 
negligible for small values of D. 063 ; Ten 
The results may be summarized as follows. For values of D é 
up to about 0.4 in. (i.e., about four times the wake width), the & «safe 5 a whe 
effects of variation of D and d are hardly detectable. With a 062 F- (1’) 
tube having an external diameter as large as 0.7 in., the wake ? === F 
profile is considerably distorted if the internal hole is correspond- 
ingly large (D/d = 1.4), but with a small internal hole even this 0.61 L 
large tube gives almost the same wake profile as a much smaller -0.15 -O0.10 —-0.05 ° 0.05 O10 0.15 it st 
tube. It was not possible to investigate much larger values of D y inches deci 


without choking the tunnel Fic. 2. 


Bae 
j 
Fr 
poise 
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In these experiments no displacement effect as found by Young 
and Maas could be detected. Hence, any effect of this kind must 
beno larger than the scatter of the experimental points and, thus, 
j;much less than the displacement at low speeds. If a tube as 
large as Tube C were used to explore a wake of this width at low 
speeds, the Young and Maas correction to the wake width would 
be of the same order as the wake width itself. (However, the 
correction given by Young ind Maas would require modification 
in an extreme case such as this, where the tube diameter is much 
larger than the total wake width. ) 

Although the explanation of these results is not understood, a 
simple consideration of the effects of the normal shock wave in 
front of the pitot tube shows that the problem is quite different 
from the corresponding one at low speeds. Immediately behind 
the shock wave there must be a peak of high velocity at the center 
of the wake. At the same time there must be considerable 
curvature of the stream lines associated with the variations of 
static pressure behind the shock wave. Thus, the initial condi- 
tions for the (subsonic) flow into the mouth of the pitot tube are 
quite different when the free stream is supersonic. 
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Computing Critical Loads by Splitting the 
Rigidity of Composite Systems 


F. Buckens 
Professor of Applied Mechanics, University of Louvain, Belgium 


July 18, 1952 


[' SEEMS INTERESTING, when noting Bijlaard’s papers on the 
stability of sandwich plates,' to point out that the splitting 
considered there occurs, in fact, in two different ways. The 
first stage is really a splitting of rigidities; in partial system (1), 
the core retains its rigidity, as well as the external plates, except 
for its flexural rigidity, while the reverse is assumed for the partial 
system (2). If P; and P, are the first buckling loads of these par- 
tial systems, their sum provides an inferior limit for the first 
critical load of the real system, 


Per = Pi + Pr (1) 


This is a particular case of Temple’s generalization of South- 
well’s well-known method.? The second stage of splitting in the 
papers of reference 1 consists of dividing the ‘‘elasticity’” (and 
not the rigidity) of the partial system (1) into two subsystems. 
In the first (1’), only shear elasticity of the core is considered, 
while in the second (1”), longitudinal elasticity of the external 
plates alone is taken into account (flexural rigidity of the core 
being neglected). Still, without rigorous proof, an analogous 
procedure had been used by Bijlaard in a previous paper® to re- 
establish Engesser’s approximate formula for the critical load of 
alatticed structure. 

As can be found along the lines of a general theorem giver: by 
Temple,? one gets the following inequality: 

When P,;’ and P," are the critical loads in the partial subsystems 
(1’)and (1"). 

For example, if the differential equation takes the form 


(d2y/dx?) + |Py/a(x)] = 0 


it suffices to admit that the inverse of the rigidity a(x) can be 


decomposed : 


1/a(x) = [1/a’'(x)] + [1/a"(x)] 
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In this manner, combining Eqs. (1) and (2), one obtains 
Foe + + (3) 


By virtue of Rayleigh’s theorem, inferior limits (1) and (2) usually 
provide good approximations at the same time. 

Besides including Dunkerley’s formula for the decomposition of 
the mass in problems of natural frequencies or whirling speeds, 
Temple’s method applies tu systems whose “‘internal’’ elasticity 
can be subdivided into several partial systems and yields 


(4) 


where \, and A, represent, respectively, the fundamental eigen- 
values of the real system and of the kth partial system; for in- 
stance, \ is the square of the natural frequency (w?) of the critical 
load (P). 

Eg. (4) had been independently established by the author, with 
some more generality, by a decomposition of “influence coeffi- 
cients.”"* By this method, “external elasticity’’—i.e., elasticity 
of the boundaries—can also participate in the decomposition in 
order to form partial systems that are easier tosolve. The papers 
in reference 4 give several examples for vibration and stability 
problems; more generally, when the eigenvalue problem can be 
put into the form of an integral equation, 


wx) = (x, dé (5) 
where A(x,) is a symmetric kernel that characterizes the given 
system, if pattial systems with kernels A;(x,¢) can be derived such 


that 


= K (xe) (6) 

then the fundamental eigenvalue \; of the real system can be 
assigned an inferior limit of the form of Eq. (4). 
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On Heat Transfer, Recovery Factors, and Spin 
for Laminar Flows* 


Arthur N. Tifford and Sheng To Chu 
The Ohio State University, Columbus 
July 25, 1952 


nN Marcu, 1950, AN O.A.R.t-SPONSORED RESEARCH PROGRAM 
I on the effect of surface pressure and temperature gradients on 
the characteristics of the laminar boundary layer was begun at 
The Ohio State University. The initial effort was directed at 
completing the computational coverage of the analytic results for 
wedges. The Computation Section of the O.A.R. at Wright-Pat- 
terson Air Force Base did the computing on their Reeves analog 

computers. 
About 25 per cent of the computations of references 1 and 2 
The two sets of 


deal with the same values of the parameters. 
Thus, the two 


function values obtained agree within | per cent. 
systematic computation programs conveniently supplement each 


other. 


* Work supported by The Flight Research Laboratory, U.S AF 
t Office of Air Research. U.S A F 
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TABLE 1* 
TEMPERATURE RECOVERY ON INSULATED SURFACES 


o = 0.5 0.7 1.0 2.0 
B (Ta — 
—0.1988 0.706 0.844 1.004 1.428 
—0.14 0.714 0.850 1.016 1.443 
0 0.718 0.846 1.018 1.442 
0.5 | 0.712 0.846 1.012 1.418 
1.0 0.710 0.844 1.008 1.424 
2.0 | 0.708 0.840 1.006 1.424 
a'/2 = 0.707 0.837 1.000 1.414 
*Tq = temperature acquired by insulated surface 
7: = local free-stream temperature 
ui = local free-stream velocity 
Cp = specific heat at constant pressure 
g, J = dimensional constants 
= Prandt! Number 
8 = pressure gradient parameter 


We were interested to note in the report of reference 3 that the 
calculated recovery factor for a separating laminar boundary layer 


is significantly lower than for other laminar boundary layers. 
In Table 1, we have reproduced our earlier analog-computed re- 
sults. On the basis of these, we had concluded that Pohlhau- 
sen’s square-root rule is just as valid for gas flows involving local 
surface pressure gradients as for the classical flow along a flat plate 
aligned with the gas stream. No exception was made of the sep- 
arating boundary layer. 

In attempting to resolve this conflict, we reviewed the numerical 
integration method used in reference 3. It is simply a direct ap- 
plication of the automatic integration method of reference 4, 
supplemented by the procedure, described in these pages in refer- 
ence 5, for solving two-point boundary value, linear, total dif- 
ferential equations without recourse to trial and error. The in- 
tegration procedure of reference 4, however, does not allow an a 
priori specification of the boundary value at infinity of the sup- 
plementary function, g. A study of the properties of the primary 
sfunction, 6, near the outer edge of the boundary layer is neces- 
Uary to determine just what that boundary value properly is. 

nly after it has been shown that @ approaches zero faster than 
exp(—o fof dn), as n approaches infinity, is it certain that g ap- 
proaches zero there. 

In studying this question, we have related the g function asymp- 
totically to the Whittaker function.6 The boundary value, zero, 
has been confirmed. We conclude that nothing is wrong with 
the integration method used in reference 3. The possibility al- 
ways exists, however, that the coding of the problem for the ma- 
chine—or the machine itself—is in error.* 

As a by-product of the above study, we noticed that a differ- 
encing method of integration with a precision comparable to that 
of the method of reference 4, and without the disadvantage at 
infinity noted above, is readily obtainable. Thus, by adding and 
subtracting the two Taylor series, 


On = On-1 + + (€?/2) + (3/3!) + 
+ (€?/2) — (€7/3!) + 
(e4/4!) 


On-2 = On-1 — €On 


we obtain 


= (1/€®) — 20n—1 + On-2 — — .. .] 
= (1/2e) [On — On-2 — — ...] 


respectively, where e is the size of the step in the integration. 


* In order to throw more light on the subject, we have performed some 
numerical integrations with a precision of about '/2 per cent and obtained a 
recovery factor of 0.836 for the separating boundary layer, Thus, we check 


the analog computation of the O.A.R. within 1 per cent. 
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Substituting these equivalents for @,—1” and @,-1’ into the Original 
second-order differential equation, we obtain a difference equa- 
tion for @, in which the dominant extrapolation error is propor. 
tional toe‘. The application of this method to the computations 
of reference | involves no additional complication and reduces 
the dominant error from to 
The improvement is wholly due to the more correct difference rep- 
resentation of the first derivative, 0,—1’. 

We should like to report here, too, one of our earlier analytic 
findings: As an interesting application of the order-of-magnj. 
tude argument, note that, for large values of y, the second term 
in the equation 


(1/0)0” + fo’ — (2 — = 0 


may be neglected as compared with the third. It is now evident 
that the effect of large values of 7 on the thermal boundary-layer 
function is the same as that of large Prandtl Numbers. In other 
words, the thermal boundary layer is thin as compared with the 
velocity boundary layer. Accordingly, the velocity profile within 
the former may be taken as linear, in general, and as parabolic jp 
the case of a separating boundary layer. 

The analytic solutions for @ are easily obtained under these cir- 
cumstances. They involve modified Bessel functions of the 
second kind (see reference 7). The solid curve of Fig. 1 repre- 
sents the universal solution for the nonseparating boundary layer 
for large positive values of y. Superimposed are data obtained 
from the calculations of Chapman and Rubesin® for the case of the 
flat plate. It is seen that, for y values near 10 (and larger), the 
Bessel function representation of the temperature profile is excel- 
lent. 

Similar analyses for large negative values of y also lead to Bes- 
sel function solutions. However, in these cases @ is found to oscil- 
late between positive and negative values as one moves away 
from the surface. Furthermore, there is no well-defined thermal 
boundary layer. It has been concluded that the solutions for 
large negative values of y are not physically realizable. It ap- 
pears in fact, in contradiction to the assertions of reference 1, 
that, when the value of y (2 — 8) is assumed more negative 
than —1, the solutions are physically unreal. In such cases the 
analysis would require an infinite heat source at the origin. 

Another possible application of some of the computational re- 
sults is the determination, by means of them, of the frictional re- 
sistance to spin of axisymmetric bodies. This extension of the 
analysis of reference 9 will be elaborated upon in a forthcoming 
U.S. Air Force Technical Report. 


Fig 1. Universal Besse! Funchon Representalion of Thermal 
Boundory Layer Characterislics at Large Values of "¥" 
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Remarks on ‘‘On the Flight Dynamics of 
Slender Special-Purpose Aircraft’’ 


Niedenfuhr 
Vibration and Flutter Unit, Aerophysics Laboratory, North 
American Aviation, Inc., Downey, Calif. 


June 16, 1952 


ie AUTHORS OF REFERENCE | consider a nonrolling aircraft 
and remark in a footnote: ‘It is not an entirely trivial mat- 
ter to show the nonexistence of the roll equation because Eq. (4) 
shows that a roll acceleration does exist even when p(t) = 0, 
and there are undoubtedly nonvanishing forces due to many of 
the variables that contribute to the roll torque.” 
Eq. (4) of the paper, when considered in detail, is seen to be in 
error, however. Repeating the essentials of the derivation, one 
has: . 
The velocity of a point ?, located at 7 is 

= (dr/dt) + a7 (a) 


(not? = so that 


v = (2+ + 
xr +6 » (b) 
—xq +r 


7] 


The angular momentum \, is given by 
h =7dmi (c) 
so that neglecting products of vibrational displacements, one has 


+ (—xqg + 7) — (s + 7) (xr + dm 
= + 7) + — (y + — x(—xg + 7)} dm} (d) 
= {x(xr + 6) — (y +o) [(s + — (y + dm 


Then, 


dL = (ddz/dt) + — (e) 
dL, = | —xyq + yr — — — 2xf — 26 — — 
+ gx(xr + 6) — (o + y) (2 + + (9 + 

(s+ r)%q + (o + y)(s + + rx(—xqg +7)} dm (f) 
dL. = | —xq(y +o) + 7) + — + 
{y+ o)(s + (r? — gq?) + rgl(y + — (2 + r)?|} dm (g) 


L = fia (h) 
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The elements of dZ given by Eq. (g) are seen to be the same as 
obtained by the authors, except for the presence of y7 and sé and 
the absence of gxo and rx#. Integration over the length of the 
airplane of yr and c@ is zero by the authors’ Eq. (2), and gxo@ and 
rxr are cancelled out completely in Eq. (f) above. For an air- 
craft that is truly slender and symmetrical, C = B, so that Eq. 
(4) should read 

a, = 0 (i) 

The question of whether or not to include the roll degree of 
Roll torques do exist because 
If the rigid 


freedom remains to be examined. 
of the other degrees of freedom, as the authors assert. 
body coordinates, p, g, r, “, v, and w are restricted to be quanti- 
ties of perturbation order (whose squares and products can be 
neglected) as is done in the conventional stability analysis, Eq. 
(4) can be written 

ag (ji) 
and the gyroscopic couplings and terms involving Mathieu func- 
tions referred to in the first part of the article will all be seen to 
be of higher order. To the first order in small quantities, then, 
the equations of motion for all degrees of freedom (roll and twist 
included) of a flexible aircraft will be a set of linear differential 
equations with constant coefficients. 


REFERENCE 
! Rosenberg, R. M., and Stoner, George, On Flight Dynamics of Slender Spec- 
ial-Pur pose Aircraft, Journal of the Aeronautical Sciences, Vol. 19, No. 1, p 


24, January, 1952. 


Authors’ Reply 
R. M. Rosenberg and George Stoner 


We wish to thank Mr. Niedenfuhr for his illuminating re- 
marks. Insofar as they pertain to Eq. (4) and the footnote of 
our paper they are quite correct. In particular, we agree that 
the acceleration term in the roll equation vanishes identically for 
the case considered by us. However (and Mr. Niedenfuhr 
agrees), roll torques due to degrees of freedom other than roll re- 
main. In view of this, we should still have added a footnote 
(reflecting essentially the sense of the one published) in order to 
warn against ignoring an equation that appears to have nonzero 
members. Yet, the actual demonstration of deleting the roll 
equation from the system is so simple as to deserve no more at- 
tention than a remark in a footnote. Fortunately, Mr. Nieden- 
fuhr’s correct observation does not invalidate the results of our 
paper. 

While agreement exists between the discusser and us regarding 
the footnote in question, we believe that a further point raised by 
Mr. Niedenfuhr merits more detailed examination. The terms 
involving ‘“‘gyroscopic couplings and . Mathieu functions’”’ 
(the latter should read ‘‘periodic coefficients”) are said to be of 
higher order. As a matter of fact, the indicated terms are of 
higher order on/y when rolling displacements of perturbation 
magnitude are considered. If the investigation is concerned with 
perturbations in roll superimposed on steady roll—i.e., if we are 
searching for solutions in the neighborhood of steady-state solu- 
tions which include roll—then first-order terms with periodic 
coefficients do arise (on the simplified premise that the plane of 
vibrations remains fixed in space while the missile rolls). How- 
ever, let us assume with Mr. Niedenfuhr that the terms with 
periodic coefficients are of higher order. We believe it to be 
dangerous advice to conclude that these terms may be ignored in 
a conventional stability analysis. 

Consider, for instance, the equation 

+ a’x = x? sin Dt (1) 


It will be noticed that the equation contains a 


(a, b=constants). 
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second-order term with a periodic coefficient of frequency b. 
If we substitute 


x =x + << |xo| 
for all ¢, where xo is some value of x and ¢€ is a time-dependent 


variable in a small neighborhood of xo, then the equation becomes 


where higher order terms in £ have been neglected. 
Now, the related equation 
+ = xo? sin bt — a2xo (3) 
has the steady-state solution 
= [xo?/(a? — b?)] sin bt — xo (4) 


Since Eq. (4) is the solution to Eq. (3), it cannot satisfy Eq. (2). 
Let us search for solutions to Eq. (2) of the form 


= E(t) + n(t) (5) 


Substitution of Eq. (5) in Eq. (2) will show that n(¢) is the solu- 
tion of the Mathieu equation 


7 + (a? — 2xo sin bt)n = [x03/(a? — b?)] X 
(1 — cos 2ht) — 2x? sin bt (6) 


Since the solution of Eq. (1) is now 
x(t) = xo + E(t) + n(t) (7) 


it follows that x(t) is stable or unstable with (¢). This observa- 
tion does not include the possibility that x(t) is also unstable with 
§(t) for the possible case (by no means excluded) that a = 6. 
Had we neglected the higher order term in Eq. (1), we should 


have found 
x = Ai sin (at + ¢) (8) 


A, ¢ = constants, which is never unstable. Therefore, the de- 
letion of higher order terms with periodic coefficients is not ad- 
visable in conventional stability analyses. 

Lest it be thought that this conclusion depends on the implied 
condition x) # 0, we indicate a similar result that does not de- 
pend on x. Let it be supposed that the solution to Eq. (1) lies in 
the neighborhood of that to the mutilated equation which results 
when the second-order term of Eq. (1) is discarded. Then, the 
solution to the actual equation may be written as 


x = A sin at + ¢(t) (9) 


and powers of ¢ higher than the first may be neglected. ( Diffi- 
culties arising at or near a zero of A sin at are easily overcome by a 
shift of the ¢-axis.) It can be shown without difficulty that, 
now, ¢(¢t) must satisfy the Mathieu equation 


+ {a* — 2A [cos (a — b)t — cos (a + 5) = 
(1/2)A?{sin bt — (1/2) [sin (2a + b)t — sin (Qa — b)t}} 
(10) 


Again, the solution of Eq. (1) is seen to be stable or unstable 
with that of a Mathieu equation. 

Returning to the case with the solution as given in Eq. (7), 
another significant difference emerges between Eqs. (7) and (8) 
if our concern is the general behavior of the solution in the neigh- 
borhood of x9 rather than merely its stability. The former re- 
veals the presence in the solution of a periodic term of frequency 
b having the large amplitude x»9?/(a? — b?), while the latter, ob- 
tained by a “linearization” that simply ignores higher order 
terms, is of frequency a. 

These considerations throw serious doubt on the view that the 
solution of a set of linear differential equations with constant 
coefficients reflects the behavior and stability of a flexible air- 
craft that is permitted to roll. These same thoughts have led us 
to consider a system lacking periodic coefficients altogether rather 
than one where their presence must be ignored if the equations 
are to remain tractable. 
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Reply to Fettis’ Comments on ‘‘Aerodynami 
Coefficients of an Oscillating Airfoil in Two. 
Dimensional Subsonic Flow”’ 


R. Timman, A. |. van de Vooren, and J. H. Greidanus 
National Aeronautical Research Institute, Amsterdam, Holland 
July 22, 1952 


numerical data published in the subject paper.?  Fettis re. 
marks that it is difficult to believe that our coefficients are more 
accurate than those of Dietze for the following reasons: 


ig THE May Issue of the JOURNAL, Fettis comments! on the 


(1) Dietze’s results have been confirmed independently by 
Fettis’ method. 


(2) The quantity [LzJe = ka’/w* should, according to Fettis, 
approach the value 1 for w ~ ©, while our coefficient has () as 
limiting value. 


(3) According to Fettis, the relation 
ka = tw(mq + Rp) 
should be satisfied, and this is not the case for our coefficients. 


With regard to the first remark, it may be emphasized that both 
Dietze’s and Fettis’ results are obtained by an approximate solu- 
tion of Possio’s integral equation, while our coefficients haye 
been computed by aid of an exact analytical theory. As neither 
Dietze nor Fettis gives exact estimates of the errors in his results, 
the accuracy of his coefficients remains doubtful. In fact, it is 
known that Dietze’s iteration procedure does not converge if 
w exceeds 0.7 (and JJ = 0.7), while the approximation of the 
kernel used in the evaluation of certain integrals may very well 
generate convergence to an erroneous value if w is smaller than 
The different approximation of the kernel by Fettis may 
have a similar effect. 

Concerning the second remark, it must be noted that the 
asymptotic value |Lale = 1 ifw = © (corresponding to the vir- 
tual mass of the airfoil) only holds in incompressible flow. In 
compressible flow, however, an oscillating airfoil behaves like a 
moving acoustic radiator, and a mathematical approximation 
method,’ valid for high frequencies and based on common assump- 
tions used in the diffraction theory of acoustic waves, confirms 
the fact that [Lale ~ Oifw —~ ©. That the concept of virtual 
mass is essentially modified in compressible flow has also been 
shown by Miles! and Longhorn.’ Curves of against 
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for different values of the Mach Number (Fig. 1) clearly show 
that [Lalrk > 0 if w > ~ and M # 0. The approach of the 
curves toward the origin is oscillatory, which is in accordance 
yith reference 3. The curve for M = 0.35 closely follows the 
curve for WV = 0 if the latter is supplemented by the [Ly |p-axis 
between 1 and 0. For the larger values of M/, the deviations 
fom the incompressible case increase, and it is obvious that Fet- 
tis’ extrapolation of Dietze’s values is incorrect. 

Thirdly, with regard to the identity kg = t#(mq + kp), our 
exact analytical expressions for the coefficients show that this rela- 
tion holds only in incompressible flow. As Fettis’ derivation of 
this “identity” for compressible flow is not yet available to us, 
no further comment can be given. 

In conclusion, we cannot admit that any coincidence of errors 
in our calculations has been brought to light. We maintain our 


daim that the accuracy of our results amounts to a few units in 


Gas-Turbine Power Plant for Aircraft 


Ali Bulent Cambel 

Assistant Professor of Mechanical Engineering, State University of 
lowa, lowa City 

July 24, 1952 


ECENT IMPROVEMENTS IN PROPELLER EFFICIENCIES once 
R again bring to the foreground propeller-driven aircraft. 
At least theoretically, a power plant suggests itself which might 
be used to some advantage in such aircraft. 

If one analyzes the performances of the various gas-turbine 
cycles as summarily shown in Figs. la and 1b, one notices that, 
for the same temperature limits, a regenerative cycle in which 
both the compression and the expansion take place isothermally 
has the highest thermal efficiency and that the efficiency of this 
plant is not affected by the pressure ratio. Further, one notices 
from Fig. lb that the highest specific output parameter is achieved 
ina plant having isothermal compression and expansion. The 
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the last digit supplied, a claim based on analytical checks for the 
limits of all errors committed. 
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analysis of this cycle indicates that, in the absence of losses, the 
need for the combustion chamber is obviated. Thus, the only 
external heat that must be supplied to the plant is that neces- 
sary to maintain isothermal expansion in the turbine. 

Although the advantages of isothermal compression and ex- 
pansion are evident, in practicality their achievement is diffi- 
cult. One method that suggests itself is to bleed gas from the 
compressor and, after raising its temperature, introduce it into 
the turbine. Of course, the air entering the turbine will have to 
be carefully apportioned, but this should not pose undue diffi- 
culty. The bleeder gas may be heated suitably with an atomic 
reactor, which may also be used to preheat the throttle gases 
entering the turbine. A schematic layout of the power plant is 
shown in Fig. 2, in which the isothermal compression is accom- 
plished by multistaging and intercooling. Because there is no 
need for combustion in the cycle, the use of air as the working 
medium is not necessary. A light gas, highly desirable because it 
would reduce the size of the plant, may therefore be used. In 
passing, it might be mentioned that, everything else remaining 
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Efficiency curves of certain gas-turbine cycles. 


Fic. la (left). 
Combustor temperature, 1,500°F. 


turbine cycles. 


compression and isothermal expansion with regeneration ; 


Fic. lb (right). 
No loss in available energy is considered. 
expansion; (3) isothermal compression and isothermal expansion; (5) adiabatic compression and isothermal expansion; (7) iso 
thermal compression and adiabatic expansion; (2) adiabatic compression and adiabatic expansion with regeneration; (4) isothermal 
(6) adiabatic compression and isothermal expansion with regeneration ; 
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Specific output parameter curves of certain gas 
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(8) isothermal compression and adiabatic expansion with regeneration. 
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Fic. 2. Arrangement of reactor gas-turbine power plant. 
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the same, the characteristics of the gas do not affect the efficiency 
of the plant. Helium suggests itself as a reasonable Possibility be. 
cause it may be obtained conveniently and is safe. Besides being 
advantageous as far as the size of the plant is concerned, heligg 
is also desirable because of its high thermal conductivity as com. 
pared to air. No complexity is expected with the design of the 
heat exchanger. In fact, it will be considerably smaller than the 
one that would have to be used with an open cycle. Probably 
the greatest difficulty will lie in the shielding arrangement. Itjg 
hoped that the geometry of the reactor shield will not offset the 
reduction in size of the plant which should be achieved by the 
thermodynamical choice of cycle and working medium. Beeause 
of the absence of combustion products, the life of the proposed 
plant may be longer. Finally, the plant will have fewer auxilig- 
ries, which should make its control and maintenance more trouble. 
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